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Introduction 



Finite size effects play an important role in modern statistical physics and quantum 
field theory. From a statistical point of view, it is known, that no phase transitions take 
place in a finite volume system, and the specific heat c(T) that is divergent at the critical 
point in infinite volume, becomes finite if the system has a finite size. Moreover, it is only 
in an interval around the critical temperature that finite size effects are relevant. Away 
from this interval, they are negligible because it is only near Tc that the correlation length 
is comparable with the size of the system. An important fact is that the specific heat 
and other critical quantities have a scaling behaviour as a response to varying the size L, 
which is fixed by the critical exponents of the infinite volume system [T], [2]. 

Interesting phenomena occur in quantum field theories as well. The most prominent 
example is the Casimir force between two neutral macroscopic bodies in vacuum [3]. This 
force can be determined from the volume dependence of the ground state energy, which is 
usually obtained by summing vacuum fluctuations. The crucial point is that besides the 
extensive bulk and boundary contributions (which are proportional to the volume of the 
domain and the surface of the boundaries, respectively) there are additional terms, which 
(in massive theories) decay exponentially with the distance of the boundary plates. In 
analogy with the Casimir energy, masses of stable particles also get exponentially small 
corrections in a finite volume [4] and the two-body interaction is also modified [5]. 

The knowledge of the properties of finite volume QFT is of central importance in at 
least two ways. On one hand, numerical approaches to QFT necessarily presume a finite 
volume box, and in order to interpret the results correctly a reliable theoretical control 
of finite size corrections is needed. Consider for example lattice QCD, where in order to 
approach the continuum limit the lattice spacing has to be adjusted as small as possible. 
The number of lattice points is limited by the computational resources, therefore it is 
typically not possible to choose a volume L^, which is large enough to neglect finite size 
effects. 

On the other hand, working in finite volume is not necessarily a disadvantage. On 
the contrary, the volume dependence of the spectrum can be exploited to obtain (infinite 
volume) physical quantities like the elastic scattering phase shifts [SI E] or resonance 
widths [TIE]. 
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In this work we investigate finite size effects in 1+1 dimensional integrable field the- 
ories. The main subjects of interest are finite volume form factors (matrix elements of 
local operators on eigenstates of the finite volume Hamiltonian) and expectation values 
and correlation functions at finite temperature. 

1+1 dimensional integrable models attracted considerable interest over the past thirty 
years. This interest is motivated by the fact that low dimensional integrable models rep- 
resent nontrivial interacting field theories that are exactly solvable: numerous physical 
quantities can be exactly calculated in these models. Phenomena qualitatively known in 
higher dimensions, such as universality, duality, etc. can be quantitatively investigated in 
1+1 dimensions. 

There is also a second motivation which is of a more practical nature. Several systems 
of condensed matter physics can be described by 1+1 dimensional effective field theories 
that in many cases lead to integrable models. Some important examples are the Kondo 
effect [9l[l0], and the effect of impurities in transport processes in general [TTl [T2l [T3l fT4]. 
behaviour of the edge states in the fractional quantum Hall effect, problems concerning 
spin chains, polymers, carbon nanotubes, etc. 

Besides the applications mentioned above two dimensional field theories also play a 
central role in string theory: they describe the dynamics of the string on the world sheet. 

In statistical physics, two dimensional critical systems can be described by conformal 
field theories (CFT) which are fixed points of the renormalization group fiow. In this 
sense, off-critical theories correspond to perturbations of CFT's by some of their relevant 
operators. Off-critical models do not possess the infinite dimensional symmetry algebra of 
the fixed point CFT, in particular scale invariance is lost. However, it was observed by A. 
B. Zamolodchikov [T5l fT6l fT7] that in certain situations there remains an infinite number 
of conserved quantities. In this case the resulting theory is still integrable and can be 
described by a factorized scattering theory. The integrals of motion restrict the possible 
bound state structure and mass ratios in the theory. Assuming further the bootstrap 
principle, ie. that all bound states belong to the same set of asymptotic particles, it is 
possible to construct the S-matrix with only a finite number of physical poles. 

Integrability can be exploited to gain information about the off-shell physics as well. 
Applying the ideas of analytic S-matrix theory and the bootstrap principle to form factors 
one obtains a rather restrictive set of equations which they have to obey [H]. These equa- 
tions can be considered as axioms for the form factor bootstrap, and supplied with the 
principles of maximum analyticity and the cluster property they contain enough informa- 
tion to determine the form factors completely [T8l I19| . Once the form factors are known, 
it is possible to construct correlation functions trough the spectral expansion, although 
the explicit summation of the series is only possible in some simple cases (20, EU [22j . 
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Integrability also offers powerful methods to explore the finite size properties of these 
models. It is possible to obtain the exact Casimir energy by means of the Thermodynamic 
Bethe Ansatz (TBA), both with periodic boundary conditions [23] or in the presence of 
non-trivial (integrable) boundary reflection factors [24]. A great deal of information is 
known about the excited state spectrum as well. Exact methods include the excited state 
TBA [25l[26] and nonlinear integral equations derived from lattice regularizations [271 [28] . 

However, less is known about finite size corrections to off-shell quantities. By Euclidean 
invariance correlation functions in a finite volume invariance correspond to the evaluation 
of thermal correlations. These objects can be compared directly to experiments, and it is 
not surprising that they attracted quite a lot of interest over the last decade [29l [30l [Ml [321, 
[331 [Ml [351 [361 [37l [381 [391 Ho] . Oiie possible way to approach finite temperature correlation 
functions is by establishing an appropriate spectral representation in finite volume. This 
constitutes the first motivation for our work, because finite volume form factors play a 
central role in this approach. Prior to our work only semiclassical results were known 
[4T1 [42] and also an exceptional exact result in the case of a free theory [43l [lH [45] . 

Besides being a promising tool to obtain correlation functions, finite volume form 
factors provide means to verify the bootstrap approach to form factors. The connection 
between the scattering theory and the Lagrangian (or perturbed CFT) formulation is 
rather indirect, however, it is generally believed, that the solutions of the bootstrap ax- 
ioms correspond to the local operators of the field theory. Evidence includes a direct 
comparison of the space of the form factor solutions with the space of primary operators 
of the CFT [46l [471 [48] ; the evaluation of correlation functions both through the spectral 
representation and from perturbation theory around the CFT [HI [49]; the evaluation of 
sum-rules like Zamolodchikov's c-theorem [50] and the A-theorem [5l], both of which can 
be used to express conformal data as spectral sums in terms of form factors. These tests 
concern quantities which are constructed from various integrals over the form factors, 
whereas the multi-particle matrix elements themselves have not been accessible. On the 
other hand, finite volume form factors can be obtained numerically in the perturbed CFT 
setting by means of the Truncated Conformal Space Approach (TCSA), and establishing 
the connection to infinite volume form factors provides a way to directly test the solutions 
of the form factor bootstrap. This constitutes the second motivation for our work. 

The first part of this thesis (chapters 3, 4 and 5) is devoted to the study of finite volume 
form factors in diagonal scattering theories. We show that these objects can be obtained 
in terms of the infinite volume form factors; in fact they are related by a proportionality 
factor which can be interpreted as a density of states of the finite volume spectrum. 

To our best knowledge, the first result on finite volume form factors in an interacting 
QFT was derived by Lellouch and Liischer in [52], where they considered the finite size 
dependence of kaon decay matrix elements. It was pointed out by Lin et.al. in [53] that 
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the final result of |52| can be explained simply by a non-trivial normalization of states in 
finite volume; they also derived a formula for arbitrary two-particle matrix elements. Our 
arguments are similar to the ones used in [53]. However, in integrable models it is possible 
to extend the calculations to arbitrary multi-particle states, because the factorized nature 
of the S-matrix allows of an analytic description of the whole finite size spectrum. 

Having established the connection to infinite volume form factors, we can directly 
access the form factor functions along certain one-dimensional sections of the rapidity- 
space parameterized by the volume L. The choice of the section corresponds to which 
multi-particle states we pick from the finite volume spectrum and it is only limited by 
the increasing numerical inaccuracy for higher lying states. This procedure provides a 
non-trivial check of our analytical calculations on finite volume form factors, and also a 
direct test of the bootstrap approach to form factors, as explained above. 

In the second part of the work (chapter 6) we develop a method to evaluate correlation 
functions at finite temperature: we introduce finite volume as a regulator of the otherwise 
ill-defined Boltzmann sum. We develop a systematic low-temperature expansion using 
finite volume form factors and show that the individual terms of this series have a well 
defined L ^ oo limit. In fact, they can be transformed into integral expressions over the 
infinite volume form factors. 

There have been previous attempts to attack the problem of finite temperature cor- 
relations. LeClair and Mussardo proposed an expansion for the one-point and two-point 
functions in terms of form factors dressed by appropriate occupation number factors con- 
taining the pseudo-energy function from the thermodynamical Bethe Ansatz [30]. It was 
argued by Saleur |31| that their proposal for the two-point function is incorrect; on the 
other hand, he gave a proof of the LeClair-Mussardo formula for one-point functions pro- 
vided the operator considered is the density of some local conserved charge. In view of the 
evidence it is now generally accepted that the conjecture made by LeClair and Mussardo 
for the one-point functions is correct; in contrast, the case of two-point functions (and 
also higher ones) is not yet fully understood. We contribute to this issue by comparing the 
low-temperature expansion obtained using finite volume form factors to the corresponding 
terms in the LeClair-Mussardo proposal. 

The thesis is organized as follows. In the next chapter we give a brief summary of 
factorized scattering theory and its connection to CFT, with a section devoted to the finite 
size properties of these models. Chapter 2 includes our analytic calculations concerning 
finite volume form factors, with two separate sections devoted to form factors including 
disconnected pieces. In chapter 3 we present our numerical studies: we test the analytic 
results of the previous chapter using TCSA in the Lee- Yang model and the Ising model 
in a magnetic field. In chapter 4 we study the leading exponential corrections (the /i- 
term) to scattering states and form factors. Chapter 5 deals with the evaluation of a 
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low-temperature expansion for correlation functions and also with the comparison of our 
results to previous proposals. The last chapter contains our conclusions. 

The material presented in this thesis is based on the papers 

• Balazs Pozsgay and Gabor Takacs: Form factors in fimte volume I: form fac- 
tor bootstrap and truncated conformal space Nucl. Phys. B788 (2007) 167-208, 
arxiv : 0706 . 1445 [hep-th] 

• Balazs Pozsgay and Gabor Takacs: Form factors in finite volume II: disconnected 
terms and finite temperature correlators Nucl. Phys. B788 (2007) 209-251, 
arxiv : 0706 . 3605 [hep-th] 

• Balazs Pozsgay: Liischer's mu-term and finite volume bootstrap principle for scatter- 
ing states and form factors Nucl. Phys. B802 (2008) 435-457, arxiv: 0803. 4445 [hep-th] 
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Chapter 1 



Integrable Models and Conformal Field 
Theories 

1.1 Exact S-matrix theories 

In this section we give a brief summary of the basic properties of factorized scattering 
theory. For an introduction to this field see the review article of Mussardo [54] . 

Let us consider a massive scattering theory with n particle types Aa with masses nia- 
One-particle states are denoted by |6')^, where 9 is the rapidity variable. We apply the 
following Lorentz-invariant normalization: 

a{ea\eb), = 2'K6ab6{ea-e^) 

Asymptotic states of the theory are defined as tensor products of one-particle states 
and are denoted by \6i, 62, ■ ■ ■ , Gn)a^a2 a„- The ordering of the rapidities is defined as 

• 9i > ■ ■ ■ > 9n ior in states 

• < ■ ■ ■ < for out states 

1.1.1 Conserved charges and factorization 

In integrable models there are an infinite number of conserved quantities, which can 
be represented on the quantum level as mutually commuting operators Qs, which act 
diagonally on the one-particle states: 

The quantity xl is the spin-s conserved charge of particle Aa. The spin represents the 
behavior under Lorentz-transformations. The charges for s = 1 and s = — 1 correspond 
to a linear combination of the energy and momentum, therefore = rria- 



7 



8 CHAPTER 1. INTEGRABLE MODELS AND CONFORM AL FIELD THEORIES 



The conserved charge operators are local in the sense, that they act additively on 
multi-particle asymptotic states: 

n 

Q,\e„ e,,..., = ( J2 xl) \0u 02,..., ^„),,.,.„.„ 

i=l 

The conservation of all higher spin charges constrains the S-matrix: 

• All scattering processes are elastic: the number of incoming and outgoing particles 
is the same. 

• The set of conserved charges (including the energy and the momentum) is the same 
for the in and out states. 

• The S-matrix factorizes: the amplitude for an arbitrary scattering process can be 
obtained by consecutive two-particle scattering processes. 

The basic object is therefore the two-particle S-matrix, which is defined for real ra- 
pidities by 

\0l,d2)ab = ^abi^^l - ^2)1^2, 

The factorization property forces the S-matrix to obey the Yang-Baxter equations [55l 

EEIEZ] 

where the summation over the repeated indices is understood. 

If there are no two particles with the same set of conserved charges, the S-matrix is 
diagonal: 

^abi^) = Sab{d)^ac^bd 

In this case the Yang-Baxter equations are automatically satisfied. In this work we only 
consider diagonal scattering theories. 

The S-matrix can be analytically continued to the whole complex plane and it satisfies 

Sab{e)s:,{e) = 1 (1.1) 

Sab{^n-9) = S^iiO) (1.2) 

due to unitarity and crossing symmetry, where Ai is the charge-conjugate of A},. 
The unitarity condition can be expressed as 

Sab{9)Sab{-9) = 1 (1.3) 

if one assumes real analyticity: 

Sab{e) = Sab{-e*) = s:,{-e) (i.4) 
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Equations (ll.ip . dl.Sp and (11.41) are rather restrictive and a general solution can be 
written as 

Sabie) = H t,{9) (1.5) 

with 

^ tanh i(6' + i7r7) sinh(6') + z sin(7r7) 

tanh ^{9 — inj) sinh(^) — i sin(7r7) 

and Qab C M is a finite set. 

1.1.2 Bound states - The bootstrap principle 

The S-matrix (11.51) may contain zeros and simple or multiple poles for purely imaginary 
rapidities. In a consistent theory, the interpretation for a simple pole is a formation of a 
bound state. If particle Ac appears to be a bound state of Aa and Ab, then one has 

Res,=,„c^ = ^{Kb? (1-6) 

with r;jj, G M and 

ml = ml + ml + 2mamb cos{ulb) 
The one-particle state 16*)^ may be identified with the formal two-particle state 

\e)^r^\e-iul,9 + tu:b),, (1.7) 

where the imaginary parts of the rapidities are determined by the conservation of energy 
and momentum under the fusion process (Fig. II. ll ): 



ma sin(M^J = mb sin«J m„ cos(m^J + cos{ulb) = 

As a consequence, one also has u^^b = + ^ab- For the higher spin conserved charges the 
identification (11.71) requires 

= X:e-^^^- + x^e^^^- (1.8) 

Scattering processes including Ac can be computed due to the factorization property of 
the S-matrix as 

Sdc{e) = s,a{e - iuic)Sdb{e + m^j (1.9) 

Equations (II. Sp and (II. 9p represent strong restrictions for the possible set of conserved 
charges and the sets Qab describing the scattering processes. The solution of these equations 
with a finite number of particles is called the bootstrap procedure. 



^Higher order poles in the S-matrix are specific to two dimensions and can be explained by the 
Coleman-Thun mechanism [MIIMIEQI. 
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(a) Fusion angles (b) The mass triangle 

Figure 1.1: Pictorial representation of particle fusions. 

1.2 Integrable theories in finite volume 

Let us consider the theory defined in a finite box of size L. One of the most character- 
istic feature of a field theory is the evolution of its discrete energy levels -Ej(L), i = 1 . . . oc 
as a function of the volume. 

In quantum field theories (including realistic 3+1 dimensional models) low-lying one- 
particle and two-particle states [HlHllT] can be described to all orders in 1/L in terms of the 
infinite volume spectrum (particle types and masses) and the elastic scattering amplitudes 
of the theory. In non-integrable theories higher lying energy levels are far more difficult 
to access due to the complicated structure of the S-matrix and only partial results are 
available [Ml [62l [63j . On the other hand, the S-matrix of integrable models factorizes and 
there are no inelastic processes present. These properties enable us to describe the whole 
finite volume spectrum to all orders in 1/L. 

In this section we briefiy describe the main characteristics of the finite volume spectra 
of diagonal scattering theories. We also introduce some notations which will be used 
throughout this work. We assume periodic boundary conditions for simplicity. 

1.2.1 One-particle states 

One-particle states of particle Aa in finite volume can be characterized with the mo- 
mentum Pa, which (due to the periodicity of the wave function) is constrained to 

Therefore, the particle type and the momentum quantum number / (the Lorentz-spin) 
uniquely determine the one-particle state in question. For convenience we introduce the 
notation 

\{I})a,L 

where the subscript L denotes, that this state is an eigenvector of the finite volume 
Hamiltonian. 
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Figure 1.2: The diagram to the left (the /x-term) shows a particle splitting in two virtual, 
on-shell particles, traveling around the cylinder and recombining. The diagram to the right 
(the F-term) shows a virtual particle going around the circumference of the cylinder. 

The energy is given to all orders in 1/L by the relativistic formula 

^=^^^^1+^ (1.10) 

There are additional correction terms to ( ll.lOp . which decay exponentially with the 
volume. These terms can be attributed to virtual processes which are absent in infinite 
volume. M. Liischer calculated the finite size mass corrections of stable particles (energy 
corrections for 1 = 0); they consist of the so-called /x-term and F-term [4] and are associ- 
ated to diagrams depicted in figure [L2l The explicit formulas for one dimensional diagonal 
scattering theories read [64] 

Ami'^) = -J^^^i^l-l^l-^lDMabcf^atciK.fe-^^^^'^ (1.11) 

b,c 

/oo jn 
-m,LcoshW^^^^g^^^) (^a,(^ + m/2) - 1) (1.12) 

where A4abc = 1 if is a bound state of Aa and Ab and zero otherwise; fiabc is the 
altitude of the mass triangle with base (see figure II. ip and (F^^)^ is the residue of 
Sab{0) corresponding to the formation of the bound state. 

1.2.2 Many-particle states 

The Bethe-Yang equations [65] serve as quantization conditions for a many-particle 
state consisting of n particles Ai^ with rapidities 9^: 

Qk{Oi, . . . ,9n)i^..A„ = m,,Lsmh9k + '^6i^i^{9k - Oi) = 271^ , k = l,...,n (1.13) 
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where Sij{9) = —i\ogSij{9). 

For each set of particle types and momentum quantum numbers Ik there is a unique 
solution of the equations above|^ The energy of the many-particle state in question is 
given to all orders in 1/L by the additive formula 

n 

E = J2^hCosh{9k) (1.14) 

k=l 

The effect of the interaction between the particles is the shift of the rapidities from the 
values they would assume in a free theory. 

Since many-particle states are completely characterized by the quantum numbers Ik, it 
is convenient to introduce the following implicit notation for a n-particle state consisting 
of particles ii,i2, ■ ■ ■ ,in with rapidities 61,62, ... ,6n- 

|{/l,/2, . . • ,-^n})i,i2,...i„,L 

where it is understood that 61,62, ... ,6n are solutions of the Bethe-equations with quan- 
tum numbers Ji, /2, . . . , /„. 

To have an unambiguous definition of the quantum numbers Ik, it is convenient to 
define phase shift functions 6ab which are continuous and odd functions of the rapidity 
difference 6; we achieve this using the following convention: 

where Sab is uniquely specified by continuity and the branch choice 

SabiO) = 

and it is an odd function of 6 due to the following property of the scattering amplitude: 

Sab{6)Sabi-6) = 1 

(which also implies Sab{0) = ±1). For a generic complex rapidity one has 

5abi6*) = Sabi6y 

Note, that with the convention above the quantum numbers Ik can be both integer and 
half-integer. Two particle states of the type AaAa have generally half-integer quantum 
numbers, because the particles have fermionic statistics (S'aa(O) = —1) in all known the- 
oriej^ . 

Similar to the one-particle energies, there are exponentially small finite size corrections 
to p.l4p . The /i-term associated to multi-particle energies is investigated in section [421 

^The solutions with two or more coinciding momenta for a given particle type do not correspond to 
eigenstates of the Hamiltonian due to Pauli principle. 
^The only exception is the free boson with S — I 
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1.3 Scattering theories as perturbed CFT's 

Conformal Field Theories are two-dimensional Euclidean field theories, which possess 
invariance under conformal transformations, including scale-invariance. In the space of 
all possible QFT's they represent fixed points under the renormalization group flow. In 
statistical physics, they describe fluctuations of critical systems in the continuum limit 
[66] ■ Conformal invariance highly constrains the behavior of the correlation functions, 
and even the operator content of the theory p7]. In this section we quote some basic 
properties of CFT's and introduce two particular models, which will serve as a testing 
ground for our results: the scaling Lee- Yang model and the critical Ising model. For a 
general introduction to conformal held theories see [68l [69] . 

Scale invariance constrains the energy-momentum tensor T^i^{x,y) to be traceless. 
Introducing complex coordinates z = x + iy and z = x — iy this condition can be expressed 
as Tzz = Tgz = 0. Conservation of the energy-momentum tensor on the other hand requires 
dzTzz = dzTzz = 0, therefore it is possible to introduce the left- and right-moving (chiral 
and anti-chiral) components T{z) = Tzz{z) and T{z) = Tzz{z). 

T{z) may be expanded into its Laurent-series around 2; = as 

T{z) = Y,z—^L^ 

neN 

The operators L„ satisfy the Virasoro algebra: 

[Ln, Lm] = {n- m)Ln+m ^ ^^'^^ ~ ^)^n+m,0 

which is a central extension of the symmetry algebra of the classical conformal group, c 
being the central charge of the theory. 

The value of c includes a great deal of information about the theory. It restricts the 
possible representations of the Virasoro-algebra, therefore also the operator content and 
the spectrum of the theory. The simplest theories are the minimal models [70], which 
contain a finite number of primary fields and possess no additional symmetries. They can 
be characterized by two coprime integers p and q and the central charge is given by 

Q{p-q)^ 

c = 1 

pq 

These theories are in addition unitary (there are no negative-norm states in the spectrum) 
if q = p + 1. 

1.3.1 CFT in a cylindrical geometry 

A cylinder with spatial circumference L can be mapped to the complex plane by the 
conformal transformation 

z = exp ^— (r — ix) j z = exp (^—{t + ix)^ (1-15) 
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where x and r are the spatial and the imaginary time coordinates, respectively. The trans- 
formation properties of the energy-momentum tensor determine the Hamilton-operator 
(generator of translations in the time-direction): 

H= f dxTrr{x,T = 0) = '^(^Lo + lo-^) (1.16) 

In minimal models the Hilbert-space is given by 

n = ^Vh^Vh (1.17) 

h 

where Vh (Vh) denotes the irreducible representation of the left (right) Virasoro algebra 
with highest weight h. 

1.3.2 Perturbing CFT's 

Conformal field theories represent statistical physical or quantum systems at criticality. 
However, they can be also used to approach noncritical models. 
Let us consider a theory defined by the action 



A = Ac FT + ^9j j d^x (x) 



where Acft is the action of some CFT and $j are its relevant operatorcl. This theory 
is no longer conformal, the coupling constants gj describe the deviation from criticality. 
However, if there is only one perturbation present, which only brakes a subset of the 
conformal symmetries, the theory may still possess an infinite number of conservation 
laws, and it may remain integrabl^ [ISl [EJ [III C2]- Let us therefore consider that the 
theory defined by the action 

A = Acft + g j d^x $(x) (1.18) 

is integrable. Scale-invariance is broken by the perturbation, the energy scale in the per- 

1 

turbed theory is set by (7 2-^* , where A$ is the scaling dimension of the field Depending 
on the original CFT and the perturbing operator the action p. 181) may define a massive 



In the sense of renormalization group flow, relevant operators describe the departure from the critical 
system. The perturbation leaves the UV behaviour of the theory unchanged, ie. it is still governed by 
the CFT. If the scaling dimension of the operator is greater or equal to 1, a renormalization procedure 
is required to properly define the perturbed theory[7l] On the other hand, perturbing with irrelevant 
operators results in an ill-defined theory in the UV, with an infinite number of counter-terms occurring 
during renormalization }72j . 

^Perturbations with multiple couplings are not expected to lead to integrable theories. For a short 
discussion see section 2.3 in 154). 
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[74] or a massless [75l [76] scattering theory. In this work we only consider the massive 
case. 

The question of which perturbed CFT corresponds to which scattering theory may 
be answered by different methods. First of all, Zamolodchikov's counting argument [I5] 
provides a sufficient condition for the existence of higher spin conservation laws. A second 
independent method is the Thermodynamic Bethe Ansatz [23l [77l [78], which provides 
integral equations to determine the finite size dependence of the vacuum energy in terms 
of the S-matrix of the theory. It can be used to analytically predict the central charge of 
the CFT, the non-universal bulk vacuum energy density, and also the scaling dimension 
of the perturbing field. 

A third method is the Truncated Conformal Space Approach (see 11.3.51) . which can 
be used to numerically determine the low-lying energy levels of the finite size spectrum. 
Matching the numerical data with the predictions of the Bethe- Yang equations fll.lSp one 
can identify multi-particle states and thus directly test the phase shifts Sab{,0)- For the 
details of this procedure see section 13.1.51 

1.3.3 The Lee- Yang model 

The non-unitary minimal model M2,5 has central charge c = —22/5 and a unique 
nontrivial primary field $ with scaling weights A = A = — 1/5. The field $ is normalized 
so that it has the following operator product expansion: 

^{z, z)^{0, 0) = C(2z)^/^$(0, 0) + {zzf^ + ... (1.19) 

where I is the identity operator and the only nontrivial structure constant is 

C = 1.911312699- •• X i 

The Hilbert space of the conformal model is given by 

nLY= Vh^Vk 

h=0,-l/5 

where Vh (V^) denotes the irreducible representation of the left (right) Virasoro algebra 
with highest weight h. 

The off-critical Lee- Yang model is defined by the Hamiltonian 

H^^^ = HcFT + iX [ «(0,x) (1.20) 
Jo 

where Hcft is the conformal Hamiltonian. This theory is related to the Lee- Yang edge 
singularity of the Ising model in an imaginary magnetic field [79l [80l [81] . Despite its lack 
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of unitarity, it is also relevant in condensed matter physics: it is related to the theory of 
non-intersecting branched polymers in 4 dimensions |82| . 

When A > the theory above has a single particle in its spectrum with mass m that 
can be related to the coupling constant as [23] 

A = 0.09704845636 ■ • • x m^^/^ (1.21) 

and the bulk energy density is given by 

B = -^m'' (1.22) 

12 ^ ^ 

The S'-matrix reads [83] 

^ ^ sinh 9 + isin^ , ^ 

Sly{0) = . . . I (1.23) 

smh (f — t sm — 

and the particle occurs as a bound state of itself at = 2ni/3 with the three-particle 
coupling given by 

= -2^3 

where the negative sign is due to the non-unitarity of the model. According to the con- 
vention introduced in section 11.21 we define the phase-shift via the relation 

SLvie) = -e'^(') (1.24) 

so that 5(0) = 0. 



1.3.4 The Ising-model in a magnetic field 

The critical Ising model is the described by the unitary conformal field theory M3 4 with 
c = 1/2 and has two nontrivial primary fields: the spin operator a with = = 1/16 
and the energy density e with = = 1/2. The magnetic perturbation 

H = Hl + h [ dxa{0,x) 
Jo 

is massive (and its physics does not depend on the sign of the external magnetic field h). 
The spectrum and the exact S matrix is described by the famous Eq factorized scattering 
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theory [84], which contains eight particles Ai, i = 1, . . . ,8 with mass ratios given by 

vr 

1712 = 2mi cos — 
vr 

mcj = 2mi cos — 
30 

771 

rriA = zmo cos — 
30 

27r 

ms = 2m2 cos — 
15 

TT 

me = 2m2 cos — 
30 
vr 

my = 2m4 cos — 
5 

TT 

ms = 2m5 cos — 
5 

and the mass gap relation is [85] 

mi = (4.40490857. . 

or 

h = K/,mf , Kh = 0.06203236 . . . (1.25) 
The bulk energy density is given by 

5 = -0.06172858982 •■■ X (1.26) 
We also quote the scattering phase shift of two Ai particles: 

All other amplitudes Sab are determined by the S matrix bootstrap [84]; the only one we 
need later is that of the Ai — A2 scattering, which takes the form 

5}, {15}, {5}, {15 

Our interest in the Ising model is motivated by the fact that this is the simplest model 
in which form factors of an operator different from the perturbing one are known, and also 
its spectrum and bootstrap structure is rather complex, both of which stands in contrast 
with the much simpler case of scaling Lee- Yang model. 

1.3.5 Finite size spectrum from Conformal Field Theory - the 
Truncated Conformal Space approach 

Let us assume that the infinite volume scattering theory corresponds to a perturbed 
conformal theory defined by the action 

A = AcFT + 9 (fx $(x) (1.28) 



Si2{e) 



18 CHAPTER 1. INTEGRABLE MODELS AND CONFORM AL FIELD THEORIES 



If this connection is already established, one can use the conformal data to obtain the 
finite size spectrum of the theory. The perturbed Hamiltonian is obtained by mapping the 
cylinder with circumference L to the complex plane. Using the transformation properties 
of the primary field one has 



where Hcft is the conformal Hamiltonian (II. 16^ . Both the spectrum Hcft and the matrix 
elements of $ between its eigenstates can be calculated, therefore the complete and exact 
finite size spectrum can be obtained in principle by diagonalizing ( 11.29^ in the infinite 
dimensional Hilbert-space (ll.lTp . 

Instead of the impossible task of diagonalizing an infinite matrix Yurov and Zamolod- 
chikov proposed the Truncated Conformal Space Approach (TCSA) [86] as a numerical 
approximation scheme. The method consists of truncating the Hilbert-space (ll.lTp to the 
states which have scaling dimension (eigenvalue under Lq + Lq) smaller than some thresh- 
old Ccut- In this finite Hilbert-space the diagonalization procedure can be performed on a 
computer, and one obtains eigenvalues Ei{L, ecut), which are approximations to the exact 
energy levels Ei{L), and 



The convergence is faster for the low-lying states. 

Besides being an effective tool to obtain finite size quantities, TCSA can also serve as 
a regulator of the UV divergences of the vacuum energy if A$ > 1. 

In this work (chapter [S]) we use TCSA to determine the finite volume form factors of 
local operators in the scaling Lee- Yang and Ising models. The numerical data is then used 
to confirm the analytic results of chapter 2. The details of the TCSA methods we used 
are explained in section 13. 1[ 




(1.29) 



lim Ei{L,ecut) = Ei{L) 



■oo 



Chapter 2 

Finite Volume Form Factors - Analytic 
results 

In this chapter we determine the finite volume matrix elements in terms of the infinite 
volume form factors and the S-matrix of the theory. 

In section 12.11 we briefiy summarize the ingredients of the so-called form factor boot- 
strap program, which leads to explicit solutions for the infinite volume form factor func- 
tions. In section 12.21 we present our results for finite volume form factors which do not 
include disconnected pieces. 

Disconnected terms occur in the case of diagonal form factors and matrix elements 
including zero-momentum particles. These two special situations are discussed in sections 
12.31 and [2?4l respectively. 



2.1 The form factor bootstrap program 

The (infinite volume) form factors of a local operator 0{t, x) are defined as 

-^mn(^l' • • • 5^ml^l' • • • > ^n)ji...j„;n...j„ = ji...>„(^l, • • • > ^'^ I C'lO, 0) 1 6*1, . . .,6n)h...i„ (2.1) 

In first place they are defined for real rapidities, however they can be analytically continued 
to the whole complex plane. The properties of the form factor functions are explained in 
detail in Smirnov's review |18| . 

With the help of the crossing relations 

-^mnX^l' • • • ' ^ml^l' • • • ' ^n) ji...jm;ii-..in ~ 

-^m-ln+l(^l' • • • ' ^m-ll^m + ^1; • • ■ ; dn) ji...jm-i;jmil-..in 

n fe— 1 

k=l 1=1 
-^m-ln-l(^l) • • • ) ^■m-l\^li ■ ■ ■ t ^k~l, • • • , dn) ji...jm-i;jmh...ik-iik+i---in ) (2-2) 
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all form factors can be expressed in terms of the elementary form factors 

F^ie^, . . . , = (o|0(o, 0)1^^1, . . . , 



which satisfy the following axioms [87| : 
I. Exchange: 



-F„ (^1, . . . , 6k, Ok+l, . . . , On)ii...ikik+i---in ~ 

Sikik+ii^k - 9k+l)F^{9i, . . . , 9k+l, Ok,..., 9n)h...ik+iik-ir. (2-3) 

II. Cyclic permutation: 

F^{e, + 2m, 92,..., = Ff (^2, . . . , ^n, ^i) (2.4) 

III. Kinematical singularity 

- i Res F^^^{9 + it:, 9', 9^,..., ^„),,.,....„ =[\~ 6,, f[ S,,,{9 - 9^)^ F^{9,, 9^),,...,„ 

(2.5) 



k=l 



IV. Dynamical singularity 

- z Res F^_,2iO + iu)j2, 9' - %u{j2, 9^,..., e„)., = V%F^^^{9, 9,,..., 9n)kn...^ (2.6) 
e=e' 

whenever k occurs as the bound state of the particles i and j. 

The essence of the form factor bootstrap program is to obtain explicit analytical solu- 
tion to the above set of recursive equations [87l EHl [HI [89] . Axioms I-IV are supplemented 
by the assumption of maximum analyticity (i.e. that the form factors are meromorphic 
functions which only have the singularities prescribed by the axioms) and possible further 
conditions expressing properties of the particular operator whose form factors are sought. 

2.2 Finite volume form factors without disconnected 
pieces 

2.2.1 Elementary form factors 

Here we determine finite volume form factors of the form 

{Q\0\{h,h,...,In})n...i^,L 

by comparing the Euclidean two-point functions of arbitrary local operators 

(C»(x)C>'(0,0)) and (C(x)O'(O, 0))l (2.7) 
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defined in the infinite volume theory, and in a finite but large volume L, respectively. 
Given that one uses the same renormalization prescriptions for the operators, the finite 
size correction to the correlation function decays exponentially with L: 

{O{x)O'{0)) - {Oix)O'{0))L ~ 0(e~^^) (2.8) 

where is some characteristic mass scale. In 12.2.31 we discuss this relation in detail; here 
we only use the fact, that the finite size correction decays faster than any power of 1/L. 

Choosing x = (r, 0) and inserting a complete set of states one obtains the spectral 
representations 

n=0 h...i„ \A:=1 / 

Ff (01, 02,..., 0„)+ exp (^-T cosh 9^ (2.9) 

and 

oo 

(O(r,0)O'(0,0))i = E(0|^(0'0)l{^i'^2,...,/n})n....,LX (2.10) 

n=0 i\...in h-.-In 

i,..,A{h, h,..., In}\O'i0, 0)\0)l exp I -T J2 ^iu cosh^, 



k=l 



where we used multi-particle states and form factors of the infinite volume and finite 
volume theories, respectively, and in (I2.10p it is understood that the rapidities 6k are 
solutions of the Bethe-Yang equations 

Qfe(^i, . . . ,i9„)j,...j„ = mj^Lsinhi9fc + ^(5i^i,(6'fc - fii;) = 27r4 , A; = l,...,n (2.11) 
The second form factor in (12. 9p is defined as 

Ff (01, 02, ... , ^n) + ....„ = 0n\O'{0, 0)|0) = Ff (^i + ZTT, 9^ + 171, 0„+^7r),,...,„ 

, which is just the complex conjugate of F^' for unitary theories. 

To relate the finite and infinite volume form factors a further step is necessary, because 
the integrals in the spectral representation (12. 9p must also be discretized. Let us consider 
this problem first for the case of free particles: 

where 

Pk = rrii^^ sinh 0k , Ek = rrii^^ cosh 0k 
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are the momenta and energies of the particles. In finite volume 

2lTlk 

and it is well-known (as a consequence of the Poisson summation formula, cf. [5]) that 

5: , . . . , ^) = (^) Yn r dp^ 9{p,, . . . + o{L--) (2.12) 

provided the function g and its first derivatives are integrable. Recalling that form 
factors are analytic functions for real momenta, in our case this is true for derivatives of 
any order, due to the exponential suppression factor in the spectral integrals, provided 
the form factors grow at most polynomially in the momentum, i.e. 

|i^n(^i + ^,^2 + ^,...,^n + ^)| ~e"l'l as 1^1 ^oo 

This is true if we only consider operators which have a power-like short distance singularity 
in their two-point functions [89j : 

(0|O(x)O(0)|0) = ^ 

Therefore the discrete sum differs from the continuum integral only by terms decaying 
faster than any power in i.e. by terms exponentially suppressed in L. Taking into 
account that (12.71) is valid for any pair O, O' of scaling fields, we obtain 

(O|0(O, 0)|{/i, . . . , = , ^ \ =F^{Qi. • . . , ^~n)n....„ + 0(e-^'^) 



(2.13) 



where 



sinh Qi, ^ 



and 

n 

pfl^S^x, . . . ,^~„) = n^H.^cosh^~fc (2.14) 

fc=i 

is nothing else than the Jacobi determinant corresponding to changing from the 
variables l-nlk to the rapidities Qk. The term 0{q~^'^) signifies that our considerations 
are valid to all orders in 1/L, although our argument does not tell us the value of /x': to 
do that, we would need more information about the correction term in the discretization 

In the case of interacting particles a more careful analysis is necessary. To generalize 
(I2.13P one has to employ the Jacobian of the mapping {Ji, /2, • • • , /n} — ^ {^i, 6*2, • • • , Qn\ 
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described by the Bethe-Yang equations. Therefore we define 
p,,...,„(ei,...,e„) = det:r(")(^i,---,^nk...i„ (2.15) 

j{n)in ay . _ dQkjOl, ■ ■ ■,On)h...i„ kl-^ n 

J^l [Ol, . . . ,On)i^...t„ — , K, < — i,...,n 

The finite volume form factors are then given by 

(0|O(0, 0)|{/i, . . . , In})n...^,L = ] ^ =F^iOi, en)n...„ + 0(e-^'^) 

y Ph...in{(^l, • • • ,^n) 

(2.16) 

where 9k are the solutions of the Bethe-Yang equations (12.111) corresponding to the state 
with the specified quantum numbers Ji,...,/^ at the given volume L. Note, that the 
Bethe-Yang equations are correct to all orders in 1/L, therefore there can be no error 
terms in (12.16p analytic in 1/L. In subsection 12.2.31 we consider the error exponent n' in 
detail. 

The quantity Pi^...iS^i, ■ ■ ■ ,9n) is nothing else than the density of states in rapidity 
space. It is also worthwhile to mention that relation (12. 16^ can be interpreted as an 
expression for the finite volume multi-particle state in terms of the corresponding infinite 
volume state as follows 

. . . , In})ii...i„,L = I ^ • • • , On)h...i„ (2-17) 

YPn...j„(6'i, ■■■,0n) 

This relation between the density and the normalization of states is a straightforward 
application of the ideas put forward by Saleur in [31] . 

To get an idea of the structure of Pii...j„(6'i, . . . ,6n) it is instructive to consider the 
simplest cases. The quantization rule for the one-particle states is identical in the free and 
interacting cases, therefore 

Pi{6) = EiL with Ei = mjCosh(6') 

In the two-particle case one has 



EiL + ipi^i^ (6*1 - 6*2) {di - O2) 

-Vhi2 (^1 - ^2) E2L + (6*1 - 62] 



where we used the symmetric function ipi^i^{6) = d/dO 5i^i^{6) and Ei^2 = "^ii,2 cosh(6'i^2)- 
The two-particle density is therefore given by 

PniM.d2) = E^E2L^ + (El + E2)L ^,,,,{6^ - 62) 

The leading term in the expression above is simply the product of two one-particle den- 
sities, just like in the free case. This is a general rule: the leading 0{L^) term in the 
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n-particle density is given by the product of one-particle densities. The effect of the inter- 
action is the appearance of sub-leading terms proportional to ipi^i^{6). These terms turn 
out to be crucial in the calculations of thermal correlations (see section 15. 2p and in direct 
numerical tests of the form factors. The necessity of introducing the full densities (12.151) 
was first shown in [H] and it is demonstrated in great detail in section [3l 

Note, that there is no preferred way to order the rapidities on the circle, since there 
are no genuine asymptotic in/ out particle configurations. This means that in relation 
( 12.16P there is no preferred way to order the rapidities inside the infinite volume form 
factor function F^ . Different orderings are related by S'-matrix factors according to the 
exchange axiom (12.31) . which are indeed phases. Such phases do not contribute to corre- 
lation functions (cf. the spectral representation (12.91) ). nor to any physically meaningful 
quantity derived from them. In section [3^ we show that relations like ( 12. 16^ must always 
be understood to hold only up to physically irrelevant phase factors. 

We also remark, that there are no finite volume states for which the quantum numbers 
of any two of the particles are identical. The reason is that 

^.^(0) = -1 

(with the exception of free bosonic theories) and so the wave function corresponding to 
the appropriate solution of the Bethe-Yang equations (12.111) vanishes. We can express this 
in terms of form factors as follows: 

(0|O(0,0)|{/i,/2,...,/„}).,...„,L = 

whenever Ik = h and ik = ii for some k and /. Using this convention we can assume that 
the summation in ( 12.101) runs over all possible values of the quantum numbers without 
exclusions. Note that even in this case the relation ( I2.16p can be maintained since due to 
the exchange axiom (12. Sp 

F^{9i, . . . , 9n)i^...ir, = 
whenever 6k = 9i and ik = k for some k and I. 



2.2.2 Generic form factors without disconnected pieces 

Using the crossing formula (12. 2p . eqn. (I2.17P allows us to construct the general form 
factor functions (12. ip in finite volume as follows: 

,,..,J{I[, IL}\0{0, 0)\{h, In})n...^,L = 

FZ+ni9'rn + ZTT, . . . , + ZTT, , . . . , On) ^ 0(6"'^^) (2 18) 



\l Ph...iS9u ■ ■ ■ ,9n)Pn...jmi9l, ■■■,9'^) 
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provided that there are no rapidities that are common between the left and the right states 
i.e. the sets jfi'i, . . . ,9n^ and |i9^, . . . , i9^| are disjoint. The latter condition is necessary 
to eliminate disconnected pieces. The rapidities entering (12.180 are determined by the 
Bethe-Yang equations f ll.lSp : due to the presence of the scattering terms containing the 
phase shift functions S, equality of two quantum numbers Ik and // does not mean that 
the two rapidities themselves are equal in finite volume L. It is easy to see that there are 
only two cases when exact equality of some rapidities can occur: 

1. The two states are identical, i.e. n = m and 

{jl---jm} = {^l...^„} 

in which case all the rapidities are pairwise equal. The discussion of such matrix 
elements is presented in section 12.31 

2. Both states are parity symmetric states in the spin zero sector, i.e. 

= {— In, . . . , — Ji} 

and the particle species labels are also compatible with the symmetry, i.e. in+i-k = ik 
and jm+i-k = jk- Furthermore, both states must contain one (or possibly more, in 
a theory with more than one species) particle of quantum number 0, whose rapidity 
is then exactly for any value of the volume L due to the symmetric assignment of 
quantum numbers. The discussion of such matrix elements is presented in section 

m 

We stress that eqns. (12.161 12.18P are exact to all orders of powers in 1/L; we refer to 
the corrections non-analytic in 1/L (eventually decaying exponentially as indicated) as 
residual finite size effects, following the terminology introduced in [8]. 

We would like to remark, that similar arguments that led us to (12.181) were previously 
used to obtain the finite size dependence of kaon decay matrix elements by Lin et al. [53] . 
The idea of normalizing finite volume form factors with particle densities also appeared 
in [90]; however, they used the unsatisfactory free-theory densities. 

2.2.3 Estimation of the residual finite size corrections 

There are two sources of error terms which may contribute to (12.161) : 
• Finite size corrections to the correlator (C(r, 0)O'(0, 0))l. 
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• Discretization errors introduced in (12 .121) . 

Here we show, that both types of error terms behave as 0{e~^'") where is some charac- 
teristic mass scale of the theory. 

One can use Liischer's finite volume expansion introduced in ^ to study finite size 
effects to the correlator. According to Liischer's classification of finite volume Feynman 
graphs, the difference between the finite and infinite volume correlation function is given 
by contributions from graphs of nontrivial gauge class, i.e. graphs in which some propa- 
gator has a nonzero winding number around the cylinder. Such graphs always carry an 
exponential suppression factor in L, whose exponent can be determined by analyzing the 
singularities of the propagators and vertex functions entering the expressions. In a massive 
theory, all such singularities lie away from the real axis of the Mandelstam variables, and 
the one with the smallest imaginary part determines /x. It turns out that the value of /i is 
determined by the exact mass spectrum of the particles and also the bound state fusions 
between them [5],[64]- Therefore it is universal, which means that it is independent of the 
correlation function considered. In general jj, < m where m is the lightest particle mass 
(the mass gap of the theory), because there are always corrections in which the lightest 
particle loops around the finite volume L, and so the mass shell pole of the corresponding 
exact propagator is always present. Contributions from such particle loops to the vacuum 
expectation value are evaluated in subsection 13.21 (for a graphical representation see figure 
13. 3p , while an example of a finite volume correction corresponding to a bound state fusion 
(a so-called /x-term) is discussed in subsection 13.2.21 (figure IXSl) . The /i-terms associated 
to scattering states and form factors are analyzed in chapter [H 

The estimation 

(0(r, O)0'(O)) - (0(r, 0)(9'(0))i ~ 0(6"^^) (2.19) 

is in fact a pessimistic one. Making use of the Euclidean invariance of the theory and 
performing the rotation (r, x) — > (— x,r), the finite-size effects in the relation above can 
be obtained by calculating finite temperature corrections to the correlations function with 
the temperature given by T = 1/L. In section [531 we consider such contributions; they 
are of order e""^^, ie. /i = m with m being the mass gap of the theory. The absence of 
correction terms with /i < m can also be explained by the fact, that it is impossible to 
draw a relevant finite volume diagram (possibly including a particle fusion) which would 
carry the suppression factor e~^^. 

Having discussed finite size effects to the the correlator, we now prove that the dis- 
cretization procedure introduces error terms which behave as e~^'^ , where in fact /x' = /i. 

Recall that the Poisson formula gives the discrete sum in terms of a Fourier transform: 
the leading term is the Fourier transform of the summand evaluated at wave number 
(i.e. the integral) and the corrections are determined by the decay of the Fourier transform 
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at large wave numbers. The function we need to consider is 

h{Pl^ ■■■■,Pn) = ^ Pii...i„{Ol, • • • , dn)~^F^{Oi, 6^2, ... , 9n)i^...ir,Fn (^1' 6^2, . . . , 6^, 



«l...i„ 



n 



exp I — r mj^^ cosh ^^jt ) (2.20) 

k=l 

where pk = sinh 9k are the momentum variables. Due to the analyticity of the form 
factors for real rapidities, this function is analytic for physical (real) momenta, and to- 
gether with all of its derivatives decays more rapidly than any power at infinity. Therefore 
its Fourier transform taken in the momentum variables has the same asymptotic property, 
i.e. it (and its derivatives) decay more rapidly than any power: 



/i(k1, . . . , K„) 



for large k. As a result, discretization introduces an error of order e~^'^. The asymptotic 
exponent /i' of the Fourier transform can be generally determined by shifting the contour 
of the integral transform and is given by the position of the singularity closest to the real 
momentum domain (this is essentially the procedure that Liischer uses in [4]). Singularities 
of the form factors are given by the same analytic structure as that of the amplitudes 
which determine the exponent in eqn. (12.81) . Thus we find that fi' = /i. 

This argument is just an intuitive reasoning, although it can be made a little more 
precise. First of all, we must examine whether the determinant Pji...i„(^i, . . . , On) can have 
any zeros. It can always be written in the form 



Pi 



Xe,,...,en) = I l[m,^L cosheA {I + o{i/L)) 



\k=l 



The leading factor can only be zero when 0^ = ^ some k, which corresponds to 
Pk = iruk, giving fi' = ruk in case this is the closest singularity. That gives a correction 

which is the same as the contribution by an on-shell propagator wound around the finite 
volume, and such corrections are already included in the term of (12. Sp . Another possi- 
bility is that some phase-shift function 5{6) in the 0(1/L) terms contributes a large term, 
which balances the 1/L pre-factor. For that its argument must be close to a singularity, 
and then according to eqn. ( 11. 6p we can write 

5{e) ~ \og{e - m) ~ o{L) 



where u is the position of the singularity in the phase-shift. This requires that the singu- 
larity is approached exponentially close (as a function of the volume L), but the positions 
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of all these singularities are again determined by singularities of the vertex functions, so 
this gives no new possibilities for the exponent /i'. 

A further issue that can be easily checked is whether the Fourier integral is convergent 
for large momenta; the function (12 .20^ is cut off at the infinities by the factor 

n 

Y\ exp(— mjj.r cosh 6'^) 

k=l 

which can only go wrong if for some k 

3fJe cosh 9k <0 

but that requires 

'^mOf. > — 
fe 2 

which is already farther from the real momentum domain then the position of the on-shell 
propagator singularity. 

2.2.4 Comparison with an exact result 

As far as we know, the only exact result in the literature was obtained in [45| and 
independently in [43l[44], where the authors considered the finite volume form factors of 
the spin-field operator a in the Ising-model at zero external field. Although the model in 
question is a free theory (it can be described in terms of free massive Majorana fermions), 
it is interesting to compare the results with the general rule presented in the previous 
subsections. On one hand, this confirms our calculations. On the other hand, it might be 
used in further work as a starting point to study sub-leading exponential corrections (for 
the leading term, see section [4731 ). 

The Hilbert space of the theory consists of two sectors: the Neveu-Schwarz (NS) and 
the Ramond (R) sectors. The fermionic fields are antiperiodic or periodic, respectively, 
leading to the quantization rule for the particle momenta 

p = msmh.{6) = 
p = msinh(^) = 
Multi-particle states are denoted by 

. . . , k]\f)-^^ 

for states of the NS and R sectors, respectively. For simplicity we only consider the T > T^ 
phase of the theory, where spin-fiip symmetry is unbroken. In this phase the allowed 
excitations have an even number of particles in both the NS and R sectors. 



L^ 

— n n G N 



ken + - (NS) 



(R) 



(2.21) 
('2.221 



and 



|ni, . . . ,n7v) 



R 
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The spin operator has non- vanishing matrix elements only between states of different 
sectors. The result for a generic finite volume matrix element reads (eq. 2.12 in [45]) 



NS 



{h, h 



kK\cr{0,0)\ni,n2, - ■ ■ ,nN)R = 

K N 

S{L)\[~g{9,^)\[g{er,^) Fk. 



■ 9k,A9' 



A' rni5 



,0n^), (2.23) 



i=l 



where On {Ok) stand for the finite-size rapidities related to the integers n (half-integers k) 
by the equations (I2.21I) - (I2.22I) . Fk,n is the spin-field form factor in infinite-space [88]; the 
explicit form of this function is not needed for the present purposes. The overall factor 
S{L) is given by 

(mL)2 ff°° dOid02 



S{L) = exp 



sinh^i sinh^2 



log 



2 yy_oo (27r)2 sinh(mL cosh 6'i) sinh(mL cosh 6*2 
The momentum-dependent leg factors g and g are 

g(0) = e""^^^ /VmL cosh^ , g{0) = e-""^^^ /VmL cosh^ 

where 



coth 



Oi — O2 



-mL cosh 9' 



+ e- 



-mL cosh I 



(2.24) 
(2.25) 
(2.26) 



27r cosh(^ - 0' 

It is straightforward to show, that the exact result (12.231) coincides with our formula 
(I2.18P to all orders in 1/L; the error terms are of order e""*^. First of all, the normalization 
factor S{L) does not depend on the form factor in question and for large L its behaviour 
is given by 

S{L) ^ 1 + 0(e-2™^) 
The function n{0) decays exponentially with L, therefore one has 

1 -~mL\ 1 



+ 0(e 



m 



+ 0(e-' 



rnL\ 



where pi{0) = EL = ■mcosh{0)L is the one-particle density derived from eqs. ( 12. 21^ - 
( I2.22p . Being a free theory, the products of one-particle densities equals the multi-particle 
densities entering (12.181) and one has indeed 



Ns(^i,^2,--- ,fcx|o-(0,0)|ni,n2, 

FK,N{Oki, ■ ■ ■Okj^.lOni, ■ 



R 



5 ^"jvy 



y/pK{Ok,, - ■■Okj,)p 



N[Vni 



5 ^"jvy 



-0{e 



'TnL\ 



It was shown in 12.2.31 the leading correction term to (I2.18P is at most of order e~^^ 
where p < m with m being the mass of the lightest particle. In chapter [4] we show that 
the leading term (the /z-term) is connected with the inner structure of the particles under 
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the bootstrap procedure. In fact, in most theories ^ < m, whereas in the present case 
the leading exponential correction is only of order e~™^. The reason for this is simple: 
the model in question is a free theory, there are no particle fusions. Therefore in this case 
there is no /i-term present. 

2.3 Diagonal matrix elements 

In this section we consider diagonal matrix elements of the form 

n...i„({/l---4}|^|{/l.../n})n...^„,L 

First we derive the first two (n = 1 and n = 2) cases using form factor perturbation 
theory [9Tj. Built on these rigorous results we conjecture the general formula in l2.3.2[ 

Two additional subsections are devoted to the discussion of our results. In 12.3.31 we 
examine the connection between two evaluation methods of the infinite volume diagonal 
form factors. In 12.3.41 we prove that a conjecture made by Saleur in [31] exactly coincides 
with ours. 

2.3.1 Form factor perturbation theory and disconnected contri- 
butions 

In the framework of conformal perturbation theory, we consider a model with the 
action 

A{n, \) = AcFT - fJ' j dtdx^{t,x) - \ J dtdx^{t,x) (2.27) 

such that in the absence of the coupling A, the model defined by the action ^(/i, A = 0) is 
integrable. The two perturbing fields are taken as scaling fields of the ultraviolet limiting 
conformal field theory, with left /right conformal weights hq> = h^ <1 and /i.^ = < 1, 
i.e. they are relevant and have zero conformal spin, resulting in a Lorentz-invariant field 
theory. 

The integrable limit ^(/i, A = 0) is supposed to define a massive spectrum, with 
the scale set by the dimensionful coupling /x. The exact spectrum in this case consists of 
massive particles, forming a factorized scattering theory with known S matrix amplitudes, 
and characterized by a mass scale M (which we take as the mass of the fundamental 
particle generating the bootstrap), which is related to the coupling ^ via the mass gap 
relation 

^ = kM^-^'** 

where k is a (non-perturbative) dimensionless constant. 

Switching on a second independent coupling A in general spoils integrability, deforms 
the mass spectrum and the S matrix, and in particular allows decay of the particles which 
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are stable at the integrable point. One way to approach the dynamics of the model is the 
form factor perturbation theory proposed in [9l]. Let us denote the form factors of the 
operator \i/ in the A = theory by 

= (o|^(o, 0)1^^1... e„)^=l 

Using perturbation theory to first order in A, the following quantities can be calculated 

m-- 

1. The vacuum energy density is shifted by an amount 

6S,ac = A (0| ^ |0)^^o . (2.28) 

2. The mass (squared) matrix M^^ gets a correction 

5Mf, = 2AF2^(z7r,0),5 5„ (2.29) 

(where the bar denotes the antiparticle) supposing that the original mass matrix 
was diagonal and of the form M^^ = m\5ab ■ 

3. The scattering amplitude for the four particle process a + h — > c + c? is modified by 

iS-((-.A) = -a^i(Hil±i£LMkM , 9 = (2.30) 

niamb smh (J 

It is important to stress that the form factor amplitude in the above expression must 
be defined as the so-called "symmetric" evaluation 

limFf {tir + e, ^ + m + e, 0, e),Sab 

(see eqn. (12.351) below). It is also necessary to keep in mind that eqn. (12.301) gives 
the variation of the scattering phase when the center-of-mass energy (or, the Man- 
delstam variable s) is kept fixed [91]. Therefore, in terms of rapidity variables, this 
variation corresponds to the following: 

SS^t {9, A) = 9S:t{e X = 0) ^^ ^ ^ dS:t {6, A) 



09 dX 

where 



A=0 



- _ madnia + madnia + imbdnia + niadmi,) cosh 6 

ou = 

niamb sinh 6 

is the shift of the rapidity variable induced by the mass corrections given by eqn. 
(I229j). 
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It is also possible to calculate the (partial) decay width of particles [92], but we do not 
need it here. 

We can use the above results to calculate diagonal matrix elements involving one 
particle. For simplicity we present the derivation for a theory with a single particle species. 
Let us start with the one-particle case. The variation of the energy of a stationary one- 
particle state with respect to the vacuum (i.e. the finite volume particle mass) can be 
expressed as the difference between the first order perturbative results for the one-particle 
and vacuum states in volume L: 

Am{L) = XL (({0}|^|{0})l - (0|^|0)l) (2.31) 

On the other hand, using Liischer's results [5] it only differs from the infinite volume mass 
in terms exponentially falling with L. Using eqn. (12.291) 

Am(L) = —F'^Un, 0) + O (e^^^) 
m ^ ' 

Similarly, the vacuum expectation value receives only corrections falling off exponentially 
with L. Therefore we obtain 

mmmL = ^ (F*(*vr, o) + mL{o\m) + ■■■ 

with the ellipsis denoting residual finite size corrections. Note that the factor mL is just the 
one-particle Bethe-Yang Jacobian pi{0) = mL cosh. 9 evaluated for a stationary particle 
^ = 0. 

We can extend the above result to moving particles in the following way. Up to residual 
finite size corrections, the one-particle energy is given by 

E{L) = ^m^ + p2 

with 

27r/ 



where / is the Lorentz spin (which is identical to the particle momentum quantum num- 
ber). Therefore 

EAE = niAm 

whereas perturbation theory gives: 

AE = \L{{{im{i}),-{o\mL) 

and so we obtain 

{{im{i})L = ^ (F*(^vr, 0) + p,{emm) + ■■■ (2.32) 
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where 

27r/ 



sinh^ = ^ pAO) = VnWT^n^ 

niL 

One can use a similar argument to evaluate diagonal two-particle matrix elements in 
finite volume. The energy levels can be calculated from the relevant Bethe-Yang equations 

m^jL sinh Oi + - 62) = 27r/i 
rrii^L sinh 62 + 6(62-61) = 27r/2 

and (up to residual finite size corrections) 

E2{L) = E2pt{L) — Eq{L) = rrii^ cosh6'i + rrii^ cosh ^2 

where ii and 12 label the particle species. After a somewhat tedious, but elementary 
calculation the variation of this energy difference with respect to A can be determined, 
using (l2:29|) and (l23nD : 

L 



AE2iL) = X 

Pi 



— 7 V (^f f^2 + «vr, 61 + zvr, 6u 62) + m^^L cosh 6iF^{in, 0) 

+mi2L cosh 612^2* (zTT, 0)^,^, 



«2«2 



where all quantities (such as Bethe-Yang rapidities 6i, masses and the two-particle 
state density P2) are in terms of the A = theory. This result expresses the fact that 
there are two sources for the variation of two-particle energy levels: one is the mass shift 
of the individual particles, and the second is due to the variation in the interaction. On 
the other hand, in analogy with ( 12.3ip we have 

AE2iL) = \L (,,,,({/i,/2}|vI/|{/i,/2}),,,,,,.-(0|vI/|0)z.) 

and so we obtain the following relation: 

,,,,({/l,/2}|^|{/l,/2}k.,,L = ^^^(f! (62+t7C,6l+t'K, 61,62) 

Pn^,[0^,62y ^ 
+mi-^L cosh 61F2 {iix, Q)i2i2 

+mi2L cosh d2F^ (in, 0)i^,^ + (0|^|0)) + ... (2.33) 



where the ellipsis again indicate residual finite size effects. The above argument is a 
generalization of the derivation of the mini-Hamiltonian coefficient C in Appendix C of 

m. 

We wish to remark that the result (I2.33P extends to two-particle states in non- 
integrable models below the inelastic threshold, as the Bethe-Yang equations remain valid 
in this case [5]. 
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2.3.2 Generalization to higher number of particles 

Let us now introduce some more convenient notations. Given a state 

|{Jl . . . In})h...i„ 

we denote 

p({/Ci, . . . , K})l = P^,^...^,^ek„. ..Ar) (2.34) 

where 6'/, / = 1, . . . , n are the solutions of the n-particle Bethe-Yang equations ( I2.1ip at 
volume L with quantum numbers Ii, . . . ,In and p({fci, • • • , K}, L) is the r-particle Bethe- 
Yang Jacobi determinant (I2.15p involving only the r-element subset 1 < ki < ■ ■ ■ < < n 
of the n particles, evaluated with rapidities Ok^, ■ ■ ■ , 9k^. Let us further denote 

^{{kl, . . . , kr})L = -F2r(^fci, • • • , ^fc, ) ifc^ ...ifc, 

where 

F^^iOu . . . , = hraFKOn + in + e, . . . ,9, + nr + e,9,, . . . , ^„)i,...i„i„...i, (2.35) 

is the so-called symmetric evaluation of diagonal n-particle matrix elements, which we 
analyze more closely in the next subsection. Note that the exclusion property mentioned 
at the end of subsection 2.1 carries over to the symmetric evaluation too: (I2.35P vanishes 
whenever the rapidities of two particles of the same species coincide. 

Based on the above results, we conjecture that the general rule for a diagonal matrix 
element takes the form of a sum over all bipartite divisions of the set of the n particles 
involved (including the trivial ones when A is the empty set or the complete set {1, . . . , n}): 

,,..,„({/i...4}|vl>|{/i.../J),^,..,„,^ = -— ^ — X (2.36) 

p{{l,...,n})L 

J2 -F(A)^p({l,...,n}\A)^ + 0(e-^^) 

Ac{l,2,...n} 

2.3.3 Diagonal matrix elements in terms of connected form fac- 
tors 

Here we discuss diagonal matrix elements in terms of connected form factors. The 
results will be used in the following subsection to prove that a conjecture made by Saleur 
in [31] exactly coincides with our eqn. (12. 36^ . To simplify notations we omit the particle 
species labels; they can be restored easily if needed. 
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Relation between connected and symmetric matrix elements 

The purpose of this discussion is to give a treatment of the ambiguity inherent in 
diagonal matrix elements. Due to the existence of kinematical poles (12.5p the expression 

F2n{0i + iT!-, O2 + On + ih , 9n, 6*2, 6*1) 

which is relevant for diagonal multi-particle matrix elements, is not well-defined. Let us 
consider the regularized version 

F2n{0i +i7i + ei, 92 + in + €2, 6'„ + in + e„, 6'„, 62, 9i) 

It was first observed in [91] that the singular parts of this expression drop when taking 
the limits ei ^ simultaneously; however, the end result depends on the direction of the 
limit, i.e. on the ratio of the parameters. The terms that are relevant in the limit can 
be written in the following general form: 

F2n{0i + ivr + ei, 62 + iir + €2, ^„ + ivr + e„, ft„, ...,62, 6^) = (2.37) 

n ^ n n n 
W — ■ ^ ^ ■■• ^ «iii2...in(^l5 • • • 5 ^n)^ii^i2---^i„ + ■ ■ ■ 

i = l * «! = ! 42 = 1 in = l 

where aij^i^„A„ is a completely symmetric tensor of rank n and the ellipsis denote terms 
that vanish when taking — > simultaneously. 

In our previous considerations we used the symmetric limit, which is defined by taking 
all ej equal: 

F2*„(^i, 62, On) = limF2„(^i + in + €,02 + in + e, On + in + e, 6^, ^2, ^1) 

It is symmetric in all the variables 9i, . . . ,9n. There is another evaluation with this sym- 
metry property, namely the so-called connected form factor, which is defined as the 
independent part of eqn. (I2.37p . i.e. the part which does not diverge whenever any of the 
ei is taken to zero: 

F^^{9i,e2,...,9n)=n\ai2...n (2.38) 
where the appearance of the factor n\ is simply due to the permutations of the e^. 

The relation for n < 3 

We now spell out the relation between the symmetric and connected evaluations for 
n = 1, 2 and 3. 

The n = 1 case is simple, since the two-particle form factor F2{6i, 62) has no singular- 
ities at 9i = 02 + in and therefore 



Fm=F^i9) = F2izn,0) 



(2.39) 
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It is independent of the rapidities and will be denoted F2 in the sequel. 
For n = 2 we need to consider 

F4{9i +i7r + ei, 6'2 + «vr + £2, O2, Oi) ~ (2.40) 

which gives 

^4 (6*1, 6*2) = aii + 2ai2 + a22 
^^=(^1,^2) = 2ai2 

The terms an and 022 can be expressed using the two-particle form factor. Taking an 
infinitesimal, but fixed €2 

ResF4(0i +iTc + d, 6^2 + ivr + 62, ^2, ^1) = a22e2 

ei=0 

whereas according to (12.61) 

ResF4(^i+«7r+ei, ^^2+^vr+e2, ^2, ^1) = ^ (1 - - ^2)^(^1 -62-111- £2)) F2(^2+^vr+e2, 62) 

£1=0 

To first order in €2 

Si9i -02-171- €2) = S{92 -9, + €2) = S{92 - 9{){l + i^{92 - 9{)e2 + . . . ) 

where 

^[9) = -i^^\ogS{9) 
is the derivative of the two-particle phase shift defined before. Therefore we obtain 

a22 = ^{92 - 9i)F^ 

and similarly 

an = ^{ei - 92)F^ 

and so 

Fl{9n 92) = Fl{9n 92) + 2^{9, - 92)F2{t7r, 0) (2.41) 

In the case of the trace of the energy-momentum tensor B the following expressions are 
known [SH 



F2® = 2TTm^ 



Ff = 87rmV(^i - ^2) cosh^ 



2 

Ff'" = 47rmV(^i - ^2) cosh(ei - ^2) 

and they are in agreement with f l2.4ip . 

For n = 3, a procedure similar to the above gives the following relation: 

F^{9i, 92, 9s) = F^{9u ^2, ^3) + [FM, ^2)((p(^i - ^3) + ^{02 - Os)) + permutations] 
-f3F2 [v5(6'i - 92)^{9i - 6*3) + permutations] (2.42) 

where we omitted terms that only differ by permutation of the particles. 
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Relation between the connected and symmetric evaluation in the general case 

Our goal is to compute the general expression 

F2n{0i, . . .,9n\ei, . . . ,e„) = F2n{0i+in + ei,92 + in + e2, 6'„+27r + e„, 6'„, ...,62,61) (2.43) 

Let us take n vertices labeled by the numbers 1,2, ...,n and let G be the set of the 
directed graphs Gi with the following properties: 

• Gi is tree-like. 

• For each vertex there is at most one outgoing edge. 
For an edge going from i to j we use the notation Eij. 

Theorem 1 f l2.43p can be evaluated as a sum over all graphs in G, where the contribu- 
tion of a graph Gi is given by the following two rules: 

• Let Ai = {ai,a2, . . . ,am} be the set of vertices from which there are no outgoing 
edges in Gi. The form factor associated to Gi is 



Note that since cannot contain cycles, the product of the ej/e^ factors will never be trivial 
(except for the empty graph with no edges). 

Proof The proof goes by induction in n. For n = 1 we have 



This is in accordance with the theorem, because for n = 1 there is only the trivial graph 
which contains no edges and a single node. 

Now assume that the theorem is true for n — 1 and let us take the case of n particles. 



Consider the residue of the matrix element (12.43P at e„ = while keeping all the finite 

R= ResF2n{6i..6n\ei..en) 




(2.44) 



• For each edge Ej^ the form factor above has to be multiplied by 



^{6, - 6k) 




F2^(0i) = F2(m,O) 
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12 12 12 

(a) (b) (c) 

Figure 2.1: The graphs relevant for n = 2 

According to the theorem the graphs contributing to this residue are exactly those for 
which the vertex n has an outgoing edge and no incoming edges. Let Rj be sum of the 
diagrams where the outgoing edge is Enj for some j = l,...,n — 1, and so 

n-1 

R = ^Rj 

j=i 

The form factors appearing in Rj do not depend on On- Therefore we get exactly the dia- 
grams that are needed to evaluate F2(n-i){Gi. .6n~i\^i - -^n-i) , apart from the proportionality 
factor associated to the link E^j and so 

Rj = — — ^n)-^2(n-l)(^l--^n-l|ei--en-l) 

and summing over j gives 

R = {eM^l - dn) + e2<^(^2 -On) + --- + e„_iV7(^„_l - ^n))F2(n-l)(^l-^n-l|ei-en-l) (2.45) 

In order to prove the theorem, we only need to show that the residue indeed takes this 
form. On the other hand, the kinematical residue axiom (12. Sp gives 

= ^ ^1 - n - ^^■)'^(^« - - '^^ - j ^2(n-i)(^i..^n-i|ei..e„_i) 

which is exactly the same as eqn. f l2.45p when expanded to first order in e^. 

We thus checked that the theorem gives the correct result for the terms that include 
a l/e„ singularity. Using symmetry in the rapidity variables this is true for all the terms 
that include at least one l/e^ for an arbitrary i. There is only one diagram that cannot 
be generated by the inductive procedure, namely the empty graph. However, there are no 
singularities (l/ci factors) associated to it, and it gives F|„(6'i, . . . , On) by definition. Qed. 

We now illustrate how the theorem works. For n = 2, there are only three graphs, 
depicted in figure 12. 1[ Applying the rules yields 

F,{e,, ^^2|ei, £2) = Fl{e,, 02) + ^{9, - O^) + |^ F^ 

which gives back (12.411) upon putting ei = €2. For n = 3 there are 4 different kinds 
of graphs, the representatives of which are shown in figure 12.21 all other graphs can be 
obtained by permuting the node labels 1, 2, 3. The contributions of these graphs are 
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12 3 12 3 

(a) (b) 

12 3 12 3 

(c) (d) 

Figure 2.2: The graphs relevant for n = 3 



(a) : ^6^(^1,02,^3) 

(6) : %i9i-92)F^ie2,es) 

ei 

^2 Cl 

ei ^3 

Adding up all the contributions and putting ei = €2 = 63 we recover eqn. (12.421) . 



2.3.4 Consistency with Saleur's proposal 

Saleur proposed an expression for diagonal matrix elements in terms of connected 
form factors in [3l], which is partially based on earlier work by Balog [93] and also on 
the determinant formula for normalization of states in the framework of algebraic Bethe 
Ansatz, derived by Gaudin, and also by Korepin (see [94] and references therein). To 
describe it, we must extend the normalization of finite volume states defined previously 
to the case when the particle rapidities form a proper subset of some multi-particle Bethe- 
Yang solution. 

It was shown in subsection 12.2.21 that the normalization of a finite volume state is 
given by 

. . . ,/n})L = , ^ ^ l^'l, . . .,6n) 

\l pn{Oi, . . . ,6'„) 

We again omit the particle species labels, and also denote the n-particle determinant by 
Pn- Let us take a subset of particle indices A G {1, . . . ,n} and define the corresponding 
sub-determinant by 

where J^'^ is the sub-matrix of the matrix ^7^"^ defined in eqn. (12. 15^ which is given by 
choosing the elements whose indices belong to A. The full matrix can be written explicitly 
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as 



J{n) ^ 



I EiL + (pi2-\ h (pin -^12 ■ ■ ■ -^In \ 

-^12 E2L + + V523 H \-^2n ■■■ -^2n 



\ -<^ln -<^2n ■ ■ ■ EnL + H h V^n-l.n / 

where the following abbreviations were used: Ei = rrii cosh 6^4, ipij = (pji = ip{6i — Oj). Note 
that pn depends on all the rapidities, not just those which correspond to elements of A. 
It is obvious that 

Pn{9i, ...,6n)= Pn{Ol, ■ ■ • ,^n|{l, ■ ■ ■ ,n}) 

Saleur proposed the definition 

{{ek}k^A\{ek}k^A)L = Pn{Oi, ...,9n\A) (2.46) 

where 

\{Ok}keA)L 

is a "partial state" which contains only the particles with index in A, but with rapidities 
that solve the Bethe-Yang equations for the full n-particle state. Note that this is not a 
proper state in the sense that it is not an eigenstate of the Hamiltonian since the particle 
rapidities do not solve the Bethe-Yang equations relevant for a state consisting of |A| 
particles (where \A\ denotes the cardinal number - i.e. number of elements - of the set 
A). The idea behind this proposal is that the density of these partial states in rapidity 
space depends on the presence of the other particles which are not included, and indeed 
it is easy to see that it is given by pn{Oi, . . . ,9n\A). 

In terms of the above definitions, Saleur's conjecture for the diagonal matrix element 

is 

,„„,„{{h...Inm{h...In})n...„,L = ~ ^ ^. X (2.47) 

PnKfl, ■ ■ ■ ,t>n) 

E ^A\m}keA)p{L . . .,6^) + 0(6-^^^) 
Ac{l,2,...n} 

which is just the standard representation of the full matrix element as the sum of all the 
connected contributions provided we accept eqn. (12.46p . The full amplitude is obtained 
by summing over all possible bipartite divisions of the particles, where the division is into 
particles that are connected to the local operator, giving the connected form factor F'^ 
and into those that simply go directly from the initial to the final state which contribute 
the norm of the corresponding partial multi-particle state. 

Using the results of the previous subsection, it is easy to check explicitly that our rule 
for the diagonal matrix elements as given in eqn. f l2.36[ ) is equivalent to eqn. ( 12. 47^ . We 
now give a complete proof for the general case. 
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Theorem 2 



Yl F^\A\m}keA)p{eu...,en\A) = J2 ^A\m}keA)pm} 



keN\A) 



(2.48) 



AcN 



AcN 



where we denoted N = {1,2, ... ,n}. 

Proof The two sides of eqn. ( 12.48^ differ in two ways: 

• The form factors on the right hand side are evaluated according to the „symmetric" 
prescription, and in addition to the connected part also they contain extra terms, 
which are proportional to connected form factors with fewer particles. 



The densities p on the left hand side are not determinants of the form (12.151) written 
down in terms of the particles contained in \ A: they contain additional terms 
due to the presence of the particles in A as well. 



Here we show that eqn. (12.48^ is merely a reorganization of these terms. 

For simplicity consider first the term on the left hand side which corresponds to A 
{m + l,m + 2, . . . ,n}, i.e. 

-^2m(^m+l) • • • > ^n)p(^l, ■ ■ ■ , On\A) 



We expand p in terms of the physical multi-particle densities p. In order to accomplish 
this, it is useful to rewrite the sub-matrix jTJ^yy^ as 



i=m+l 



J^''''\n\a = J'^{9,,...,9^) + 



\ 



E f2i 

i=m+l 



i=m+l 



where JT""* is the m-particle Jacobian matrix which does not contain any terms depending 
on the particles in A. The determinant of J7^\^ can be written as a sum over the subsets 
of N \ A. For a general subset -B C \ A let us use the notation B = 62, • • • , b\B\}- 
We can then write 



p{9i,...,9n\A) = detJ^''^ 



N\A 



B 



piNMAuB)) 



B 



i=l \ci=m+l 



(2.49) 



where p{N\{AUB)) is the p-density ( I2.15p written down with the particles in N\{AUB). 
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Applying a suitable permutation of variables we can generalize eqn. (12.491) to an arbi- 
trary subset A (Z N: 

\B 

p{9,, ...,9^\A)= det U\A = 5^ p(iV \ (A U S)) ^( J] <^b.,cj (2.50) 

B C i=l 

where the second summation goes over all the sets C = {ci, C2, . . . , c\b\} with \C\ = \B\ 
and Ci G A. The left hand side of eqn. (12.481) can thus be written as 



E^2|A|(WfceA)p(^i,---,^n|A) = Yl p(iV\(^Ui?))5^F(^,s42.51) 

A,B C 

An5 = 

where F(^a,b,c) = F2\A\{{9k}keA) '■Ph. 



A^BCN 



B 



i=l 



We now show that there is a one-to-one correspondence between all the terms in ( 12.511) 
and those on the right hand side of (I2.48P if the symmetric evaluations are expanded 
according to Theorem 1. To each triplet {A, B, C) let us assign the graph G(^a,b,c) defined 
as follows: 

• The vertices of the graph are the elements of the set AU B. 

• There are exactly \B\ edges in the graph, which start at hi and end at Cj with 
i = l,...,\B\. 

The contribution of (j(a,b,c) to -F2*(|A|+|_B|)({^fc}'i-'6AuB) is nothing else than F{^a,b,c) which 
can be proved by applying the rules of Theorem 1. Note that all the possible diagrams 
with at most n vertices are contained in the above list of the G(^a,b,c)^ because a general 
graph G satisfying the conditions in Theorem 1 can be characterized by writing down the 
set of vertices with and without outgoing edges (in this case B and A) and the endpoints 
of the edges (in this case C). 

It is easy to see that the factors p{N \{AU B)) multiplying the F(^a,b,c) in ( 12.511) are 
also the correct ones: they are just the density factors multiplying -F2*(|yi|_^|^|)({6'fc}fceyiuB) 
on the right hand side of (12.48p . Qed. 



2.4 Zero-momentum particles 
2.4.1 Theories with one particle species 

In theories with only a single particle species, there can only be a single particle of zero 
momentum in a multi-particle state due to the exclusion principle. For the momentum to 
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be exactly zero in finite volume it is necessary that all the other particles should come 
with quantum numbers in pairs of opposite sign, which means that the state must have 
2n + 1 particles in a configuration 

\{h, ■ ■ ■ ,ln,0, —In, • • • , — 

Therefore we consider matrix elements of the form 

({/(, ...,4,0, -4, -/(}|<i>|{/i, ...,/,, 0, -lu . . . , -h})L 

(with = or / = corresponding to a state containing a single stationary particle). We 
also suppose that the two sets {/i, . . . , Ik] and . . . , J/} are not identical, otherwise we 
have the case of diagonal matrix elements treated in section I2.3[ 
We need to examine form factors of the form 

F2k+2i+2{iT!- + 9[, . . . ,iiT + 9'i„i7!- - e'l^, . . . ,iTT - 9[,i7!- + e,0,9i, . . . , 9i, -9i, -9i) 

where the particular ordering of the rapidities was chosen to ensure that no additional S 
matrix factors appear in the disconnected terms of the crossing relation (12.21) . Using the 
singularity axiom (12.51) . plus unitarity and crossing symmetry of the S-matrix it is easy 
to see that the residue of the above function at 6' = vanishes, and so it has a finite limit 
as 6* ^ 0. However, this limit depends on direction just as in the case of the diagonal 
matrix elements considered in section 4. Therefore we must specify the way it is taken, 
and just as previously we use a prescription that is maximally symmetric in all variables: 
we choose to shift all rapidities entering the left hand state with the same amount to 
define 

Tk,i{9[, . . . , ^fc|6'i, ...,9i) = 

lim F2k+2i+2{iTr + 9[ + e, . . . ,iTT + 9'j^ + e^in - 9[ + e, . . . ,in - 9[ + e, 

^7l + e,0,9u...,9l,-9l,...,-9^) (2.52) 

Using the above definition, by analogy to ( 12.36p we conjecture that 

{{I[, ...,4,0, -4, . . . , -I[m{h, ...,/,, 0, -Ii, -h})L = (2.53) 

1 ^ 

\/ P2k+ii9[, ...,9[,0, -9[, . . . , -9[)p2i+i{9i, . . . , 6^;, 0, -9i, -9i) 
[j^k,i0i. • • • > ^'k\^u ...A)+mL F2fc+2/(^7r + 9[, . . . ,11: + 9'^,, 
m - . . . , m - 9[A, • • • , k -k • • • , -^~i)) + 0(e-'^^) 

where 9 denote the solutions of the appropriate Bethe-Yang equations at volume L, p„ is a 
shorthand notation for the n-particle Bethe-Yang density ( I2.15p and equality is understood 
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up to phase factors. We recall from section [2.2.1l that relative phases of multi-particle states 
are in general fixed differently in the form factor bootstrap and in numerical methods 
(TCSA). Also note that reordering particles gives phase factors on the right hand side 
according to the exchange axiom f l2.3p . This issue is obviously absent in the case of 
diagonal matrix elements treated in the previous section, since any such phase factor 
cancels out between the state and its conjugate. Such phases do not affect correlation 
functions, or as a consequence, any physically relevant quantities since they can all be 
expressed in terms of correlators. 

There is some argument that can be given in support of eqn. (12.531) . Note that the 
zero-momentum particle occurs in both the left and right states, which actually makes 
it unclear how to define a density similar to p in (12. 46^ . Such a density would take into 
account the interaction with the other particles. However, the nonzero rapidities entering 
of the two states are different and therefore there is no straightforward way to apply 
Saleur's recipe (12.471) here. Using the maximally symmetric definition (I2.52P the shift e 
can be equally put on the right hand side rapidities as well, and therefore we expect that 
the density factor multiplying the term F2k+2i in ( I2.53p would be the one-particle state 
density in which none of the other rapidities appear, which is exactly mL for a stationary 
particle. This is a natural guess from eqn. (I2.36P which states that when diagonal matrix 
elements are expressed using the symmetric evaluation, only densities of the type p appear. 

Another argument can be formulated using the observation that eqn. ( I2.53p is only 
valid if Tk,i is defined as in (I2.52P ; all other possible ways to take the limit can be related 
in a simple way to this definition and so the rule (12.530 can be rewritten appropriately. 
Let us consider two other natural choices 

^tii^i^ ■ ■ • • • • ,^«) = 

limF2fc+2/+2(^7r + e[,...,i7i + 6*^,^ -9'^,...,in- 9[,i7i,e,9i, ...,9i, -9i, -9i) 
^k,A9[, . . .,9'f,\9i, ...,9i) = 

limF2fc+2i+2(«vr + 9[, . . . ,i7i + 9'i„iTr - 9';^, . . . Jtt - 9[,iiT + e, 0,611, ■■■,9i, -9i, -9i) 

e— >0 

in which the shift is put only on the zero-momentum particle on the right /left, respectively. 
Using the kinematical residue axiom ( 12. 5p . JF^ can be related to T via 

J^k,i{9i, ■ ■ ■ , 9'k\9i, . . .,9i) = J^k,i{9'i-, . . . ■,9'J^^9l, . . . ,9i) 
I 

+2 J2 v{9i)F2k+2i{n^ + 9[, . . . ,17, + 9^171 - 9'k, . . . ,171 - 9[,9i, . . . ,9i, -9i, -9i) 

i=l 

J^k,i{9'i^ . . . , 9'k\9i, . . . ,9i) = J^k,ii9'i, ■ ■ ■ ,9'^.\9i, . . . ,9i) 

k 

-2Y,M)F2k+2i{ii^ + 9[,..., 17, + 91x7,-91..., 17, -9[,9^,...,9u-9u..., -9,) 

i=l 
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With the help of the above relations eqn. (12.53P can also be rewritten in terms of The 
way JF and therefore also eqn. (12.531) are expressed in terms of JF^ shows a remarkable 
and natural symmetry under the exchange of the left and right state (and correspondingly 
JF+ with JF~), which provides a further support to our conjecture. 

The above two arguments cannot be considered as a proof; we do not have a proper 
derivation of relation p.53p . On the other hand, as we show in section [3771 it agrees very 
well with numerical data which would be impossible if there were some additional (p terms 
present; such terms, as shown in section [2. 2. II would contribute corrections of order 1// in 
terms of the dimensionless volume parameter / = mL. 

2.4.2 Generalization to more than one particle species 

In a general diagonal scattering theory, there can be more than one particles of zero 
momentum in a multi-particle state. This can occur if the following two conditions hold: 

• The sets of the rapidities are "parity-symmetric" for each particle species separately, 
ie. for every particle with non-zero rapidity 6 there is a particle of the same species 
with rapidity —6. The same holds for the momentum quantum numbers. It is con- 
venient to denote such a state with 2n + m particles by 

-Ji, /a, -h, ...,In, -In, 0, 0, ... , 0}).^.^,^.^ ii...j™,L 

• The scattering phases between the zero-momentum particles satisfy 

n ^.u,(0) = l k=l...m 

1=1. ..m 

This condition is necessary for the Bethe-Yang equations to hold. 

Here we consider matrix elements between two different states \ip) and \ip') which 
satisfy the previous two conditions. The generalization of (I2.53P is straightforward: one 
has to include a disconnected term for every proper subset of those zero-momentum 
particles, which are present in both states. To precisely formulate this rule, once again 
(and in this section for the last time ^ ) we need to introduce new notations. We wish to 
remark that the result presented below is unpublished material. 

Let 1^) and \ip') be a scattering state with n + o and n' + o particles, respectively, 
where o denotes the number of those zero-momentum particles which are present in both 
states. The sets of the remaining particles consist of "parity-symmetric" pairs or they may 
also contain zero-momentum particles that are not present in the other state. For the 
rapidities of these particles we use the standard notation 6*1 . . . 6'„ and 9[ . . . 9'^, and in 
the following we do not distinguish whether a particular rapidity belongs to a pair or is 
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zero. Particle types of these first n and n' particles are denoted hy ii, ... in and i'l, . . . i^/, 
whereas for the remaining o particles we use ji, . . .jo- 
For every subset Ac {ji, ■ ■ - jo} we define 

^nn' (^) = Jim Fn+n'+2\A\ (ivT + 6*1 +£,..., ivT + 6'„ + £, jlT + E, . , ilT + E , 0^^^^^, e[, . . . , 6'^,) B 

\A\ times \A\ times 

where the label of particle types is given by 

B = {ii. ..in}AA{i[ . . . v} 

with A consisting of the elements of A listed in reverse order. 
With these notations, the conjectured general formula reads 

^/Pn+o{Ol, . . . ,9n,0, . . . , O)pn'+o{(^l, ■ ■ ■ , 0, . . . , 0) 

E J'L'{A)l[{m,L) + 0{e-'^'^) (2.54) 

Ac{ji,...jo} j<^B 

where B = {ji, . . .jo} \ A. Note, that an exchange of any two rapidities in Oi . . .6^ or 
6[. . . 6'^, yields the same phase factor in each term; the expression (I2.54p is therefore 
well-defined and is to be understood up to an overall phase factor, similar to (12.530 and 

(EED. 



Chapter 3 



Finite Volume Form Factors - 
Numerical analysis 

Here we compare the analytic predictions of the previous chapter to numerical data 
obtained by the Truncated Conformal Space Approach (TCSA). 

In 13.11 we explain in detail the methods we used to numerically determine the finite 
volume matrix elements in the Lee- Yang model and in the Ising model in a magnetic 
field. In sections 13.2113.51 we present results on form factors without disconnected terms. 
The special cases of diagonal matrix elements and form factors with zero-momentum 
particles are investigated in 13.61 and 13.71 respectively. 

3.1 Form factors from truncated conformal space 

3.1.1 Truncated conformal space approach for scaling Lee- Yang 
model 

We use the truncated conformal space approach (TCSA) developed by Yurov and 
Zamolodchikov in [86] . 

Due to translational invariance of the Hamiltonian (ll.20p . the conformal Hilbert space 
Ti can be split into sectors characterized by the eigenvalues of the total spatial momentum 

the operator Lq — Lq generates Lorentz transformations and its eigenvalue is called Lorentz 
spin. For a numerical evaluation of the spectrum, the Hilbert space is truncated by im- 
posing a cut in the conformal energy. The truncated conformal space corresponding to a 
given truncation and fixed value s of the Lorentz spin reads 

^Tcs(s, ecut) = [\i^)&'H \ (Lo - lo) \^) = s\^), (Lq + ^ - ^) 1^) = : e < e^t} 
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On this subspace, the dimensionless Hamiltonian matrix can be written as 

where energy is measured in units of the particle mass m, / = mL is the dimensionless 
volume parameter, 

= m (3.2) 

is the conformal inner product matrix and 

Bff = mzrzm.=,=, (3.3) 

is the matrix element of the operator $ at the point z = z = 1 on the complex plane 
between vectors \j) from ?iTcs('5, Ccut)- The natural basis provided by the action of 
Virasoro generators is not orthonormal and therefore G^'^'^~^ must be inserted to trans- 
form the left vectors to the dual basis. The Hilbert space and the matrix elements are 
constructed using an algorithm developed by Kausch et al. and first used in [95]. 

Diagonalizing the matrix hij we obtain the energy levels as functions of the volume, 
with energy and length measured in units of m. The maximum value of the cutoff Ccut 
we used was 30, in which case the Hilbert space contains around one thousand vectors, 
slightly depending on the spin. 



3.1.2 Exact form factors of the primary field $ in the Lee- Yang 
model 

Form factors of the trace of the stress-energy tensor were computed by Al.B. 
Zamolodchikov in [19], and using the relation 



e = iA7r(l-A)$ (3.4) 

we can rewrite them in terms of $. They have the form 

" " f(0 — ) 

Fn{ei,...,e„) = {<^)Hr.Qn{Xu...,Xn)l[ ^ (3.5) 

1=1 j=i+l * J 



with the notations 



V Jo tsinh nt 

H =(^^ 
" ' 21/2^(0) 
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The mass-gap of the theory is related to the coupling constant as 

A = 0.09704845636 ■ ■ ■ x m^^/^ (3.6) 
The exact vacuum expectation value of the field $ is 

($) = 1.239394325 ■■■xi mr'^l^ 

which can be readily obtained using (|3.6l 13.41) and also the known vacuum expectation 
value of e IH 



(6) 



vrm^ 



4^3 

The functions Qn are symmetric polynomials in the variables Xj. Defining the elementary 
symmetric polynomials of n variables by the relations 

n n 

^{x\x^) = ^x''~'of\xx,...,Xn) , aj"^ = for i > n 

they can be constructed as 

Qi = 1 , Q2 = o^? , Qz = ofaf 
Qn = o-j^V^^^Pn , n > 3 

P„ = detA<(") where M\f = a^^l,^^, , j = 1, . . . , n - 3 
Note that the one-particle form factor is independent of the rapidity: 

Ff = 1.0376434349 ■ ■ ■ x im'^^^ (3.7) 

3.1.3 Truncated fermionic space approach for the Ising model 

The conformal Ising model can be represented as the theory of a massless Majorana 
fermion with the action 

Alsing = ^ y (^^^ + ^^^) 

On the conformal plane the model has two sectors, with the mode expansions 



X]rez+i ^ Neveu-Schwarz (NS) sector 

brZ^''''^^'^ Ramond (R) sector 



and similarly for the anti-holomorphic field ip. The Hilbert space is the direct sum of a 
certain projection of the NS and R sectors, with the Virasoro content 
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The spin field a connects the NS and R sectors, and its matrix elements B^j' in the 
sector with a given conformal spin s (cf. eqn. (13.31) ) can be most conveniently computed 
in the fermionic basis using the work of Yurov and Zamolodchikov [96], who called this 
method the truncated fermionic space approach. The fermionic basis can easily be chosen 
orthonormal, and thus in this case the metrics G^^^ on the spin subspaces (cf. eqn. (|3.2I) ) 
are all given by unit matrices of appropriate dimension. Apart from the choice of basis 
all the calculation proceeds very similarly to the case of the Lee- Yang model. Energy and 
volume is measured in units of the lowest particle mass m = nii and using relation (11.250 
one can write the dimensionless Hamiltonian in the form (13. ip . The highest cutoff we 
use is Ccut = 30, in which case the Hilbert space contains around three thousand vectors 
(slightly depending on the value of the spin chosen). 

We remark that the energy density operator can be represented in the fermionic lan- 
guage as 

e = ipip 

which makes the evaluation of its matrix elements in the fermionic basis extremely simple. 

3.1.4 Ising model - Form factors of the energy density operator 

The form factors of the operator e in the Eg model were first calculated in [97] and 
their determination was carried further in [92]. The exact vacuum expectation value of 
the field e is given by [98] 

(e) = e;,|/i|«/i^ , e;, = 2.00314... 

or in terms of the mass scale m = mi 

(e) = 0.45475 ■ ■ • x m (3.8) 

The form factors are not known for the general n-particle case in a closed form, i.e. no 
formula similar to that in (13.51) exists. They can be evaluated by solving the appropriate 
polynomial recursion relations derived from the form-factor axioms. We do not present 
explicit formulae here; instead we refer to the above papers. For practical calculations we 
used the results computed by Delfino, Grinza and Mussardo, which can be downloaded 
from the Web in Mathematica format [99] . 

3.1.5 Evaluating matrix elements of a local operator O in TCSA 

Identification of multi-particle states 

Diagonalizing the TCSA Hamiltonian (13. ip yields a set of eigenvalues and eigenvectors 
at each value of the volume, but it is not immediately obvious how to select the same 
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state at different values of the volume. Therefore in order to calculate form factors it is 
necessary to identify the states with the corresponding many-particle interpretation. 

Finding the vacuum state is rather simple since it is the lowest lying state in the spin-0 
sector and its energy is given by 



where the ellipsis indicate residual finite size effects decaying exponentially fast with 
volume L and B is the bulk energy density which in the models we consider is exactly 
known (11.221 11.261) . One-particle states can be found using that their energies can be 
expressed as 



again up to residual finite size effects where s is the spin of the sector considered and i is 
the species label (every sector contains a single one-particle state for each species). 

Higher multi-particle states can be identified by comparing the measured eigenvalues 
to the levels predicted by the Bethe-Yang equations. Fixing species labels ii, . . . and 
momentum quantum numbers /i, . . . , /„, eqns. (12.111) can be solved to give the rapidities 
6i,...,6n of the particles as function of the dimensionless volume parameter / = mL. 
Then the energy of the multi-particle state in question is 



which can be compared to the spectrum. 

For each state there exists a range of the volume, called the scaling region, where L 
is large enough so that the omitted residual finite size effects can be safely neglected and 
small enough so that the truncation errors are also negligible. More precisely, the scaling 
region for any quantity depending on the volume can be defined as the volume range in 
which the residual finite size corrections and the truncation errors are of the same order of 
magnitude; since both sources of error show a dependence on the state and the particular 
quantity considered (as well as on the value of the cutoff) , so does the exact position of 
the scaling region itself. 

In the scaling region, we can use a comparison between the Bethe-Yang predictions and 
the numerical energy levels to sort the states and label them by multi-particle quantum 
numbers. An example is shown in figure [3Tl where we plot the first few states in the spin-0 
sector of the scaling Lee- Yang model and their identification in terms of multi-particle 
states is given. In this case, the agreement with the predicted bulk energy density and 
the Bethe-Yang levels in the scaling region is better than one part in 10^ for every state 
shown (with the TCSA cutoff taken at Ccut = 30). 



Eo{L) = BL + ... 




n 



E^::^^ (L) = BL + J2 cosh e, + ... 



k=l 
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one-particle 



vacuum 

i . i . 1 . t . 1 . 1 . 1 . 1 . t . 1 . 1 . 1 . t • 1 . 1 . • • • • a . • • 1 L 

- ■ 10 20 30 40 

Figure 3.1: The first 13 states in the finite volume spectrum of scaling Lee- Yang model. We 
plot the energy in units of m (with the bulk subtracted): e(Z) = {E{L) — BL)/m, against 
the dimensionless volume variable / = mL. ra-particle states are labeled by . . . ,/„), 
where the Ik are the momentum quantum numbers. The state labeled |2, 1, —3) is actually 
two- fold degenerate because of the presence of | —2, —1, 3) (up to a splitting which vanishes 
as e~', cf. the discussion in subsection 3.3). The dots are the TCSA results and the 
continuous lines are the predictions of the Bethe-Yang equations (12. lip . The points not 
belonging to any of the Bethe-Yang lines drawn are two- and three-particle states which 
are only partly contained in the first 13 levels due to line crossings, whose presence is a 
consequence of the integrability of the model. 
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Evaluation of matrix elements 

Suppose that we computed two Hamiltonian eigenvectors as functions of the volume 
L (labeled by their quantum numbers in the Bethe-Yang description f l2.1ip . omitting the 
particle species labels for brevity): 

\{h,...,In})L = J2^^iIu...Jn;m) 

i 

\{I[,...,I'k})L = Y."^, {![,..., I',; L)\j) 
j 

in the sector with spin s and spin s', respectively. Let the inner products of these vectors 
with themselves be given by 

X = J2'^,{h,...,In;L)G\f^,{h,...,In;L) 

id 

M' = 5^vI;,(/(,...,/^;L)G;fvI/,(/(,...,/^;L) 

id 

It is important that the components of the left eigenvector are not complex conjugated. In 
the Ising model we work in a basis where all matrix and vector components are naturally 
real. In the Lee- Yang model, the TCSA eigenvectors are chosen so that all of their compo- 
nents are either purely real or purely imaginary depending on whether the basis vector 
\i) is an element of the h = h = or the h = h = —1/5 component in the Hilbert space. It 
is well-known that the Lee- Yang model is non-unitary, which is reflected in the presence 
of complex structure constants as indicated in p.lOp . This particular convention for the 
structure constants forces upon us the above inner product, because it is exactly the one 
under which TCSA eigenvectors corresponding to different eigenvalues are orthogonal. We 
remark that by redefining the structure constants and the conformal inner product it is 
also possible to use a manifestly real representation for the Lee- Yang TCSA (up to some 
truncation effects that lead to complex eigenvalues in the vicinity of level crossings [86]). 
Note that the above conventions mean that the phases of the eigenvectors are fixed up to 
a sign. 

Let us consider a spinless primary field O with scaling weights Ao = Aq, which can 
be described as the matrix 

Oi;'''^ = {i\0{z,z)\j)\^^^^, , |z) GHTcs(s',e™t) , |j) eHTCs(s,e™t) 

between the two truncated conformal space sectors. Then the matrix element of O can 
be computed as 

m-'^- {{![,..., r,}\OiO,0)\{h,...,In})L = 

^) /,.;!) (3.9) 
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Figure 3.2: The vacuum expectation value of $ in finite volume. The dashed line shows 
the exact infinite volume value, while the continuous line corresponds to eqn. (13.101) . 

where the volume dependent pre-factor comes from the transformation of the primary 
field O under the exponential map ( ll.lSp and we wrote the equation in a dimensionless 
form using the mass scale m. The above procedure is a generalization of the one used by 
Guida and Magnoli to evaluate vacuum expectation values in |100| : it was extended to 
one-particle form factors in the context of the tricritical Ising model by Fioravanti et al. 
in \m\ . 

3.2 Vacuum expectation values and one-particle form 
factors 

3.2.1 Scaling Lee- Yang model 

Before the one-particle form factor we discuss the vacuum expectation value. Let us 
define the dimensionless function 

(f)(1) = -im2/5(0|$|0)L 

where the finite volume expectation value is evaluated from TCSA using (13.91) . We per- 
formed measurement of as a function of both the cutoff Ccut = 21 ... 30 and the volume 
/ = 1 ... 30 and then extrapolated the cutoff dependence fitting a function 

0(/,eeut) = 0(O + ^(Oec"uf' 
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(where the exponent was chosen by verifying that it provides an optimal fit to the data) . 
The data corresponding to odd and even values of the cutoff must be extrapolated sepa- 
rately [8], therefore one gets two estimates for the result, but they only differ by a very 
small amount (of order 10~^ at / = 30 and even less for smaller volumes). The theoretical 
prediction for 0(/) is 

0(/) = 1.239394325 ■ ■ ■ + 0(6"') 

The numerical result (after extrapolation) is shown in figure 13.21 from which it is clear 
that there is a long scaling region. Estimating the infinite volume value from the fiattest 
part of the extrapolated curve (at I around 12) we obtain the following measured value 

0(/ = oo) = 1.23938... 

where the numerical errors from TCSA are estimated to affect only the last displayed 
digit, which corresponds to an agreement within one part in 10^. 

There is also a way to compute the leading exponential correction, which was derived 
by Delfino [33| : 

= ($) + - V F^izTi, 0)uKo{mir) + ... (3.10) 

TT ^ — ' 

i 

where 

/"OO 

Ko{x) = dO cosh^e-^^""'^^ 
Jo 

is the modified Bessel-function, and the summation is over the particle species i (there is 
only a single term in the scaling Lee- Yang model). This agrees very well with the numerical 
data, as demonstrated in table 13.11 and also in figure 13.21 Using Liischer's finite-volume 
perturbation theory introduced in [4], the correction term can be interpreted as the sum of 
Feynman diagrams where there is exactly one propagator that winds around the cylinder, 
and therefore eqn. (I3.10p can be represented graphically as shown in figure 13.31 On the 
other hand, by Euclidean invariance this finite size correction coincides with the first term 
in the low-T expansion for one-point functions at finite temperature (see subsections lS.l.ll 
andE2J). 

To measure the one-particle form factor we use the correspondence (12.160 between the 
finite and infinite volume form factors to define the dimensionless function 

= -^m'/^ (/2 + (2ns)Y' (0|$|{s})i 

where is the finite volume one-particle state with quantum number I = s i.e. from 

the spin-s sector. The theoretical prediction for this quantity is 



fl{l) = 1.0376434349 ■ ■ ■ + 0(e"') 



(3.11) 
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Figure 3.3: Graphical representation of eqn. (13.101) . 
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Table 3.1: Comparison of eqn. (|3.10p to TCSA data 
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Figure 3.4: One-particle form factor from sectors with spin s = 0,1,2. The continuous 
line shows the exact infinite volume prediction. 

The numerical results (after extrapolation in the cutoff) are shown in figure 13.41 The 
scaling region gives the following estimates for the infinite volume limit: 

/0(/ = oo) = 1.037654... 
/i(/ = oo) = 1.037650... 
/2(/ = oo) = 1.037659... 

which show good agreement with eqn. (13.111) (the relative deviation is again around 10~^, 
as for the vacuum expectation value). 

3.2.2 Ising model in magnetic field 

For the Ising model, we again start with checking the dimensionless vacuum expecta- 
tion value for which, using eqn. (13.81) we have the prediction 

0(/) = -{e)L = 0.45475 ■■■ + (e"') 
m 

where m = mi is the mass of the lightest particle and I = mL as before. The TCSA 
data are shown in figure 13.51 Note that there is substantial dependence on the cutoff Ccut 
and also that extrapolation in Ccut is really required to achieve good agreement with the 
infinite volume limit. Reading off the plateau value from the extrapolated data gives the 
estimate 

— (e) = 0.4544... 

m 



58 CHAPTER 3. FINITE VOLUME FORM FACTORS - NUMERICAL ANALYSIS 



46 , 

44 h I ' " ' ♦ • ♦ ♦ J ♦ ♦ I ♦ ; ; I 



0.42 
0.4 

0.38 



0.34 
0.32 



■ extrapolated 
• Ccut = 30 
Ccut = 28 

h ^ e,,t = 26 



"MM,,,,. 

' ▼ T 



10 15 20 25 30 35 40 



Figure 3.5: Measuring the vacuum expectation value of e in the Ising model 

for the infinite volume vacuum expectation value, which has 8 ■ 10~^ relative deviation 
from the exact result. Our first numerical comparison thus already tells us that we can 
expect much larger truncation errors than in the Lee- Yang case. It is also clear from figure 
I3.5l that in order to attain suitable precision in the Ising model extrapolation in the cutoff 
is very important. 

Defining the function 

m = {e)L/{e) 

we can calculate the leading exponential correction using eqn. (13.101) and the exact two- 
particle form factors from [99]. It only makes sense to include particles i = 1,2,3 since 
the contribution of the fourth particle is sub-leading with respect to two-particle terms 
from the lightest particle due to m4 > 2mi. The result is shown in figure [3T6t we do not 
give the data in numerical tables, but we mention that the relative deviation between the 
predicted and measured value is better than 10~^ in the range 5 < / < 10. 

From now on we normalize all form factors of the operator e by the infinite volume 
vacuum expectation value (13.81) . i.e. we consider form factors of the operator 

^ = (3.12) 

which conforms with the conventions used in [92l [99]. We define the dimensionless one- 
particle form factor functions as 

i^^(0= (^(^)' + (2^-)') ^ (oi^lMkx 

In the plots of figure 13.71 we show how these functions measured from TCSA compare to 
predictions from the exact form factors for particles i = 1, 2, 3 and spins s = 0, 1, 2, 3. 

It is evident that the scaling region sets in much later than for the Lee- Yang model; 
therefore for the Ising model we do not plot data for low values of the volume (all plots 
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Figure 3.6: The volume dependence of the vacuum expectation value of e in the Ising 
model, showing the extrapolated value and the prediction from eqn. ( 13.10p . normalized 
by the value in the infinite volume limit. 

start from / ~ 10... 15). This also means that truncation errors in the scaling region are 
also much larger than in the scaling Lee- Yang model; we generally found errors larger by 
an order of magnitude after extrapolation in the cutoff. We remark that extrapolation 
improves the precision by an order of magnitude compared to the raw data at the highest 
value of the cutoff. 

Note the rather large finite size correction in the case of ^3. This can be explained 
rather simply as the presence of a /i-term. Based on Liischer's finite- volume perturbation 
theory [4] it is expected, that similar to the /i-term for the masses of stable particles (II. lip , 
the leading contribution (which is associated to the diagram depicted in figure 13. 8p has 
the volume dependence 



Therefore we can expect a contribution suppressed only by e'^ . A numerical fit of the 
/-dependence in the s = case is perfectly consistent with this expectation. As a result, 
no scaling region can be found, because truncation errors are too large in the volume 
range where the exponential correction is suitably small. We do not elaborate on this 
issue further here: chapter [4] is devoted to the study of exponential corrections. 
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Figure 3.7: One-particle form factors measured from TCSA (dots) compared to the 
infinite-volume prediction from exact form factors. All numerical data have been extrap- 
olated to Ccut = oo and s denotes the Lorentz spin of the state considered. The relative 
deviation in the scaling region is around 10""^ for Ai and A2; there is no scaling region for 
/I3 (see the discussion in the main text). 
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Figure 3.8: Leading finite size correction (a so-called /x-term) to the one-particle form 
factor of ^3, which results from the process of splitting up into two copies of Ai which 
then wind around the cylinder once before recombining into A^ again. 

3.3 Two-particle form factors 

3.3.1 Scaling Lee- Yang model 

Following the ideas in the previous subsection, we can again define a dimensionless 
function for each two-particle state as follows: 



Relation (I2.16P gives the following prediction in terms of the exact form-factors: 



/2(0/i/2 = 



(0|$|{/i,/2})l 



/ = mL 



—im 



2/5 



/2(0/l/2 




Ffm)A{l)) + 0{e-^) 



(3.13) 



where 6'i(/), 6*2 (/) solve the Bethe-Yang equations 



Zsinh^i + 5(^1 
I sinh h + 5{e2 



02) 



2ttL 



27r/2 



and the density of states is given by 



pni9i,92) = cosh 6*1 cosh 6^2 + / cosh 6'iV5(6'2 - 6i) + / cosh 6'2V5(6'i - 62) 



the phase shift 6 is defined according to eqn. fll.24p and 
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There is a further issue to take into account: the relative phases of the multi-particle states 
are a matter of convention and the choice made in subsection 13.1.51 for the TCSA eigen- 
vectors may differ from the convention adapted in the form factor bootstrap. Therefore in 
the numerical work we compare the absolute values of the functions f2{l) computed from 
TCSA with those predicted from the exact form factors. Note that this issue is present 
for any non-diagonal matrix element, and was in fact tacitly dealt with in the case of 
one-particle matrix elements treated in subsection 13.2. 1[ 

The prediction (13.131) for the finite volume two-particle form factors is compared with 
spin-0 states graphically in figure [33] and numerically in table 13.21 while the spin-1 and 
spin-2 case is presented in figure 13.101 and in table 13.31 These contain no more than 
a representative sample of our data: we evaluated similar matrix elements for a large 
number of two-particle states for values of the volume parameter / running from 1 to 30. 
The behaviour of the relative deviation is consistent with the presence of a correction of 
e~' type up to / ~ 9 ... 10 (i.e. the logarithm of the deviation is very close to being a linear 
function of /), and after / ~ 16 ... 18 it starts to increase due to truncation errors. This 
is demonstrated in figure [3TT] using the data presented in table 13.31 for spin-1 and spin-2 
stateJl], but it is equally valid for all the other states we examined. In the intermediate 
region / ~ 10 ... 16 the two sources of numerical deviation are of the same order, and so 
that range can be considered as the optimal scaling region: according to the data in the 
tables agreement there is typically around 10"'^ (relative deviation). It is also apparent 
that scaling behaviour starts at quite low values of the volume (around / ~ 4 the relative 
deviation is already down to around 1%). 

It can be verified by explicit evaluation that in the scaling region the Bethe-Yang den- 
sity of states (p) given in ( 12.15p differs by corrections of relative magnitude 10^^ — 10^^ 
(analytically: of order l/l) from the free density of states (p°) in (12.140 . and therefore 
without using the proper interacting density of states it is impossible to obtain the preci- 
sion agreement we demonstrated. In fact the observed 10~^ relative deviation corresponds 
to corrections of order at I = 10, but it is of the order of estimated truncation errorj^. 

These results are very strong evidence for the main statement in (I3.13P (and thus also 
(I2.16P ). namely, that all 1/L corrections are accounted by the proper interacting state 
density factor and that all further finite size corrections are just residual finite size effects 
decaying exponentially in L. In section 13.41 we show that data from higher multi-particle 
form factors fully support the above conclusions drawn from the two-particle form factors. 



^Note that the dependence of the logarithm of the deviation on the volume is not exactly Hnear because 
the residual finite size correction can also contain a factor of some power of I, and so it is expected that 
a log I contribution is also present in the data plotted in figure 13.111 

^Truncation errors can be estimated by examining the dependence of the extracted data on the cutoff 
Ccut, as well as by comparing TCSA energy levels to the Bethe-Yang predictions. 
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Figure 3.9: Two-particle form factors in the spin-0 sector. Dots correspond to TCSA data, 
while the lines show the corresponding form factor prediction. 
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Table 3.2: Two-particle form factors |/2(0I the spin-0 sector 



64 CHAPTER 3. FINITE VOLUME FORM FACTORS - NUMERICAL ANALYSIS 




Figure 3.10: Two-particle form factors in the spin-1 and spin-2 sectors. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 





h = 3/2 , /2 = -1/2 


= 5/2 , /a = -3/2 


h = 5/2 , h = -1/2 


Ji = 7/2 , I2 = -3/2 


I 


TCSA 


FF 


TCSA 


FF 


TCSA 


FF 


TCSA 


FF 


2 


0.077674 


0.089849 


0.048170 


0.054711 


0.064623 


0.074763 


0.042031 


0.047672 


4 


0.072104 


0.073571 


0.049790 


0.050566 


0.062533 


0.063932 


0.044034 


0.044716 


6 


0.056316 


0.056444 


0.045031 


0.045100 


0.051828 


0.052009 


0.040659 


0.040724 


8 


0.042051 


0.042054 


0.039191 


0.039193 


0.041370 


0.041394 


0.036284 


0.036287 


10 


0.031146 


0.031144 


0.033469 


0.033467 


0.032757 


0.032759 


0.031850 


0.031849 


12 


0.023247 


0.023245 


0.028281 


0.028279 


0.026005 


0.026004 


0.027687 


0.027684 


14 


0.017619 


0.017616 


0.023780 


0.023777 


0.020802 


0.020799 


0.023941 


0.023936 


16 


0.013604 


0.013599 


0.019982 


0.019977 


0.016808 


0.016802 


0.020659 


0.020652 


18 


0.010717 


0.010702 


0.016831 


0.016822 


0.013735 


0.013724 


0.017835 


0.017824 


20 


0.008658 


0.008580 


0.014249 


0.014227 


0.011357 


0.011337 


0.015432 


0.015413 



Table 3.3: Two-particle form factors |/2(0I the spin-1 and spin-2 sectors 
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Figure 3.11: Estimating the error term in (13.13^ using the data in table 13.31 The various 
plot symbols correspond to the same states as specified in figure 13.101 

3.3.2 Ising model in magnetic field 

In this case, there is some further subtlety to be solved before proceeding to the 
numerical comparison. Namely, there are spin-0 states which are parity refiections of each 
other, but are degenerate according to the Bethe-Yang equations. An example is the state 
|{1, — l})i2 in figure [3?T2] (b). which is degenerate with |{ — 1, l})i2 to all orders in 1/L. In 
general the degeneracy of these states is lifted by residual finite size effects (more precisely 
by quantum mechanical tunneling - a detailed discussion of this mechanism was given 
in the framework of the fc-folded sine-Gordon model in |102j ). Since the finite volume 
spectrum is parity symmetric, the TCSA eigenvectors correspond to the states 



|{l,-l})±,i = ^ 



(|{1,-1})i2,l±|{-1,1})i2,l) 



Because the Hilbert space inner product is positive definite, the TCSA eigenvectors 
|{1,— l})f2 can be chosen orthonormal and the problem can be resolved by calculating 
the form factor matrix element using the two-particle state vectors 



V2 



(l{l,-l})^2,L±l{l,-l})r2,L) 



Because the Ising spectrum is much more complicated than that of the scaling Lee- Yang 
model (and truncation errors are larger as well), we only identified two-particle states 
containing two copies of Ai , or an Ai and an A2. The numerical results are plotted 
in figures [3TT2] (a) and (b), respectively. The finite volume form factor functions of the 
operator (|3.12p are defined as 



111 (0 



Pnm),e2mm{h,i2})u 
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where 



/ sinh 02 + 5ii(^2 - ^1) = 27r/2 

Pii(^i, ^2) = cosh6'i cosh 02 + / cosh 6iLpii{02 — Oi) +1 cosh 6'2V5ii (6*1 — 62) 



and 
with 



/12 = V/^ll(^l'^2)(0|vl/|{/i,/2})l2 



/sinh^i + (5i2(^i -^2) = 27r/i 

— / sinh 02 + (5i2(^2 - ^1) = 27r/2 
mi 

Pi2(^i, ^2) = — cosh 6^1 cosh 62 + / cosh 6'iV9i2(^^2 - ^1) + — ^ cosh 62'-Pi2{0i - 62) 
nil fT^i 

and are compared against the form factor functions 

F^{di{l)A{l))ii 

and 

F^{ei{i)M))i2 

respectively. 

Although (as we already noted) truncation errors in the Ising model are much larger 
than in the Lee- Yang case, extrapolation in the cutoff improves them by an order of 
magnitude compared to the evaluation at the highest cutoff (in our case 30). After ex- 
trapolation, deviations in the scaling region become less than 1% (with a minimum of 
around 10~^ in the AiAi^ and 10~^ in the A1A2 case), and even better for states with 
nonzero total spin. As noted in the previous subsection this means that the numerics is 
really sensitive to the dependence of the particle rapidities and state density factors on 
the interaction between the particles; generally the truncation errors in the extrapolated 
data are about two orders of magnitude smaller than the interaction corrections. 

It is a general tendency that the agreement is better in the sectors with nonzero spin, 
and the scaling region starts at smaller values of the volume. This is easy to understand for 
the energy levels, since for low-lying states nonzero spin generally means higher particle 
momenta. The higher the momenta of the particles, the more the Bethe-Yang contribu- 
tions dominate over the residual finite size effects. This is consistent with the results of 
Rummukainen and Gottlieb in |103| where it was found that resonance phase shifts can 
be more readily extracted from sectors with nonzero momentum; our data show that this 
observation carries over to general matrix elements as well. 
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Figure 3.12: Two-particle form factors in the Ising model. Dots correspond to TCSA data, 
while the lines show the corresponding form factor prediction. 
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3.4 Many-particle form factors 

3.4.1 Scaling Lee- Yang model 

We also performed numerical evaluation of three and four-particle form factors in 
the scaling Lee- Yang model; some of the results are presented in figures 13.131 and 13.141 
respectively. For the sake of brevity we refrain from presenting explicit numerical tables; 
we only mention that the agreement between the numerical TCSA data and the prediction 
from the exact form factor solution is always better than 10"^ in the scaling region. For 
better visibility we plotted the functions 

fk{l)h...i, = -tm^^'^pkiOu . . . , ^fe)(0|$|{/i, . . . , h})L , l = mL 

for which relation (I2.16P gives: 

fkil)h...i, = -tm'/'F^ie,, ...,ek) + 0(e-') (3.14) 

Due to the fact that in the Lee- Yang model there is only a single particle species, we 
introduced the simplified notation p„ for the ra-particle Jacobi determinant. 

The complication noted in subsection 3.2.2 for the Ising state |{1, — l})i2 is present in 
the Lee- Yang model as well. The Bethe-Yang equations give degenerate energy values for 
the states . . . , Ik})L and \{—Ik, ■ ■ ■ , —h})L (as noted before, the degeneracy is lifted 
by quantum mechanical tunneling). For states with nonzero spin this causes no problem, 
because these two states are in sectors of different spin (their spins differ by a sign) and 
similarly there is no difficulty when the two quantum number sets are identical, i.e. 

{/i, . . . , Jfc} = {— /i, . . . , —h} 

since then there is a single state. However, there are states in the zero spin sector (i.e. 
with Ylk = 0) for which 

We use two such pairs of states in our data here: the three-particle states |{3, —1, — 2})^, 
|{2, 1, -3})l and the four-particle states |{7/2, 1/2, -3/2, -5/2})^, |{5/2, 3/2, -1/2, -7/2})^. 
Again, the members of such pairs are related to each other by spatial refiection, which 
is a symmetry of the exact finite- volume Hamiltonian and therefore (supposing that the 
eigenvectors are orthonormal) the finite volume eigenstates correspond to 

lUi, • • • , 4})! = ^ (lUi, . . • , h])L ± . . . , -/i})l) 

and this must be taken into account when evaluating the form factor matrix elements. In 
the Lee- Yang case, however, the inner product is not positive definite (and some nonzero 




Figure 3.13: Three-particle form factors in the spin-0 sector. Dots correspond to TCSA 
data, while the lines show the corresponding form factor prediction. 

vectors may have zero "length", although this does not happen for TCSA eigenvectors, 
because they are orthogonal to each other and the inner product is non-degenerate), but 
there is a simple procedure that can be used in the general case. Suppose the two TCSA 
eigenvectors corresponding to such a pair are vi and V2. Then we can define their inner 
product matrix as 

using the TCSA inner product (13.211 . The appropriate basis vectors of this two-dimensional 
subspace, which can be identified with |{Ji, . . . , Ik})L and |{— /it, • . . , — can be found 
by solving the two-dimensional generalized eigenvalue problem 

g ■ w = \P ■ w 

for the vector {101,102} describing orientation in the subspace, with 




This procedure has the effect of rotating from the basis of parity eigenvectors to basis 
vectors which are taken into each other by spatial refiection. 
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Figure 3.14: Four-particle form factors in the spin-0 sector. Dots correspond to TCSA 
data, while the lines show the corresponding form factor prediction. 

3.4.2 Ising model in a magnetic field 

As we already noted, it is much harder to identif}!^ higher states in the Ising model 
due to the complexity of the spectrum, and so we only performed an analysis of states 
containing three Ai particles. We define 



/ill (0/, 



Pium), HI), 12, h})iu 



where 9i{l) are the solutions of the three-particle Bethe-Yang equations in (dimensionless) 
volume / and pm is the appropriate 3-particle determinant. The results of the comparison 
can be seen in figure [3T5l The numerical precision indicated for two-particle form factors 
at the end of subsection 3.2.2, as well as the remarks made there on the spin dependence 
apply here as well; we only wish to emphasize that for Ai^i^i states with nonzero total 
spin the agreement between the extrapolated TCSA data and the form factor prediction 
in the optimal part of the scaling region is within 2 x 10~^. 



^To identify AiAiAi states it is necessary to use at least Ccut = 22 or 24 and even then the agreement 
with the Bethe-Yang prediction is still only within 20%, but the identification can be made for the first few 
AiAiAi states using data up to ecut — 30. Truncation errors are substantially decreased by extrapolation 
to ecut = 00. 
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Figure 3.15: Three-particle form factors in the Ising model. Dots correspond to TCSA 
data, while the lines show the corresponding form factor prediction. 

3.5 General form factors without disconnected pieces 

3.5.1 Scaling Lee- Yang model 

In this model there is a single particle species, so we can introduce the following 
notations: 

hnii)t::j: = -^^'^'{{li, • • • , ik}\m o)i{/i, . . . , 14) l 

and also 

fkn(i)t::!i = -^^'^'\/pk{e[, . . . , h)^Pn{ei. . . . , ~en){{i[. . . . ,/ai$(o,o)i{/i, . . . ,4})^ 

for which relation (12.18^ yields 

^p^{e,,...X)Pk{e[,....e'j 

hn{l)t:± = -tm'/'F^^M + ^^,---J'i+^^Ju---Jn) + Oie-') (3.15) 

For the plots we chose to display / or / depending on which one gives a better visual 
picture. The numerical results shown here are just a fraction of the ones we actually 
obtained, but all of them show an agreement with precision 10~^ — 10~^ in the scaling 
region (the volume range corresponding to the scaling region typically varies depending 
on the matrix element considered due to variation in the residual finite size corrections 
and truncation effects). 

The simplest cases involve one and two-particle states: the one-particle-one-particle 
data in figure [3T6l actually test the two-particle form factor F2 , while the one-particle- 
two-particle plot 13.171 corresponds to F* (we obtained similar results on Ff using matrix 
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Figure 3.16: One-particle-one-particle form factors in Lee- Yang model. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 

elements /22)- Note that in contrast to the comparisons performed in subsections 3.2 and 
3.3, these cases involve the form factor solutions (13.51) at complex values of the rapidities. 
In general, all tests performed with TCSA can test form factors at rapidity arguments 
with imaginary parts or tt, which are the only parts of the complex rapidity plane where 
form factors eventually correspond to physical matrix elements. 

One-particle-three-particle and one-particle-four-particle matrix elements /13 and fu 
contribute another piece of useful information. We recall that there are pairs of parity- 
related states in the spin-0 factors which we cannot distinguish in terms of their elementary 
form factors. In subsection 3.3 we showed the example of the three-particle states 

|{3,-l,-2})i and \{2,1,-3})l 

and the four-particle states 

|{7/2, 1/2, -3/2, -5/2})l and |{5/2, 3/2, -1/2, -7/2})z. 

In fact it is only true that they cannot be distinguished if the left state is parity-invariant. 
However, using a one-particle state of nonzero spin on the left it is possible to distinguish 
and appropriately label the two states, as shown in figures 13.181 and 13.191 This can also 
be done using matrix elements with two-particle states of nonzero spin: the two-particle- 
three-particle case /23 is shown in 13.201 (similar results were obtained for 724). Examining 
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Figure 3.17: One-particle-two-particle form factors in Lee- Yang model. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 

the data in detail shows that the identifications provided using different states on the left 
are all consistent with each other. 

It is also interesting to note that the /14 (figure I3.19P and data (figure 13.200 
provide a test for the five-particle form factor solutions F5. This is important since it 
is progressively harder to identify many-particle states in the TCSA spectrum for two 
reasons. First, the spectrum itself becomes more and more dense as we look for higher 
levels; second, the truncation errors grow as well. Both of these make the identification of 
the energy levels by comparison with the predictions of the Bethe-Yang equations more 
difficult; in the Lee- Yang case we stopped at four-particle levels. However, using general 
matrix elements and the relations (13.150 we can even get data for form factors up to 8 
particles, a sample of which is shown in figures [3T2T] (/as and /44, corresponding to 6 and 
8 particle form-factors) and 13.221 {f^^ which corresponds to 7 particle form factors). 

3.5.2 Ising model in magnetic field 

In the case of the Ising model, we define the functions 

/il...j„;il...i„(0 = \/ Pii...i„(6'l, . . . , On)Pn...j^{0'l, ■ ■ ■ ,^'m)>^ h-jraiWl, ■ ■ ■ J'm}\^\{h, ■ ■ • , /n})ii...i„,L 

which are compared against form factors 

where 9i and 9j denote the rapidities obtained as solutions of the appropriate Bethe-Yang 
equations at the given value of the volume. We chose states for which the necessary form 
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Figure 3.18: One-particle-three-particle form factors in Lee- Yang model. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 
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Figure 3.19: One-particle-four-particle form factors in Lee- Yang model. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 
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Figure 3.20: Two-particle-three-particle form factors in Lee- Yang model. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 




Figure 3.21: Three-particle-tliree-particle and four-particle-four-particle form factors in 
Lee- Yang model. Dots correspond to TCSA data, while the lines show the corresponding 
form factor prediction. 
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Figure 3.22: Three-particle-four-particle form factors in Lee- Yang model. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 



factor solution was already known (and given in [99]) i.e. we did not construct new form 
factor solutions ourselves. 

One-particle-one-particle form factors are shown in figure 13.231 these provide another 
numerical test for the two-particle form factors examined previously in subsection 3.2.2. 
One-particle-two-particle form factors, besides testing again the three-particle form factor 
Ai^iAi (figure [3^241 (a)) also provide information on A1A1A2 (figure [3^241 (b)). 

Finally, one-particle-three-particle and two-particle-two-particle matrix elements can 
be compared to the AiAi^iAi form factor, which again shows that by considering general 
matrix elements we can go substantially higher in the form factor tree than using only 
elementary form factors. 

We remark that the cusps on the horizontal axis in the form factor plots correspond 
to zeros where the form factors change sign; they are artifacts introduced by taking the 
absolute value of the matrix elements. The pattern of numerical deviations between TCSA 
data and exact form factor predictions is fully consistent with the discussion in the closing 
paragraphs of subsections 3.2.2 and 3.3.2. The deviations in the scaling region are around 
1% on average, with agreement of the order of 10^^ in the optimal range. 
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Figure 3.23: One-particle-one-particle form factors in the Ising model. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 
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Figure 3.24: One-particle-two-particle form factors in the Ising model. Dots correspond 
to TCSA data, while the lines show the corresponding form factor prediction. 
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Figure 3.25: One-particle-three-particle and two-particle-two-particle form factors in 
Ising model. Dots correspond to TCSA data, while the lines show the corresponding fc 
factor prediction. 
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3.6 Diagonal matrix elements 

Here we compare the predictions of section 12.31 for the diagonal form factors to the 
numerical data obtained from TCSA. 

3.6.1 Diagonal one-particle and two-particle form factors 

Figure (I3.26h shows the comparison of eqn. (12.321) to numerical data obtained from 
Lee- Yang TCSA: the matching is spectacular, especially in the so-called scaling region (the 
volume range where residual finite size corrections are of the order of truncation errors, cf. 
|104| ) where the relative deviation is less than 10""^. Diagonal one-particle matrix elements 
for the Ising model are shown in figure 13.271 

Formula f l2.33p describing diagonal two-particle matrix elements is tested against nu- 
merical data in the Lee- Yang model in figure 13.281 and the agreement is as precise as 
it was for the one-particle case. Similar results can be found in the Ising case; they are 
shown in figure 13.291 

3.6.2 The general result 

Formula 12.361 can be tested against matrix elements with n = 3 and n = 4 in the Lee- 
Yang model, which are displayed in figures 13.301 and 13.311 respectively. The agreement is 
excellent as before, with the relative deviation in the scaling region being of the order of 
10-^ 

3.7 Zero-momentum particles 

In this section we test the predictions of sec. 12.41 for disconnected pieces associated to 
zero-momentum particles. 

3.7.1 Lee- Yang model 

First we test eqn. (I2.53P using low-lying symmetric states with a zero-momentum 
particle present. Data for matrix elements of the type 

({0}|$|{-/,0,/})^ and ({-/', 0,/'}|$|{-/,0,/})^ 

are shown in figures [3.321 and [3.331 respectively. The agreement is precise as in all previous 
cases. 

The support for eqn. ( 12. 53^ can be strengthened using 5-particle states. It is not easy 
to find them because they are high up in the spectrum, and identification using the 



3.7. ZERO-MOMENTUM PARTICLES 



81 



/ii 




5 10 15 20 25 30 



Figure 3.26: Diagonal 1-particle matrix elements in the scaling Lee- Yang model. The 
discrete points correspond to the TCSA data, while the continuous line corresponds to 
the prediction from exact form factors. 

process of matching against Bethe-Yang predictions (as described in subsection 13.1.50 
becomes ambiguous. We could identify the first 5-particle state by combining the Bethe- 
Yang matching with predictions for matrix elements with no disconnected pieces given by 
eqn. (I2.18p . as shown in figure 1X341 Some care must be taken in choosing the other state 
because many choices give matrix elements that are too small to be measured reliably 
in TCSA: since vector components and TCSA matrices are mostly of order 1 or slightly 
less, getting a result of order 10""^ or smaller involves a lot of cancellation between a large 
number of individual contributions, which inevitably leads to the result being dominated 
by truncation errors. Despite these difficulties, combining Bethe-Yang level matching with 
form factor evaluation we could identify the first five-particle level up to / = 20. 

The simplest matrix element involving a five-particle state and zero-momentum dis- 
connected pieces is the 1-5 one, but the prediction of eqn. ( 12.530 turns out to be too small 
to be usefully compared to TCSA. However, it is possible to find 3-5 matrix elements 
that are sufficiently large, and the data shown in figure [3351 confirm our conjecture with 
a relative precision of somewhat better than 10~^ in the scaling region. 

We close by noting that since the agreement is better than one part in 10^ in the 
scaling region, which is typically found in the range of volume / ~ 10 ... 20, and also this 
precision holds for quite a large number of independent matrix elements, the presence of 
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Figure 3.27: Diagonal 1-particle matrix elements in the Ising model. The discrete points 
correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 
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Figure 3.28: Diagonal 2-particle matrix elements in the scaling Lee- Yang model. The 
discrete points correspond to the TCSA data, while the continuous line corresponds to 
the prediction from exact form factors. 
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Figure 3.29: Diagonal 2-particle matrix elements in the Ising model. The discrete points 
correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 
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Figure 3.30: Diagonal 3-particle matrix elements in the scaling Lee- Yang model. The 
discrete points correspond to the TCSA data, while the continuous line corresponds to 
the prediction from exact form factors. 
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Figure 3.31: Diagonal 4-particle matrix elements in the scaling Lee- Yang model. The 
discrete points correspond to the TCSA data, while the continuous line corresponds to 
the prediction from exact form factors. 
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Figure 3.32: l-particle-3-particle matrix elements in the scaling Lee- Yang model. The 
discrete points correspond to the TCSA data, while the continuous line corresponds to 
the prediction from exact form factors. 
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Figure 3.33: 3-particle-3-particle matrix elements in the scaling Lee- Yang model. The 
discrete points correspond to the TCSA data, while the continuous line corresponds to 
the prediction from exact form factors. 
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Figure 3.34: Identifying the 5-particle state using form factors. The discrete points cor- 
respond to the TCSA data, while the continuous hue corresponds to the prediction from 
exact form factors. 
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Figure 3.35: 3-particle-5-particle matrix elements in the scaling Lee- Yang model. The 
discrete points correspond to the TCSA data, while the continuous line corresponds to 
the prediction from exact form factors. 
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Figure 3.36: Ai — Ai^i^i matrix element in Ising model with a zero-momentum particle 
additional ip terms in eqn. (12.53p can be confidently excluded. 



3.7.2 Ising model in magnetic field 

We first test our analytic results on the example of the matrix element 

i({0}|$|{-l,0,l}),_ 



Since all particles are of species Ai, formula (I2.53P is applicable here. In figure [37361 we 
plot the numerical results against the analytic prediction. 

Due to the fact that the Ising model has more than one particle species, it is possible 
to have more than one stationary particles in the same state. Our TCSA data allow us 
to locate one such state, with a stationary Ai and A2 particle, for which we have the 
prediction 

1 



1,12 



_({0}|vl/|{0,0})i2 



limF3(i7r + e, 0, 0)112 + miL ^1(0)2 



where -Fi(0)2 is the one-particle form factor corresponding to A2. This is compared to 
TCSA data in figure 13.371 and a convincing agreement is found. 

Note that in both of figures 13.361 and 13.371 there is a point which obviously deviates 
from the prediction. This is a purely technical issue, and is due to the presence of a line 
crossing close to this particular value of the volume which makes the cutoff dependence 
more complicated and so slightly upsets the extrapolation in the cutoff. 
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Figure 3.37: Ai — A1A2 matrix element in Ising model with zero-momentum particle 



Chapter 4 

Residual finite size effects 



The Bethe-Yang equations describe multi-particle energies in finite volume to all orders 
in 1/L. The same holds for the results presented in section 2: formulas (12.18p . (12.36P and 
(I2.53P determine the finite volume form factors to all orders in The aim of this 
chapter is to derive the /i-term (which is the leading exponential correction) associated 
to moving one-particle states, generic multi-particle scattering states, and finite volume 
form factors. 

The analytic results are derived using the principle of bound state quantization in 
finite volume. To confirm the calculations, all our formulas are tested against TCSA data. 

We wish to remark, that there has been a recent interest in exponential corrections also 
in the framework of AdS/CFT correspondence. The /i-term and F-term for moving one- 
particle states was derived in |105| and |106| . for multi-particle states in |107| and |108| . 
For the /i-term these results coincide with ours, but they can be considered as independent; 
they use completely different methods: generalizations of Liischer's approach, summation 
of the vacuum fluctuations and calculations based on the TBA for excited states. 

4.1 One-particle states 

The energy of a moving one-particle state is given to all orders by 



Here we determine the leading /x-term associated to (14. ip Based on the description of mass 
corrections it is expected that this contribution is associated to the fusion AaAi, — > Ac 
with the smallest /x^^. We assume that a = b, ie. the fusion in question is a symmetric one. 
This happens to be true for the lightest particle in models with the "$^-property" and 
for other low lying states in most known models. A possible extension to non-symmetric 
fusions has not yet been carried out, it is left for further research. 




Pc = 



L 



(4.1) 
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The (infinite volume) bootstrap principle for a symmetric fusion consists of the iden- 
tification 

resulting in 

me = 2ma cos«J {fil^f = ml-^ 

Smallness of fi'^^ means that rric is close to 2ma, in other words the binding energy is 
small. 

For a moment let us lay aside the framework of QFT and consider quantum mechanics 
with an attractive potential. Bound states are described by wave functions 

^(xi,X2) = e^^(^^+^2V(^i - 2:2) 

where P is the total momentum and iIj{x) is the appropriate solution of the Schrodinger 
equation in the relative coordinate. It is localized around x = and shows exponential 
decay at infinity. Except for the region Xi ~ X2, the wave function can be approximated 
with a product of plane waves with imaginary momenta pi 2 = P ±ik. The interaction 
results in the quantization of the allowed values of k. 

The theory in finite volume is described along the same lines. There are however two 
differences: 

• The total momentum gets quantized. 

• ipi^) (and therefore k) obtains finite volume corrections. 

This picture also applies to relativistic integrable theories. We consider Ac as a simple 
quantum mechanical bound state of two elementary particles and use the infinite volume 
scattering data to describe the interaction between the constituents. To develop these 
ideas, let us consider the spectrum of the theory defined on a circle with circumference 
L. We state the identification 

lAM), ~ \AMM02))l (4.3) 

where the 6^1 2 are complex to describe a bound-state; this idea also appeared in [95 l I109| . 
Relation (14.31) can be regarded as the finite volume realization of (14. 2p . The total energy 
and momentum of the bound state have to be purely real, constraining the rapidities to 
take the form 

9i = 9 + iu, 62 = 6 -iu (4.4) 
where the dependence on L is suppressed. Energy and momentum are calculated as 



E = 2ma cos(m) cosh(6') p = 2ma cos(m) sinh(^) 



(4.5) 



4.1. ONE-PARTICLE STATES 



91 



The quantization condition for an n-particle state is given by the Bethe-Yang equations 

n 

e''^'^l[S,^,^{e^-d,) = l, j = l...n 

k=l 

To quantize the bound state in finite volume, an appropriate analytic continuation of 
the above equations with n = 2 can be applied. This procedure is justified by the same 
reasoning that leads to original Bethe-Yang equations: one assumes plane waves (with 
imaginary momenta) except for the localized interaction, which is described by the S- 
matrix of the infinite volume theory. Inserting (14. 4p and separating the real and imaginary 
parts 

^inia cos{u) smh{0)L ^-rria sm{u) cosh{0)L g (2iu) = 1 (4 6) 

^ima cos{u) smh{e)L ^rua sm(u) cosh(0)L (— 2?^) = 1 (4 7) 

Multiplying the two equations and making use of S{2iu) = S{—2iu)^^ one arrives at 

27r/ 

g2im,cosHsinh(e)L _ ^ or 2m„ cos(m) sinh(0) = — (4.8) 

Ij 

which is the quantization condition for the total momentum. I is to be identified with 
the momentum quantum number of Ac. The quantization condition for u is found by 
eliminating 9 from (14. 6p : 



V Saa{2iu) = {-ly 



(4.9) 



The exponential factor forces u to be close to the pole of the S-matrix associated to the 
formation of the bound-state. For the case at hand it reads 

Saaie ~ Z<J ~ (4.10) 

with = 2m^^. Note the appearance of (—1)^ on the rhs. of (14.91) . which is a natural 
consequence of the quantization of the total momentum. This sign determines the direction 
from which the pole is approached. 

The exact solution of (14.91) can be developed into a power series in e"'^''"^, where the 
first term is found by replacing u with in the exponent: 

u-ul = {-IYI {T:S e-'^S.Wi+(^)^ + Oie-'^'^'^'^) (4.11) 

First order corrections to the energy are readily evaluated to give 

E = i?o-(-l)^(rL)^^e-^VO(e--»«-) (4.12) 
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where Eq is the ordinary one-particle energy 



Eo = \ rnl + 



271 1 ^ ^ 



L 



In the case of zero momentum the former result simplifies to the leading term in p.lip . 
For large volumes we recover 

2ttI 

u —>■ 6 —>■ arsh- 



ac 



Having established the quantization procedure we now turn to the question of momen- 
tum quantum numbers inside the bound state. For the phase shift we adopt the convention 
that was used throughout this work 

Note that 6ab{iu) is purely imaginary. 

With this choice of the phase shift the Bethe-Yang equations in their logarithmic form 

I sinh(6' + iu) + 5ii (2m) = 27rJi 
Zsinh(6' — iu) + 2m) = 27r/2 

imply Ji = I2. Quantization of the total momentum on the other hand requires Ii = 
I2 = 1/2. Note that different conventions for 6ab would result in a less transparent rule 
for dividing / among the two constituents. The only disadvantage of our choice is the 
appearance of the unphysical half-integer quantum numbers. With this convention the 
bound-state quantization can be written in short-hand notation as 



IU}Ll~IU/2,//2}) 



aa,L 



4.1.1 Numerical analysis 

It seems plausible that by exactly solving the bound state quantization condition (14.91) 
one obtains all higher order corrections that go as e~'^^'°-°-^ with n e N. In this subsection 
we present numerical evidence to support this claim. 

We investigate the Ising model in the presence of a magnetic field. The first three 
particles lie below the two-particle threshold and they all show up as AiAi bound states. 
These fusions are responsible for the leading /i-term. The corresponding parameters (in 
units of mi) are listed in the table below. The exponent of the next-to-leading correction 
(the error exponent) is denoted by 
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a 


rria 


a 

mi 






1 


1 


0.86603 


205.14 


1 


(mi) 


2 


1.6180 


0.58779 


120.80 


0.95106 




3 


1.9890 


0.10453 


1.0819 


0.20906 


(Ml) 



In [64] Klassen and Melzer performed the numerical analysis of mass corrections. The 
analytic predictions were compared to TCSA data and to transfer matrix results. They 
observed the expected behavior of mass corrections of Ai and A2; in the former case 
they were also able to verify the F-term. On the other hand, the precision of their TCSA 
data was not sufficient to reach volumes where the /z-term for ^3 could have been tested. 
This limitation is a natural consequence of the unusually small exponent fifi. the next- 
to-leading contribution is of order e"^'^"^, still very slowly decaying. 

We employ the TFCSA routines that were successfully used in the previous chapter. 
Calculations are performed for 1 = 0,1,2,3 at different values of the volume. One-particle 
states of Ai, A2 and A3 are easily identified: they are the lowest lying levels in the spec- 
trum, except for / = where the lowest state is the vacuum. We use the dimensionless 
quantities / = rriiL and e = E/rrii. 

The results are extrapolated from Ccut = 20.. 30 to Ccut = cxd. Experience from the 
previous results shows that this extrapolation technique reduces the numerical errors by 
an order of magnitude. We resign here from quantitatively monitoring the TCSA errors 
and constrain ourselves to a range of the volume parameter where it is safe to neglect 
truncation effects. 

^3 

We begin our analysis with the most interesting case of A3. At each value of / and I 
the following procedure is performed. 

• The energy correction is calculated according to (I4.12p 

• The quantization condition (14. 9p is solved for u and the energy correction is calcu- 
lated by 

Ae = 2 cosh(6') cos(m) — eo 

where 6 is determined by the total momentum quantization (14. 8p and eo is the 
ordinary one-particle energy. 

• The exact correction is calculated numerically by Ae = e^*"'^^ — eo. 

The choice for the range of the volumes is limited in two ways. On one hand, I has to 
be sufficiently large in order to reduce the contribution of the F terms and other higher 
order finite size corrections. On the other hand, numerical errors grow with the volume 
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Figure 4.1: Finite size corrections to one-particle levels (in sectors / = . . . 3) as a 
function of the volume. The TCSA data are plotted against theoretical predictions of the 
single /i-term associated to the AiAi —>■ A^ fusion (solid curve) and the exact solution of 
the bound-state quantization (dotted curve). 



and eventually become comparable with the finite size corrections, resulting in an upper 
bound on /. The window I = 30.. 40 is suitable for our purposes. 

The results are shown in figure 14. 1[ It is clear that the //-term yields the correct 
prediction in the L — ^ oo limit. However, higher order terms cause a significant deviation 
for / < 40, which is in turn accurately described by the bound state prediction. The sign 
of the correction depends on the parity of I as predicted by (14.121) . 

Based on the success of this first numerical test we also explored the region / < 30. 
Inspecting the behavior of m as a function of / reveals an interesting phenomenon. It is 
obvious from (14.111) that u{l) is monotonously increasing if / is odd, with the infinite 
volume limit fixed to m"^. However, the complex conjugate pair ^12 approaches the real 
axis as / is decreased and they collide at a critical volume / = Ic- For / < Ic they separate 
again but stay on the real line, providing a unique solution with two distinct purely real 
rapidities. The same behavior was also observed in [95 l I109| . 

The interpretation of this phenomenon is evident: if the volume is comparable to 
the characteristic size of the bound-state, there is enough energy in the system for the 
constituents to become unbound. Therefore the A3 one-particle level becomes an AiAi 
scattering state for / < Ic. We call this phenomenon the "dissociation of the bound state". 
The same result was obtained also in the boundary sine-Gordon model by a semiclassical 
analysis |110) . 
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Figure 4.2: A3 one-particle levels in sectors / = ... 3 as a function of the volume. 
Dots represent TCSA data, while the lines show the corresponding prediction of the AiAi 
bound state quantization. In sectors 1 = 1 and / = 3 the bound state dissociates at Ic and 
for I < Ic a conventional AiAi scattering state replaces ^3 in the spectrum. The values of 
Ic are shown by the two arrows. 



The value of Ic can be found by exploiting the fact that the Jacobian of the Bethe- Yang 
equations (viewed as a mapping from (6*1, 6*2) to (Ji, J2)) vanishes at the critical point. A 
straightforward calculation yields 

where (pii{0) = 6'ii{9). The numerical values for the case at hand are 

Ic = 27.887 (/ = 1) and Ic = 26.434 (/ = 3) 

We are now in the position to complete the numerical analysis. The Bethe-equations 
are solved at each value of /, providing two distinct real rapidities for / < Ic (with I being 
odd), and a complex conjugate pair otherwise. The energy is calculated in either 

e = cosh(^^i) + cosh(6'2) 

which is compared to TCSA data. The results are exhibited in figure [42l 

The agreement for the upper two curves (/ = 1 and / = 3) is not as surprising as 
it may seem because what one sees here are conventional AiAi scattering states. The 
Bethe- Yang equation determining their energy is exact up to 0{e~^'^) where /i' = /xj^^ is 
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the smallest exponent that occurs in the sequence of finite volume corrections of Ai. On 
the other hand, the energy levels are analytic functions of L, which leads to the conclusion 
that the prediction of (I4.9P is correct up to 0(e~''ii^) even for L > Lc. Comparing the 
numerical values one finds > We conclude that the bound state picture indeed 
accounts for finite volume corrections up to the first few orders in e~f^°^^ (the first 8 orders 
in the case at hand). 

Ai and A2 

Particles Ai and A2 also appear as AiAi bound states. However, there is no point 
in applying the complete bound state quantization to them, because the error terms 
dominate over the higher order contributions from (14. 9^ : the exponents of the sub-leading 
finite size corrections mi and /xf2 smaller than and Nevertheless, the leading 
/i-term can be verified by choosing suitable windows in /. 

In figures [47fl and [481 log( | Ae | ) is plotted against the prediction of (14.12P for / = 6.. 16 
and / = 6.. 22. (the sign of Ae was found to be in accordance with (14.12^ for both Ai and 
A2) 

In the case of Ai perfect agreement is observed for I = 10.. 18 in the sectors / = 
and / = 1. For 1 = 2 and / = 3 the energy corrections become too small and therefore 
inaccessible to TCSA (note that the prediction for J = 3 is of order 10~^). 

In the case of A2 precise agreement is found for / = 14.. 22 in all four sectors. 

A, 

Here we present an interesting calculation that determines the leading mass corrections 
of A^. The standard formulas are inapplicable in this case, because lies above the two- 
particle threshold. However, it is instructive to consider the composition of A^ under the 
bootstrap principle and to evaluate the /z-term prediction. 

There are two relevant fusions 

A1A3 A5 with /i?3 = 0.2079 

A2A2 -> A5 with /i^2 = 0.6581 

Numerical evaluation of Liischer's formula for the /i-term shows that the contribution 
of the second fusion is negligible for / > 30. The first fusion on the other hand yields a 
significant discrepancy when compared to TCSA data. This failure is connected to the 
two-particle threshold and it can be explained in terms of the bound state quantization. 
Experience with A^ suggests that one should first take into account the energy corrections 
of A3 and consider the A1A3 — > A5 fusion afterwards. A3 can be split into AiAi leading 
to the "triple bound state" Ai^iAi —* A^. In infinite volume one has (see also fig. II. li e.) 
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Figure 4.3: Finite size mass of as a function of the volume. The squares represent the 
TCSA data which is compared to the leading /x-term (solid curve) and the solution of the 
quantization condition for the "triple bound state" AiAiAi (dotted curve). The straight 
line shows the infinite volume mass. 

The finite volume realization of this identification is most easily carried out in the / = 
sector with 



The second equation is automatically satisfied due to unitarity and real analyticity, 
whereas the first and the third are equivalent and they serve as a quantization condi- 
tion for u. The finite volume mass of A^ is given in terms of the solution by 



In the large L limit the infinite volume mass is reproduced by m — 2ul^. Figure 14.31 
demonstrates the agreement between TCSA and the prediction of (I4.13p . 

The possibility of solving the quantization of the triple bound state in a moving frame 
looks very appealing. In the general case the rapidities are expected to take the form 
(^1 + iu, 02, Oi — iu) where 6i and 62 do not necessarily coincide. However, the numerical 
precision of our TCSA data was not sufficient to check our predictions. 

4.1.2 Comparison with the TBA for excited states 

As a conclusion of this section ( 14. lip is compared to the lowest order results of the TBA 
approach. The general discussion of excited states TBA equations in diagonal scattering 




e-"^^^^"(")^5n(zu)5n(2zu) = 1 
Sii{iu)Sii{—iu) = 1 
e'"i^'"(")^5n(-m)5n(-2m) = 1 



m^^l) = 2 cos(m) + 1 



(4.13) 
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and 



theories is not available. For simplicity we restrict ourselves to the Lee- Yang model which 
was considered in the original papers [25l [26] . 

The exact TBA equations for moving one-particle states read 

/oo ^0 
— m cosh (6') L (6') (4.14) 
-oo 2vr 

8(9) = niR cosh e + \og ^}l^~^'\ -{^^L){e) (4.15) 

where 

sinh(g) + 2 sin(7r/3) .dS{e) 
^^^^ = sinh(^^)-.sin(vr/3) ^^^^ = 

/oo inr 
^f{0 - 0')9{O') 

Here the volume is denoted by R to avoid confusion with L{9). The complex rapidity 9o 
satisfies the consistency equation 

e{eo) = mRcosheo + i-rc - log(5(2ilm^o)) - {'P^L){e) = i{2n + l)n (4.16) 

The convolution term in (I4.16p can be neglected and one obtains Im^^o = vr/6 + 5 where 
6 is exponentially small. To zeroth order one also has 

mi? cosh (Re6'o) = (4?t, + (1 — sign5))7r 

which is the ordinary one-particle quantization condition with / = 2n + |(1 — sign(5). 
Neglecting the contribution of the integral in f l4.14p and substituting 9q = 9 + iu one has 
E = 2msin(M) cosh(6'). This is exactly the energy of an AA bound state with the imaginary 
rapidities 9±iu. Separating the real and imaginary parts of (I4.16P and still neglecting the 
convolution term (which is responsible for the F-term) one obtains equations (14.91) and 
f l4.8p . thus proving the consistency of the two approaches. 



4.2 Multi-particle states 

4.2.1 Bethe-Yang quantization in the bound state picture 

Let us consider a scattering state |{/, h, . . . , In})cbi fon l composed of n + 1 particles, 
the first one being Ac. The energy of this state is calculated as 

n 

E = J2^^,cosh{9j) + ... (4.17) 

where {9,9i, . . . , 9n) is the solution of the Bethe-Yang equations 

Qj{9i, ...,9n)= sinh(%)L + ^ 5i^i,{9, - 9k) = 2nl, j = l...n (4.18) 
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We determine the leading part of the /x-term associated to (I4.17P by considering Ac as an 
AaAa bound state inside the multi-particle state. Therefore we write 

|{/, A, . . . , In})ctM - l{^/2> ^2, Ji, . . . , In}),a,M (4-19) 

The energy is then determined by analytic continuation of the n + 2 particle Bethe-Yang 
equations. They read 

n 

g-m,sm(«)cosh(e)Lgim,cosHsmh(e)L^^^^2m) Jj5„b^,(^ + m-0j) = 1 (4.20) 

J=l 

n 

i=i 

n 

e''^'^'^''''^'^^''Sab,{ej~9-iu)Sab^{e,-e + iu)\[Sb^,Se,-ek) = i (4.22) 

k=l 

The ordinary n + 1 particle Bethe-equations are reproduced in the L ^ oo limit by 
multiplying (14.20p and (14.2ip and making use of the bootstrap equation 

We now proceed similar to the previous section and derive a formula for the leading 
correction. The shift in the imaginary part of the rapidity can be calculated by making 
use of (ICTD and fliTnll as 

Au = u- ul^ = LJ^!^Le-^'^^^H(>)Le'"'^ smh(e)L/2 "Q g^^^ ^ --a^ _ Q,^ ^423) 

i=i 

Multiplying KW\ and 

n 

n 

e''^'^''''''^'^^''Sab^{e,-e-iu)Sa,,{e,-e + iu)\[s,^,^{e^-e,) = i (4.25) 

k=l 

Let us define 

Sab,{0 -lu- e,)Sab,{e + iu- e,) ^ Scb,{e - 0.)e'^-^^'^^'~'^^ 

where 

^cb, (0) = lipcb, {9 + lul,) - iipcb, {9 - iulc) with ipab{9) = 5'^^{9) 
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Using 2ma cos(m) 



rrir. 



the logarithm of (I4.24p and (14.251) can be written as 



Qo{0, 6*1, ... , 9n) 
Qj{9, 6i,. . . , On) 
The Ihs. can be expanded around the n + 1 particle solution 



2/i^„sinh(^)L- J] Vcb, 
^cb,{d - Oj)Au 



9i) Au 



,6n) to arrive at 



/ 9-9\ 
6,-6, 



/2/i^,sinh(elL-E"=i 



\ 



where 



dQk 
06, 



6„ 



Au 



(4.26) 



The final result for the energy correction reads 



AE = -2fxl^cosh{6)Au+ 



( mcSinh(^^) ^ 
sinh(^i^ 

ym6„ sinh(^„) / 



/2 



^cb,{6-6,) 



\ 



Vcbr, 



6„ 



- %A 



(4.27) 



Am 



with Am given by (14.230 . 

Based on the previous section it is expected that there is a similar contribution for 
every fusion leading to each one of the constituents of the multi-particle state. 

We wish to remark that in the case of the Lee- Yang model it is possible to derive 
(I4.27P from the excited state TBA, along the lines of subsection I4.1.2[ This carries over 
to other models with diagonal scattering, provided that the excited state TBA equations 
are available. The point of our calculation is, that the /i-term can be obtained without 
any reference to the TBA equations. Moreover, our arguments can be applied directly to 
form factors as well (see section 14. 3p . 



4.2.2 Multi-particle states - Numerical analysis 

We first consider finite size corrections to AiA^, states. They are not the lowest lying 
two-particle states in the spectrum, but they possess the largest /x-term which is connected 
to the AiAi — i> A3 fusion. Given a particular state \{Ii,h})i'iL ^^e following procedure 
is performed at each value of the volume: 

• The two-particle Bethe-Yang equation for \{Ii, h}) l solved and the /i-term is 
calculated according to ( 14.27p . 
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• The exact three-particle Bethe-Yang equation is solved for 73/2, 13/2})^^^^ ^ 

The results for different AiA^ levels are shown in figure I4.9[ The situation is similar to 
the case of the A3 one-particle levels: the bound state quantization yields a remarkably 
accurate prediction, whereas the single /i-term prediction only becomes correct in the 
L ^ 00 limit. 

In table 14.11 we present a numerical example for the dissociation of the bound state 
inside the two-particle state. In this case an A1A3 state turns into a conventional AiAiAi 
three-particle state at Ic ~ 30. 

Finite size corrections to AiAi and A1A2 states are also investigated, the leading fx- 
term given by the fusions AiAi — > Ai and AiAi A2, respectively. In the former case 
we calculate separately the contribution associated to both Ai particles and add them to 
get the total correction. Results are exhibited in figures [4T0l and [4.111 and formula (14.270 
is verified in both cases. 



4.3 Finite volume form factors 

The connection between finite volume and infinite volume form factors was derived in 
section 12.21 as 



,L{{i[,...,iL}\o{o,o)\{h,...,in}),,. 



F^je'^ + in,..., e[+ 1^,61,..., en)j^...j,i,...i„ ^ 0(e-^'^) (4 28) 

A/Pn...j„(6'l, . . . ,9n)Pj^...j^{9[, ...,6'^) 

where the rapidities 9 are solutions of the corresponding Bethe-Yang equations. Here we 
assume for simplicity, that there are no disconnected terms present, ie. it is supposed that 
9j 7^ 9'i^ whenever ij = i^. 

Based on general arguments it was shown in subsection 12.2.31 that /i' > /i where /x is 
determined by the pole of the S-matrix closest to the physical line. A systematic finite 
volume perturbation theory (Liischer's method applied to form factors) is not available. 
However, it is expected that the actual value of /x' depends on what diagrams contribute 
to the form factor in question. Apart from the insertion of the local operator they coincide 
with the diagrams determining the finite size corrections of the multi-particle state. There- 
fore ii' is associated to the bound state structure of the constituents of the multi-particle 
state. 

In this section we show that the leading correction term can be obtained by the bound 
state quantization. 
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4.3.1 Elementary one-particle form factors 

(14.280 yields a simple prediction for the elementary one-particle form factor: 

Ff(/,L) = (0|O(0,0)|{/}),^^ = -^ + 0(e-'^'^) (4.29) 

where E is the one-particle energy, and F^ = F^{9) is the infinite volume one-particle 
form factor, which is constant by Lorentz symmetry. 

The /i-term associated to (14.29^ is derived by employing the bound state quantization. 
We gain some intuition from the section 14.11 where it was found that the bound state 
AaAa may dissociate at a critical volume Lc. For L < Lc there is no one-particle level of 
type Ac in the given sector of the spectrum, however an AaAa scattering state appears 
instead. Finite volume form factors of this state are calculated using (|4.28l) as 

i^f(/,^) = ^^;^%% forL<L. (4.30) 

The generalization to L > Lc seems to be straightforward: one has to continue an- 
alytically (14.301) to the solutions of the Bethe-Yang equation with imaginary rapidities 
9i 2 = 9 ± iu. However, note that equations (14.29P and (14.301) are valid up to a phase 
factor. In order to continue analytically to imaginary rapidities we also need to fix this 
phas^. 

The two-particle form factor satisfies 

F'^{9i-, 6*2)00 = Saa{9i — 92)F'^{92, 9i)aa 

The simplest choice for the phase is therefore 

F''{9u92)aa = VSaa{9i-92)\F'^{9u92)aa\ (4.31) 

This choice is dictated by CPT symmetry [1111 [52l lll2| . and it is respected by all known 
solutions of the form factor bootstrap axioms. There is a sign ambiguity caused by the 
square root, but it can be fixed by demanding (S'aa(O))^/^ = i and continuity. Using ( 14.3ip 



,OrT r^ _ VSaai92~9i)F^{9i,92) 



and upon analytic continuation 



,0 ^Saa{-2lu)F<^{9 + lU,9-lu) 



F^{I, L) = ^ ' \ ' for L>Lc (4.32) 

\/ Paa{9 + IU,9 - lu) 



^As we noted in earlier the phase of a (non-diagonal) infinite volume form factor is unphysical in 
the sense that it may be redefined by a complex rotation of the state vectors and physical quantities, 
e.g. correlation functions, do not depend on such redefinitions. However, the bootstrap program uniquely 
assigns a phase to each form factor. 
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It is easy to see that the result (14.291) is reproduced in the L ^ oo limit. First observe 
that 

pO^e + ZU,e- Zu)aa ~ . Ff (g) 

The residue of paa is determined by (paai^iu) and it reads 

S" (2iu) 1 
+ iu,9 — iu)aa ~ '^'TfiaL cos(m) cosh(6')(— i) °° , . = cosh(6')- 



The singularities in the numerator and denominator of ( 14. 32^ cancel and indeed 



acJ 



^JnricL cosh(6') 



We emphasize that it is crucial to include the extra normalization factor ^J Saai—'^iu) to 
obtain a meaningful result. 

Expression (14.321) can be developed into a Taylor-series in m — First we use the 
exchange axiom to arrive at 



,OfT T^ _ F^{e -iu,9 + iu) 



F^{I, L) = ' ' ' = (4.33) 

y/Saa{-'2lu)paaid + lU, 9 - lu) 



The form factor axioms imply that 



hm F''{9-iu,e + iu)aa = ^F^ 

aa 



The simple pole of ipaai'^iu) in paa is canceled by Saa{—'2iu), therefore both the numerator 
and the denominator of (14.331) have continuous limits as m — m^^. 

The form factor F'-'{6—iu, 9+iu)aa only depends on u by Lorentz-symmetry. Therefore 



aa 

where 



Expanding the S-matrix element into a Laurent-series in the vicinity of the pole 

Saa{2tU) = +^:;,° + ... 

2(m - Kc) 



Expanding the denominator: 

S{—2iu)paa{6 + iu,9 — iu) = 

S{-2iu)EiE2L^ - i{Ei + E2)LS' {2iu){s{-2iu)y = 
EcL f2EiE2L^ AEcL 



+ ( ~^r:: - -^p^s'-' iu - ud + ... 



ere Y y Y (Y 

\ aal ^ \ aa) \ 



^ .4: aa I \ ac) 
aa) 
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where 



Ec = El + E2 = 2ma cos(m) cosh(6') 



Putting all this together 

F^{I,L) = ^^ + 



+ 



o 



+ 



-E,E,L' + 



(re 

V ac 



.2 aa 



(u - + 



Note that in the preceding formulas Ec does include the leading order correction to the 
usual one-particle energy E^ = a/wi^ + {2txIY/L'^. Using 



the final result is given by 



-<J+0(e-2^^) and EiE2 = ^,El-^^ 



u 



mt. 



F^{I L) — 2«r^^ {Faa) , 



u - Kc) + 



■^0 



(r 



aaJ 



(El 



o^2 



m 

m2"^ E^ 



-C)+0(e-2^^) (4.34) 



with 



4.3.2 One-particle form factors - Numerical analysis 

Let us introduce the dimensionless form factors as 

^Ok(0,0)|{/}),, 



nil 



We use the methods described in section [3TT] to determine fi{I,l) for i = 1,2,3 and 
/ = 0, 1, 2, 3. The numerical results are compared to the exact infinite volume form factors. 

We start our investigation with fs{I,l), for which relatively large exponential correc- 
tions were found in subsection 13.2. 2[ It is convenient to consider 



Ml,l) = {eolf'Ml,L) with lim/3(/,0 = i^; 



1^00 



(4.35) 



The numerical results are demonstrated in fig. 14. 4[ Note, that this is the same figure as 
fig- 13.71 but here the interpretation of the huge deviations from F3 is also provided. 

We also tried to verify the predictions for fi and f2- In the latter case reasonably good 
agreement was found with TCSA, the results are demonstrated in fig. I4.12[ In the case of 
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Figure 4.4: Elementary finite volume form factors of one-particle levels. Here the 
normalization (14.35P is applied to obtain a finite / oo limit, which is given by the 
infinite volume form factor F3 = {0\e\A3{9)). The TCSA data are plotted against the 
bound state prediction. The ordinary evaluation of /a is simply the constant F3. 



/i we encountered the unpleasant situation that the F-term decays slower than the TCSA 
errors grow, thus making the observation of the /i-term impossible. 

It is straightforward to generalize (14.321) to matrix elements between two different 
one-particle states. For b ^ c one has for example 

/rrii im\ F^{9b + i'K,9 + iu,9 -iu)baa 

\/ Pb[9b)Paa[9 + IU,9 - lU) 

Numerical examples are presented in figures [4751 (a)-(c) for c = 3 and b = 1,2. 

The most interesting case is the one shown in fig. 14.51 (d) where the matrix element 
between two different A3 one-particle states are investigated. This can be done by con- 
sidering both ^3 particles as the appropriate AiAi bound states and then calculating the 
finite volume form factor 3{{I}\s\{I'})3 1 as 

n{{I/2,I/2}\e\{I'/2,I'/2}),,,= ^0 + ^u + ^.,9 - ^u + ^n,9' + ^u' ,9' - ^u'U, 
/ ' / ^1 11 / ' / ^/ii.i ^pn(^^' + tW, 9' - iu')pii{9 + iu,9- lu) 



Once again we find complete agreement with the TCSA data. 
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Figure 4.5: One-particle-one-particle form factors, dots correspond to TCSA data. The 
solid lines represent the ordinary evaluation of the finite volume form factors, while the 
dotted lines show the bound state prediction. 

4.3.3 Elementary multi-particle form factors 



The generalization of (14.321) to multi-particle states is straightforward, the only task 
is to find the appropriate phase factor. Similar to the one-particle case one has 



i<j 



A general n particle finite volume form factor with real rapidities can thus be written as 

-F^{9i, . . . , 9n)bi...b„ 



p"(^l,...,ft„),,...b„ 

Substituting the solution of the Bethe-equation for the state |{//2, J/2, h, . . . , In})aabi b„ l 
and making use of the real analyticity condition 



\S. 



ab.yOj 



9-iu)SabA9,~9 + zu)\ = 1 



one gets 



(0|(9|{//2,//2,/i,...,/„})^^,^...,^^^ 



A/5'aa(-2m) \F'^{9 + iu,9 -iu,9i,..., 9n)aab^...bn I 



^p(n+2)(^ + m, ^ - m, ^1, . . . , 9n)aab,. 



(4.36) 
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(a) (0|e|{l,l})i3_^ 




20 22 24 26 28 30 32 
/ 

(b) (0|£|{2,0})i3^^ 



Figure 4.6: Elementary form factors of A1A3 scattering states, dots correspond to TCSA 
data. The solid lines are obtained by a "naive" evaluation of the finite volume form factors, 
while the dotted line represents the bound state prediction, (in this case Ai^iAi form 
factors at the appropriate rapidities) 



up to a physically irrelevant phase. 

It is easy to show once again that the "naive" result is reproduced in the L ^ 00 limit. 
To do so, we first quote the dynamical pole equation of the infinite volume form factor: 

F^{e + iu,d- tu, e,,..., en)aab,...K = ^y^. F^'ie, e,,..., e^),t,...K + o{i) 

2(m - Kc) 

The singularity of is given by 

ReS,^S.^p("+'H^ + iu,9- iu, 9i,..., 9n)aab^...b„ = ^P^"+'^ iO,Ol,..., 9n)cb,...b„ 

The "naive" formula is now recovered by inserting the last two equations into (14.361) . 

The leading exponential corrections can be obtained by plugging (14.23P and (14.260 
into (14.361) and expanding to first order in m — u^^. This procedure is straightforward but 
quite lengthy, therefore we refrain from giving the details of the calculations. 

In fig. 14.61 two examples are presented for the evaluation of (14.361) applied to AiA^ 
two-particle states. 
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Figure 4.7: Finite size corrections to Ai one-particle levels in sectors / = . . . 3, logioAe 
is plotted as a function of the volume. Dots represent TCSA data, while the lines show 
the yU-term corresponding to the AiAi —>■ Ai fusion. 
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Figure 4.8: Finite size corrections to A2 one-particle levels in sectors / = . . . 
is plotted as a function of the volume. Dots represent TCSA data, while the 
the yU-term corresponding to the A^Ai — > A2 fusion. 



3, logioAe 
lines show 



4.3. FINITE VOLUME FORM FACTORS 



109 



0.03 
0.025 - 

0.02 
0.015 - 

0.01 
0.005 - 



22 



26 28 30 

(a) |{l,0})i3 



0.02 
0.018 
0.016 

0.014 
0.012 
0.01 
0.008 
0.006 
0.004 



(b) |{0,2}),3 




20 25 30 

(C) |{0,l})i3 



22 24 26 28 30 32 34 



(d) |{l,l})i3 



Figure 4.9: Finite size corrections to AiA^ scattering states as a function of the volume. 
Dots represent TCSA data, the sohd hue shows the ;U-term corresponding to the AiAi 
A3 fusion. The dotted lines are obtained by the exact solution of the quantization condition 
for the AiAiAi three-particle system. 
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Figure 4.10: Finite size corrections to AiAi scattering states, logioAe is plotted as a 
function of the volume. Dots represent TCSA data, while the solid line show the sum of 
the two /x-terms corresponding to the AiAi Ai fusions. 
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(a)|{0,l})i2 (b)|{l,l})i2 

Figure 4.11: Finite size corrections to ^1^42 scattering states, logioAe is plotted as a 
function of the volume. Dots represent TCSA data, while the solid line show the /x-term 
corresponding to the A^Ai —>■ A2 fusion. 



I 


Oi 


02 




e{l) (predicted) 


e{l) (TCSA) 


22 


0.44191 


0.68235 


-0.58742 


3.51877 


3.51900 


23 


0.42574 


0.64280 


-0.55190 


3.46201 


3.46219 


24 


0.60523 


0.41155 


-0.51892 


3.41239 


3.41255 


25 


0.56934 


0.39931 


-0.48831 


3.36890 


3.36907 


26 


0.38910 


0.53475 


-0.45989 


3.33071 


3.33089 


27 


0.50092 


0.38119 


-0.43350 


3.29709 


3.29729 


28 


0.46686 


0.37634 


-0.40899 


3.26743 


3.26767 


29 


0.42947 


0.37739 


-0.38622 


3.24122 


3.24164 


30 


0.38646 + 0.02332 i 


0.38646 - 0.02332 i 


-0.36505 


3.21801 


3.21832 


31 


0.37059 + 0.04042 i 


0.37059 - 0.04042 i 


-0.34537 


3.19741 


3.19775 


32 


0.35572 + 0.05104 i 


0.35575 - 0.05104 i 


-0.32706 


3.17908 


3.17945 


33 


0.34182 + 0.05891 i 


0.34182 - 0.05891 i 


-0.31002 


3.16275 


3.16313 


34 


0.32876 + 0.06513 i 


0.32876 - 0.06513 i 


-0.29415 


3.14816 


3.14856 


35 


0.31650 + 0.07022 i 


0.31650 - 0.07022 i 


-0.27936 


3.13511 


3.13547 


36 


0.30498 + 0.07446 i 


0.30498 - 0.07446 i 


-0.26556 


3.12341 


3.12235 



Table 4.1: An example for the dissociation of the AiAi bound state inside a scattering 
state. |{2, ODgj^ ^ is identified with |{1, 1, 0})j^^^ ^ and the corresponding Bethe-Yang equa- 
tions is solved. For / < 30 there is a real A\A\A\ three-particle state in the spectrum, 
whereas at / ~ 30 two of the rapidities become complex and the two-particle state A^Aj, 
emerges. 
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Figure 4.12: Finite size corrections to the elementary form factors of ^2- Dots represent 
TCSA data, while the lines show the //-term prediction corresponding to the AiAi — > A2 
fusion. 
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Chapter 5 



Correlation functions at finite 
temperature 

5.1 The form factor expansion for correlation functions 

The finite temperature correlation function of local operators Oi and O2 is defined by 

(o.(...,)o.(.....)>«^ ^(-""''^';^^°-'--'-" ,5.1) 

where R = 1/T is the temperature dependent extension of the Euclidean time direction 
used in thermal quantum field theory and H is the Hamiltonian. 

In the following we briefly discuss the evaluation of the correlation function at zero 
temperature using the form factor expansion and show that there is no straightforward 
generalization of this method to the finite temperature case. To keep the exposition simple 
we assume that the spectrum contains a single massive particle of mass m. 

At zero temperature the two-point function of operators Oi and O2 is given by the 
vacuum expectation value 

(0|Oi(xi,ri)O2(x2,r2)|0) 
Inserting a complete set of asymptotic states and using Euclidean invariance one obtains 

{0\O,ixuh)O2ix2,h)\0) = 

00 

/ d9,...d9r, F^-{e^...e^)F^-{e^...e^r e— '^^r^i-^M^^) (5.2) 
n=0 ^' ^ 

where r = ■>/ {xi — 0:2)^ + (ti — r2)^. The essence of the form factor approach [HHl [TBI 11131 
11141 [To] is to determine the form factors appearing in (15.21) as solutions of the form factor 
bootstrap program, and then numerically integrate the resulting expressions to obtain the 
contributions to the correlator. Due to the exponential damping factor and the vanishing 
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property of the form factor functions at = the series (I5.2p is rapidly converging even 
for relatively small values of mr, therefore in most cases it is sufficient to consider only 
the first few terms in (15.2p . 

Switching on a non-zero temperature one faces severe problems. A straightforward 
application of the form factor expansion results in 

(Oi(0,0)O2(0,r))«= ^^(m|Ci|n)(n|02|m)e-(^-")^*e-"^*' (5.3) 

m,n 

where for simplicity we have chosen the displacement between the two operators to lie in 
the imaginary time direction and the summation over m and n runs over the complete 
set of asymptotic states of the theory. Expression (15.31) is ill-defined in two ways: 

• The partition function Z is infinite if one considers the infinite volume theory. 

• One has to integrate over the poles of the form factors associated to the disconnected 
contributions, whenever the rapidities of two particles entering m and n approach 
each other. 

One encounters similar problems already by the evaluation of the thermal expectation 
value of local operators. The formal form factor expansion for the one-point function 
reads 

(O)" - '^'XZh) = I E I'nlOHe-''-' (6.4) 

The diagonal form factors entering the expression above are ill-defined in the infinite 
volume theory due to disconnected terms. 

In order to obtain well-defined expressions for the correlation functions it is necessary 
to introduce a regularization scheme. 

A. LeClair and G. Mussardo proposed an evaluation scheme [30] based on the Ther- 
modynamic Bethe Ansatz. It was proven by Saleur [H] that it yields the correct results 
for expectation values of local conserved charges. However it was also argued in [3l] that 
the LeClair-Mussardo formalism fails in the case of the two-point function. The same con- 
clusion was reached in [35] studying massless limits of correlation functions. G. Delfino 
proposed a second, independent regularization scheme [33] for the one-point function; 
these two approaches are presented and compared in subsection 15.1.11 Mussardo showed 
using a toy model, that the two methods yield different results [34]; the same conclusion 
is reached by comparing the systematic low-temperature expansion of the two proposals 
(see [5XB . 

A third approach is to consider a theory in finite volume L, to evaluate the form factor 
expansion using the well-defined finite volume form factors, and to take the limit L oo. 
This method is built on first principles and it does not rely on any additional assumptions; 
moreover it is expected to work for any n-point function. 
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In this work we develop a systematic low-temperature expansion for the one-point 
function with a well-defined L ^ oo limit; this calculation is presented in section [521 We 
find complete agreement with the LeClair-Mussardo proposal up to the third non-trivial 
order. In section 15.31 we apply our method to the two-point function. However, we only 
calculate the simplest terms; the systematic evaluation of the form factor expansion is out 
of the scope of the present work. 

Finally we note, that it is also possible to numerically evaluate the form-factor expan- 
sion [HH EHl SO] or to develop a term by term regularization scheme based on the infinite 
volume connected form factors |115) . Eventually, after the completion of this thesis there 
appeared an independent work by Konik and Essler |116| which uses (in addition to a 
novel infinite volume regularization procedure) finite volume techniques similar to the 
ones presented in this work. 



5.1.1 Comparison of the LeClair-Mussardo and the Delfino ap- 
proaches 

LeClair and Mussardo proposed the following expression for the low temperature (T <^ 
m, or equivalently mR ^ 1) expansion for the one-point function (15.41) : 



n=0 



n\ (27r) 



.i=l 



(5.5) 



where is the connected diagonal form factor defined in eqn. (I2.38P and e{9) is the 
pseudo-energy function, which is the solution of the thermodynamic Bethe Ansatz equa- 
tion 

/rif)' 
— ip{e -e')\og{l + e-''^^'^) (5.6) 

The solution of this equation can be found by successive iteration, which results in 



mRcoshiO)- [ ^^(0 - + i / ^^(9-9')e 

J 2n 2 J 2n 



/\—2mRcosh9' 



^{6 - e')ip{e' - e")e 



//N „—mR cosh d'—mR cosh d" 



27r 27r 

Using this expression, it is easy to derive the following expansion from (15. 5p 

dO 
2^ 
dOi d92 



+ 



(5.7) 



(O) 



R 



(O) 



rpc ( —mRcoshO —2mRcosh9\ 

-^2 e ) 



27T 27T 



iF^ie,,e2) + 2^e,-e: 



pc\ g— mi? cosh Si g— mi? cosh 02 



2K2 



+0 (e-=^™^) 



(5. 



where {O) denotes the zero-temperature vacuum expectation value. The above result can 
also be written in terms of the symmetric evaluation (I2.35p as 
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1^0^^ 1^0^ -|- j ^g— mi? cosh 6* g-2mi?cosh6'j _|_ 

2 y 27r 27r 

where we used relations ( 12.39^ and (I2.4ip . Obtaining the third order correction from the 
LeClair-Mussardo expansion is a somewhat lengthy, but elementary computation, which 
results in 

6 / 27 27 27^6^^^'^^'^^)" 



dOi de2 

27 27 

1 f deide2 

2 



— 3r7ii? cosh 9\ 



! ^^F^X^i - ^2)e-™^(™^'^'^+2coshe.) ^5_^Q^ 
J 2tt 2tc 



where we used eqns. (12.391 12.411 I2.42P to express the result in terms of the symmetric 
evaluation. 

Delfino's proposal for the one-point function reads: 



n=0 ^ ^ li=l 



Q—mR cosh 9i 



-mR cosh ( 



(5.11) 



which gives the following result when expanded to second order: 



1 j ^^Fl{e,, ^2)e-"^^(™^'^'i+™^'^'^) + O (e-^-^) (5.12) 



Note that the two proposals coincide with each other to this order, which was already 
noted in [33]. However, this is not the case in the next order. Expanding (15.111) results in 

1 f d6id62d6^ /I \ -mi?(coshei+cosh6»2+cosh6»3) 

gJ 27 2^27^^^^^'^^'^^^" 

f '^^1 '^^2 jj^sfn n N„-miJ{cosh 6*1+2 cosh 6*2) , f ps„-3m/? cosh 6»i -tO\ 

It can be seen that the two proposals differ at this order (the last term of (15.101) is missing 
from (I5.13P ). which was already noted by Mussardo using a toy model in |34|, but the 
computation presented above is model independent and shows the general form of the 
discrepancy. 
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5.2 Low-temperature expansion for the one-point func- 
tion 

We now evaluate the finite temperature expectations value in a finite, but large volume 

L: 

where Hi is the finite volume Hamiltonian, and Tr^ means that the trace is now taken 
over the finite volume Hilbert space. For later convenience we introduce a new notation: 

• • • , On)L = \{h, ■ ■ ■ , In})L 

where 9i,. . . ,6n solve the Bethe-Yang equations for n particles with quantum numbers 
Ii, In at the given volume L. We can develop the low temperature expansion of (15.140 
in powers of e"*"^ using 



^1 J2 'e-"^^(™^^^^''+™^'^^^'^^f\^f |(9|^f\ef)i + 
6 

n(3) n(3) ^(3) 

+0(e-^™^) (5.15) 



and 



Tr^ (e"^-^^) = 1 -[- y^g-m.ficosh(e(i)) _^ j_ 'g-mR(cosh(ef' )+cosh(4^')) 

^ /9(3) fl(3) fl{3) 

The denominator of (15.141) can then be easily expanded: 



2 



_ g-mi?coshe(i) _|_ I g-m_Rcosh6l(i) I _ }_ 'g-m_R(cosh '+cosh ) 



Tri(e-^^^) ^ \^ 2 

^Hcosh6l(l) ^ , [ -mRcoshefi) \ ST^ ' -mR(coshef' +cosh6l^^') 



g-mH cosh 61(1) j _|_ j g-mRcoshefi) j ' ^ 

_^ 'g-mH(cosh6»f '+coshe^^^+cosh6»^^') _|_ -4m_R\ (5 17) 



6 

/,{3) fl(3) p(3) 
''^2 '^3 
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The primes in the multi-particle sums serve as a reminder that there exist only states for 
which all quantum numbers are distinct. Since we assumed that there is a single particle 
species, this means that terms in which any two of the rapidities coincide are excluded. All 
n-particle terms in (15.151) and (15.161) have a 1/n! pre- factor which takes into account that 
different ordering of the same rapidities give the same state; as the expansion contains 
only diagonal matrix elements, phases resulting from reordering the particles cancel. The 
upper indices of the rapidity variables indicate the number of particles in the original finite 
volume states; this is going to be handy when replacing the discrete sums with integrals 
since it keeps track of which multi-particle state density is relevant. 

We also need an extension of the finite volume matrix elements to rapidities that are 
not necessarily solutions of the appropriate Bethe-Yang equations. The required analytic 
continuation is simply given by eqn. (12.361) 

(^1, . . . , 9M0l, • • • , On)L = ^ E ^2|A| ({^^. W)Pn-|A| ({^.}.^A)L+0(e-'^^) 

(5.18) 

where we made explicit the volume dependence of the n-particle density factors. The last 
term serves as a reminder that this prescription only defines the form factor to all orders 
in 1/L (i.e. up to residual finite size corrections), but this is sufficient to perform the 
computations in the sequel. 

Using the leading behavior of the n-particle state density, contributions from the n- 
particle sector scale as L", and for the series expansions ( I5.15p . (I5.16p and (15.171) it is 
necessary that mL -C e™^. However if mR is big enough there remains a large interval 

where the expansions are expected to be valid. After substituting these expansions into 
(I5.14P we will find order by order that the leading term of the net result is 0(L°), and the 
corrections scale as negative powers of L. Therefore in (15.141) we can continue analytically 
to large L and take the L ^ oo limit. 

It is an interesting question why this limit exists, or to be more precise, why do form 
factors with low number of particles determine the leading terms of the low-temperature 
expansion. As a matter of fact, the particle density N/L has a non- vanishing L ^ oo limit, 
the thermal average is therefore dominated by thermodynamic configurations where the 
particle number grows proportional with the volume. However, the leading contributions 
come from disconnected terms and they still involve form factors with a low number of 
particles, as it was pointed out in |115| . 
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5.2.1 Corrections of order e "^^ 

Substituting the appropriate terms from (15 .17^ and f lS.lSp into (15.141) gives the result 

{0)R = {0)l + ^e-"^'^"^^^^''' - {0)l) + 0(e-2'"«) 



em 



Taking the L ^ oo limit one can replace the summation with an integral over the states 
in the rapidity space: 



and using (12.321) we can write 



so we obtain 



p^{9)mO\9)L-{0)L)=F^ + 0{e->''') 

/df) 
27r 



(5.19) 



which coincides with eqn. (15. 9p to this order. 



5.2.2 Corrections of order e 



'2mR 



Substituting again the appropriate terms from (15.170 and (I5.15P into (I5.14p gives the 
result 



9(1) 

/ \ / 



-mR cosh 



{{e^?\o\e^^^)L- {0)l 



+1 ^ 'g-™i?(coshef)+cosh^^f ) (^(^ef\ef\o\ef\ef^)L-{o)L)+o{e-'''^'') 



5(2) ^g(2) 



The 0(e ^'"^) terms can be rearranged as follows. We add and subtract a term to remove 
the constraint from the two-particle sum: 



^1 ^ e-*^^^^i''+™^^^^^'')((^^f\ef 

g(2)^g(2) 
g(2)^g(2) 



e^^) 



120 



CHAPTER 5. CORRELATION FUNCTIONS AT FINITE TEMPERATURE 



(2) (2) 

The 9\ =02 terms correspond to insertion of some spurious two-particle states with 
equal Bethe quantum numbers for the two particles (Ji = J2). The two-particle Bethe- 
Yang equations in this case degenerates to the one-particle case (as discussed before, the 
matrix elements can be defined for these "states" without any problems since we have 
the analytic formula (15.181) valid to any order in l/L). This also means that the density 
relevant to the diagonal two-particle sum is pi and so for large L we can substitute the 
sums with the following integrals 

"1,2 ''l —''2 ''l '''2 

Let us express the finite volume matrix elements in terms of form factors using (12.321) and 
(12:3311 : 

P2{0i,e,)[{ef\ef\o\df\df^)L-{o)L) 

-Pi {Oi) Pi (^2) {{ei\o\e,)L + (^2|0|^2)l - 2(0)z.) = Fi{e,, e^) + o(e-'^^) 

Combining the above relation with ( I5.19p , we also have 

P2{0,9) 

where we used that Fl{9,9) = 0, which is just the exclusion property mention after eqn. 
(I235|). Note that 

"'""■^ -1+0(L-') 



P2(«,9) 

and therefore in the limit L ^ 00 we obtain 

/ ^^-2mR cosh ps -\- - f ^^^^^^(5)^^ ^^^g-m_R{cosh 6»i+cosh ^2) 
27r ^ 2 y 27r 27r ' 

which is equal to the relevant contributions in the LeClair-Mussardo expansion (|5.9I) . 

5.2.3 Corrections of order e~^"^^ 

In order to shorten the presentation, we introduce some further convenient notations: 

Ei = m cosh 6i 

{ei,...,en\o\ei,...,e^^)L = {i...n\o\i...n)L 

Pn{Oi, . . . ,9n) = p{l . . .n) 

Summations will be shortened to 

E - E 

6*1. ..6»„ l...n 

E' - E' 

e,...e„ i...n 
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Given these notations, we now multiply (15.151) with (I5.17P and collect the third order 
correction terms: 



1 y e-^(^i+^2+i?3) (/i23|C|123)l - 
6 

123 

V 1 J 23 




L 



To keep trace of the state densities, we avoid combining rapidity sums. Now we replace 
the constrained summations by free sums with the diagonal contributions subtracted: 

E' = E-E 

12 12 1=2 

E' = e-(e + e+e)+25: 

123 123 \1=2,3 2=3,1 1=3,2/ 1=2=3 

where the diagonal contributions are labeled to show which diagonal it sums over, but 
otherwise the given sum is free, e.g. 

E 

1=2,3 

shows a summation over all triplets 6^\6^\6^^ where = ^2^'' and ^3^^ runs free (it 
can also be equal with the other two). We also make use of the notation 

F{12...n) = F^„{9i,...,9n) 

so the necessary matrix elements can be written in the form 

= F(123) + p(l)F(23) + ■ ■ ■ + p(12)F(3) + . . . 

= 2p(2)F(12) + 2p(12)F(3) + p(22)F(l) 

= 3p(lll)F(l) 

= F(12)+p(l)F(2)+p(2)F(l) 

= 2p(l)F(l) 

= F(l) (5.20) 

where we used that F and p are entirely symmetric in all their arguments, and the ellipsis 
in the the first line denote two plus two terms of the same form, but with different 
partitioning of the rapidities, which can be obtained by cyclic permutation from those 
displayed. We also used the exclusion property mentioned after eqn. (I2.35p . 



p(123)((123|C>|123)l- 


{0)l 


p(122) ((122|C|122)l - 


{0)l 


p(lll)((lll|0|lll)i- 


{0)l 


p(12)((12|0|12)z.- 


{0)l 


p(ll)((ll|0|ll)z.- 


{0)l 


p(l)((l|0|l)^- 


{0)l 
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We can now proceed by collecting terms according to the number of free rapidity 
variables. The terms containing threefold summation are 



lj2^.R(^E,+E,+E,) ((i23|0|123)i - {0)l) - ^ E E ((23|0|23)z. - 

123 1 2,3 

V 1 2 3 1,2 3 / 



Replacing the sums with integrals 

- / 

1 



1,2 



y ^ /^^^p(i23) 



1,2,3 



and using (15.20p we get 

1 r ^^^g-i?(Si+s.+S3)(^(123) + 3p(l)F(23) + 3p(12)F(3)) 
6 J 2n 2n 2tc 

- 1 1 ^55e-^^^^^^^^^^Up(l)i^(23) + 2p(l)p(2)F(3)) 

r ^^^^_.(.,,.,,.3) (p(l)p(2)F(3) - ip(12)F(3) 
J zn zn zn \ z 



+ 



where we reshuffled some of the integration variables. Note that all terms cancel except 
the one containing F(123) and writing it back to its usual form we obtain 

6 J 2t[ 2tt 2tt 

It is also easy to deal with terms containing a single integral. The only term of this form 
is 

- J2 e"'^^'''^'''^'''H(123|C|123)L- 

^ 1=2=3 
(3) (3) (3) 

When all rapidities 0\ ,62 ,9^ are equal, the three-particle Bethe-Yang equations reduce 
to the one-particle case 

mL sinh Of^ = 2nh 
Therefore the relevant state density is that of the one-particle state: 



l|^e— V(l)((lll|0|lll>.-(0>.) - V(l)^ni) 
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where we used that 
when L —>■ oo. 

The calculation of double integral terms is much more involved. We need to consider 

^.R(^E,+E,+E,) ((i23|C)|i23)i - 
^.u(E,+j,,+E,) ((23|C|23)l - {0)l) 

1 2=3 

^.RiE,+E2+Es) ({3\0\3)l - {0)l) (5.23) 

1=2 3 

We need the density of partially degenerate three-particle states. The relevant Bethe-Yang 
equations are 

ml sinh 9i + 6{9i - 62) = 27r/i 
mL sinh 62 + 26(92 - Oi) = ^nh 

where we supposed that the first and the third particles are degenerate (i.e. I3 = Ji), 
and used a convention for the phase-shift and the quantum numbers where 5(0) = 0. The 
density of these degenerate states is then given by 




p(13,2) = det 



LE^ + if{ei-e2) -^(^1-^2) 
-2^(01 - 62) LE2 + 2ip{e^ - 62) 



where we used that ip{6) = ip{—9). Using the above result and substituting integrals for 
the sums, we can rewrite eqn. (I5.23P in the form 

1 r ^^^-RiE.+2E,)^^^^2)^F{2) 

2 J 271 271 '^^ ' p{22) ^ ' 
1 f dOidOs _j^(^2Ei+E3) .^n\ 1 



+ 2 J 27 27»-"'^""-"''W^(^'m3)^<=" 
where the ellipsis denote two terms that can be obtained by cyclical permutation of the 
indices 1, 2, 3 from the one that is explicitly displayed, and these three contributions can 
be shown to be equal to each other by relabeling the integration variables: 

-'2 I §^e-«^^^^^^^^^^(2p(l)i^(12) + 2p(12)F(l)+p(ll)F(2)) 

1 r ^^,-RiE.^2E,)^^^^2)^Fi2) 

2 J 271 271 ^'^^ >(22) ^ ^ 
1 f dOidOs _j^^2Ei+E3) 1 



+ 2 J 2^ 2^^"'"^"""'"«'"^)m3)^<^) '5-^*) 
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We first evaluate the terms containing F(23) which results in 

fd9id92 s/Q /] X -miJ(coshei+2cosh6»2) (K 0K\ 

~ J 2^27^^ ^' '^ ^ ' 

using that 

^P(i)^i.o(.-., 

We can now treat the terms containing the amplitude F{1) = F{2) = F{3) = F|. 
Exchanging the variables 9i ^ 62 in the second line and redefining 6^ — > 62 in the third 
line of eqn. f l5.24p results in 

2 y 27r 27r \ p(112) ^ ^ '^^ p(ll) ^ 

The combination of the various densities in this expression requires special care. From the 
large L asymptotics 

p{{)^EiL , p{ij)^E,EjL^ , p{ijk) EiEjEkL' , p(13, 2) ~ Ei^sL^ 

it naively scales with L. However, it can be easily verified that the coefficient of the leading 
term, which is linear in L, is exactly zero. Without this, the large L limit would not make 
sense, so this is rather reassuring. We can then calculate the sub-leading term, which 
requires tedious but elementary manipulations. The end result turns out to be extremely 
simple 

so the contribution in the L —>■ 00 limit turns out to be just 

2 J 2tx 2tx 

Summing up the contributions (I5.2ip . (15.221) . (15.25P and (15.271) we obtain 

1 f ^^^^s^^^ ^^-jg-mi?(coshei+coshe2+coshe3) 

6 J 271 271 271 ^ ' ' 

/ dOi d92 Q \ -m_R{cosh9i+2coshe2) i [ ^^^s -SmRcoshSi 

27r 27r ^' J 27t ^ 

_1 r dB^d^ps^^Q^ _ ^^^^-mRiM+2cosi.e,) ^5 28) 

2 J 271 271 

which agrees exactly with eqn. (15.101) . 

This result gives an independent support for the LeClair-Mussardo expansion. The 
calculation is model independent, and although we only worked it out to order e~^™^, 
it is expected that it coincides with the LeClair-Mussardo expansion to all orders. For a 
complete proof we need a better understanding of its structure, which is out of the scope 
of the present work. 
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5.3 Low-temperature expansion for the two-point func- 
tion 

The method presented in the previous section has a straightforward extension to higher 
point correlation functions. For example, a two-point correlation function 



^ TrL(e-«^^^) 



can be expanded inserting two complete sets of states 

TvL (e-^^^Oi(x,t)C>2(0)) = 5^e-^^"(^)(n|O(x,t)i|m)i(m|O2(0)|n)i (5.29) 

m,n 

The above expression can be evaluated along the lines presented in the previous section, 
provided that the intermediate state sums are properly truncated. The explicit evaluation 
of expansion f l5.29p has not yet been carried out. As a first step, here we investigate the 
simplest terms and calculate the first nontrivial contribution; the systematic evaluation of 
higher order terms is left for further work. We wish to remark, that the material presented 
in this section is unpublished. 
We define 

(Oi(x,t)O2(0))f = JiVM (5.30) 



Z 

N,M 



where 



Inm — 

/i.../jv Ji...Jm 

{{J, . . . Jm}|O2(0)|{/i . . . /^})^e^(^-^^)^e-^^(^-*)e-^^* (5.31) 

and £^1,2 and Pi, 2 are the total energies and momenta of the multi-particle states. 

It is easy to see, that the terms = 0, M = . . . 00 and = . . . 00, M = add up 

to zero-temperature correlation functions 

00 00 

J] J0A/= (Ol(x,t)O2(0))L 5^/^o= (Ol(x,i?-t)O2(0))L 

M=0 N=0 

The first nontrivial term is therefore In, which is given by 

hi = J2 ({n\Om\{J})L {{J}\02m{I})L e^(^-^^)^e-^^(^-*)e-^^' (5.32) 
I, J 

where now Ei^2 = mcosh(6'i 2) and pi^2 = msinh(6'i,2) are finite size one-particle energies 
and momenta. Upon f l2.18p and (12.321) the two-particle matrix elements are given by 

(wio,(o)i{j}>. = V''V'':!!i +^.Ao.) 

\/pi{0i)pi{62) 

V Pi (C*!) Pi (6*2) 
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Substituting the above formulas into (15.321) one obtains 

r = V ^2^^^^! +^^.^2)Ff-(g2 +^Vr,gi) 

^ ^ Pl(^2) ^ ^ Pl(^l) 

+ 5^(Oi)(02)e-^^« (5.33) 

The first term in (15.331) can be transformed in the L —>■ oo Umit into the well-defined 
(singularity free) double integral 

J 271 271 ' ' 

(5.34) 

In the second and third terms of (15.331) one can recognize the first thermal corrections to 
the expectation values of O2 and Oi, as given by (15.190 . The last term is of 0{L) and 
thus divergent; however, it gets canceled by the 0{L) term coming from the expansion 

1/Z = 1 + 0(6"""^) 

Putting everything together 
{0^{x,t)02m'' = {Oif W + 

+ {{0^ix,t)02m - (Oi) {02)) + {{0^{x,R-t)02m - {O,) (O2) ) (5.35) 



+ / ^^F^'iOi + ZTT, e2)F^\e2 + ZTT, ^^)e^(-nh(eO-sinh(e2))m.g-Mfl-*) cosh(ei)-mt coshCe^) 

27r 271 
+0(e-'"^) 



Equation (15. 35^ is a new result of this work and it can serve as a starting point to 
calculate higher order terms. § 

5.3.1 Comparison with the LeClair-Mussardo proposal 

Based on the TBA approach LeClair and Mussardo proposed the following formula 
for the thermal two-point function (equation (3.3) in [30] ) 

(O(x,t)O(0,0))^= {{Off + Y,^, E / ^•••^ \f{UAOs)e.^[-a,{te,+rxk,) 

x|(0|O|^^i...^^^)^^_j' (5.36) 



^We wish to note, that similar techniques are used in a recent work [116) which appeared after the 
completion of this thesis. In addition to the results above, the regularized finite contributions of /12 and 
/21 are also determined in ^116j . 
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where fa,{Oj) = 1/(1 + e-'^^'^(^^)), Sj = e{9j)/R and kj = k{9j) with e{9) being the solution 
of the TBA equations (I5.6P and k{9) given by 

k{9) = msmh{9) + j d9' 5(9-9') pi{9') 

2npi{9){l + e"(^)) = mcosh(^) + j d9'^{9 - 9')pi{9') 

The form factors appearing in ( 15.360 are defined by 

{Q\0\9i . . . 9n),^,„„^ = F^{9, - mai, ...,9^- ^vraTv) a, = (1 - a,)/2 E {0, 1} 

We do not attempt here a thorough analysis of (15.36^ : nevertheless we wish to point 
out, that to order 0(e~™^) it coincides with our result (15.35p . The double integral (15.340 
is given for example by the N = 2, ai 2 = (+1, —1) term in ( 15.360 after one substitutes 
the Oth-order approximations 

e{9) = mR cosh(^) + 0(6"'"^) k{9) = m sinh(e) + 0(6"'"^) 

The validity of the LeClair-Mussardo proposal was doubted by two different authors 
[3ll[35]. We did not find a discrepancy at order ^(e""^^); it is a subject of future research 
to decide whether all terms in our low-temperature expansion are properly reproduced by 
( 1536)1 . 
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Chapter 6 
Conclusions 



In this thesis we investigated finite size eff'ects in 1+1 dimensional integrable quantum 
field theories; we considered massive integrable models with diagonal scattering. There 
were two main subjects of interest: finite volume form factors of local operators and finite 
size effects to correlation functions. 

In the first main part of the work we gave a description of finite volume form factors 
in terms of the infinite volume form factors (solutions of the form factor bootstrap) and 
the S-matrix of the theory. 

First we compared the spectral representation of correlation functions in an infinite 
system and in a finite volume L. We showed that the matrix elements in finite volume 
essentially coincide with the infinite volume form factors up to a nontrivial normalization 
factor, which is related to the particle densities of the finite volume spectrum. We showed, 
that our results capture all finite size corrections which scale as powers of 1/L. We also 
showed that the residual finite size effects are of order e"'^^ with a universal characteristic 
mass scale of the theory, which does not depend on the particular form factor in question. 

As a second step, we conjectured a formula for generic matrix elements without dis- 
connected pieces. Together with the elementary case this is a new result of the author, 
although it can be considered as a generalization of partial results which appeared in 
[8] and independently in [52l I112| . Disconnected terms occur if there is a particle with a 
given rapidity, which present in both the "bra" and the "ket" vector. We showed that in 
finite volume this only happens in the case of diagonal form factors and matrix elements 
of "parity-symmetric" states including zero-momentum particles. 

The diagonal matrix elements require special care because of the presence of various 
disconnected terms and the ambiguity of the infinite volume diagonal form factors. We 
used form factor perturbation theory to derive the first two cases (diagonal matrix el- 
ements of one-particle and two-particle states) in terms of the symmetric evaluation of 
the infinite volume diagonal form factors. Based on this result, we conjectured a general 
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formula, which we expressed with both the symmetric and connected evaluations. The 
latter expression was used to prove that our result coincides with a conjecture made inde- 
pendently by Saleur in [31]. The rigorous proofs concerning the two simplest cases (eqs. 
f l2.32p and f l2.33p ) are new results of the author, together with the determination of the 
general relation between the symmetric and the connected form factors (Theorem 1). 

To complete the description of finite volume matrix elements we also conjectured 
a formula for the case of states with zero-momentum particles. The most general result 
concerns situations with possibly more than one zero-momentum particles present in both 
states; this only occurs in theories with multiple particle species. In |117) we published the 
result for theories with only one particle type; the generalization (12.54p is an unpublished 
result of this work. 

Following the analytic study, in section [3] we turned to numerical methods to confirm 
our results and conjectures. We investigated the massive Lee- Yang model and the critical 
Ising model in a magnetic field and considered form factors of the perturbing field, and 
the energy operator, respectively. We used TCSA as a numerical tool to obtain finite 
volume spectra and form factors. First we identified finite volume states at different values 
of the volume L: we matched the numerical data with predictions of the Bethe-Yang 
equations. We then calculated the finite volume form factors and compared them to 
our analytic formulas (the infinite volume form factors were already available in both 
models). In all cases (including form factors with disconnected pieces) we observed a 
satisfactory agreement in the scaling region, where both the residual finite size effects and 
the truncation errors are negligible. 

We stress that our numerical investigation not only confirms our results about ma- 
trix elements in finite volume, but it is also the first direct test of the infinite volume 
form factors. These functions are obtained by solving the form factor axioms of the boot- 
strap program and then selecting the solutions with the desired symmetry and scaling 
properties. Although the identification of scattering theories as perturbed conformal field 
theories is well-established, there were only a few direct tests of the individual form factor 
functions prior to our work; the usual tests in the literature proceed through evaluation 
of correlation functions, sum-rules, etc. In this work, on the other hand, we directly com- 
pared the form factor functions at different values of the rapidity parameters to TCSA 
data, thus providing evidence for the applicability of the bootstrap approach to form 
factors. 

In the case of the Ising model we observed that in some situations there is no scal- 
ing region, ie. the residual finite size corrections remain relevant even in large volumes, 
which are out of the reach of our TCSA routines. The large exponential corrections were 
explained by the presence of a /z-term with a small exponent. These findings served as 
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a motivation for section [H where we investigated the /i-terms associated to finite size 
energies and finite volume form factors. 

It is well-known that the /i-term is always connected with the inner structure of the 
particles, ie. their composition under the bootstrap procedure. We applied a simple quan- 
tum mechanical picture of bound state quantization in finite volume: the momenta of 
the constituents take complex values but they are still determined by the analytical con- 
tinuation of the Bethe-Yang equations. We showed that if this principle is implemented 
properly, the /i-term can always be obtained by analytically continuing the usual formulas 
to the complex rapidities of the constituents. Moreover it was shown, that in some circum- 
stances it is possible for the constituents of a bound state to unbind. This phenomenon 
was demonstrated numerically in case of the Ising model, where moving A3 states with 
odd momentum quantum numbers dissociate in small volume into conventional AiAi 
scattering states. 

In the second part of the work (section [5|) we investigated finite temperature corrections 
to vacuum expectation values and correlation functions. We introduced finite volume 
as a regulator of the otherwise ill-defined Boltzmann sum and established the spectral 
representation of n-point functions in terms of the finite volume form factors. We showed 
that it is possible to derive a systematic low temperature expansion with a well-defined 
L — s> 00 limit. 

In the case of the one-point function we calculated the first three nontrivial terms and 
found complete agreement with the LeClair-Mussardo approach. It was pointed out, that 
the Delfino proposal differs from the LeClair-Mussardo formula exactly at this order; the 
discrepancy is easily explained in light of our results. Any evaluation scheme of thermal 
averages based on a spectral representation has to deal with infinities, which are intro- 
duced by disconnected terms of the form factors and by different contributions to the 
partition function itself. These infinities cancel as expected, however a well-defined regu- 
larization procedure is needed to obtain the left-over finite pieces. This point was missed 
in the proposals prior to our work. The LeClair-Mussardo formula is believed to be correct 
to all orders, however a general proof (possibly based on a re-summation scheme of our 
low-temperature expansion) is not available. 

In our calculations the volume L served as a cutoff to regulate the divergent discon- 
nected contributions. Those terms in the low-temperature expansion which scale with 
positive powers of L add up to zero, as it is necessary to obtain a meaningful result. We 
stress that in order to obtain the correct C(l) terms it was crucial to use the interacting 
densities p^"'\ which differ from the free-theory densities by sub-leading terms. 

In this work we also considered the two-point function. We only calculated the first 
non-trivial term, which is a new (unpublished) result of this work. Our leading-order 
calculation showed agreement with the LeClair-Mussardo approach. A study of higher 
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order terms is left for future work, which will hopefully shed more light on the validity of 
the LeClair-Mussardo proposal. 

There are several directions to generalize our results on finite volume form factors. 

The extension to boundary field theories has already been established |118| together 
with the evaluation of finite temperature expectation values of boundary operators |119| . 

The generalization to non-diagonal scattering theories would be desirable, since these 
theories (for example the 0(3) a-model) serve as effective field theories describing long- 
range interactions is real- word condensed matter systems (Heisenberg spin chains, etc.). 
We expect that the principles laid out in this work (finite volume matrix elements and 
infinite volume form factors connected by a density factor) carry over to the non-diagonal 
case. However, non-diagonal scattering poses technical difficulties. One one hand, the 
quantization of scattering states in finite volume is more complicated, since it involves a 
diagonalization of the transfer matrix. On the other hand, infinite volume form factors 
(decomposed into the different channels of the scattering) are also more difficult to obtain. 

It is an interesting question whether some of our results can be extended to massless 
scattering theories. These models can be used to study renormalization group flows with 
a nontrivial conformal IR flxed point [75 l [76 l [T20] . Although the form factor bootstrap for 
massless theories has already been established |120l I121| . it is not clear what happens in a 
finite volume. In the absence of a mass-gap one does not have control over residual finite 
size effects, which decay exponentially in a massive theory. The situation is even worse 
in the case of thermal expectation values and correlation functions: there is no suitable 
variable (such as e~"^^ in a massive theory) to perform a low-temperature expansion. 

Finally we would like to remark that an exact description of finite volume form factors 
would be desirable, since they could be used to calculate correlations at finite temper- 
ature in a different way, than the one presented in this work. Instead of dealing with a 
theory defined in infinite volume and finite temporal extent R = 1/T one can perform 
an Euclidean rotation (r, x) (— x, r) to deal with a zero-temperature system defined 
in a finite volume R. Thermal corrections then become finite size corrections which can 
be readily calculated by a conventional zero-temperature form factor expansion. Obvi- 
ously one needs the exact finite volume form factors in this approach: the summation over 
form factors not including the exponential corrections only reproduces the infinite volume 
correlation functions. 

As far as we know, the only exact results available concern the form factors of the order 
and disorder operators in Ising model with zero magnetic field [43l SH |45]. However, 
these calculations are rather model-specific and it is not clear at all how to approach 
more complicated (interacting) models. A TBA-like integral equation similar to the one 
describing excited state energies [25] would be desirable; however, it is not clear how to 
solve this problem. 
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